Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 743,260 Δ1 bp 5.6% coding (434/744 nt) ybgI → NIF3 family metal‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913743,2600G. 5.6% 242.3 / 7.6 71coding (434/744 nt)ybgINIF3 family metal‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (36/31);  new base . (2/2);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

GCATTTGCCGCTTGACGCACATCCTGAGCTGGGCAATAACGCACAGCTGGCGGCGTTACTGGGGATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTGGTGCACGGGCGGCGGGCAAAGTTTTATCGATAGCGCCGCGCGTTTTGGCGTGGATGCTTTTATTACTGGC  >  NC_000913/743114‑743405
                                                                                                                                                  |                                                                                                                                                 
gcATTTGCCGCTTGACGCACATCCTGAGCTGGGCAATAACGCACAGCTGGCGGCGTTACTGGGGATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTg                                                                                                                                                   <  1:112209/147‑1 (MQ=255)
gcATTTGCCGCTTGACGCACATCCTGAGCTGGGCAATAACGCACAGCTGGCGGCGTTACTGGGGATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGa                                                                                                                                                 >  8:96101/1‑149 (MQ=255)
 cATTTGCCGCTTGACGCACATCCTGAGCTGGGCAATAACGCACAGCTGGCGGCGTTACTGGGGATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGAt                                                                                                                                                <  2:178619/149‑1 (MQ=255)
 cATTTGCCGCTTGACGCACATCCTGAGCTGGGCAATAACGCACAGCTGGCGGCGTTACTGGGGATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGAt                                                                                                                                                <  2:77603/149‑1 (MQ=255)
      gCCGCTTGACGCACATCCTGAGCTGGGCAATAACGCACAGCTGGCGGCGTTACTGGGGATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGAATGaa                                                                                                                                            <  7:322669/148‑1 (MQ=255)
       ccGCTTGACGCACATCCTGAGCTGGGCAATAACGCACAGCTGGCGGCGTTACTGGGGATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAg                                                                                                                                           <  8:202103/148‑1 (MQ=255)
          cTTGACGCACATCCTGAGCTGGGCAATAACGCACAGCTGGCGGCGTTACTGGGGATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGcgcg                                                                                                                                       >  6:303621/1‑149 (MQ=255)
             gACGCACATCCTGAGCTGGGCAATAACGCACAGCTGGCGGCGTTACTGGGGATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCt                                                                                                                                    >  1:377435/1‑149 (MQ=255)
               cGCACATCCTGAGCTGGGCAATAACGCACAGCTGGCGGCGTTACTGGGGATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTgg                                                                                                                                  <  3:227161/149‑1 (MQ=255)
                gCACATAATGAGCTGGGCAATAACGCACAGCTGGCGGCGTTACTGGGGATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTgg                                                                                                                                  >  3:39477/1‑148 (MQ=255)
                  acaTCCTGAGCTGGGCAATAACGCACAGCTGGCGGCGTTACTGGGGATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCGGGATTGAAGCGCGTCTGGGAc                                                                                                                               >  5:42492/1‑149 (MQ=255)
                     tCCTGAGCTGGGCAATAACGCACAGCTGGCGGCGTTACTGGGGATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGt                                                                                                                             >  7:209616/1‑148 (MQ=255)
                      ccTGAGCTGGGCAATAACGCACAGCTGGCGGCGTTACTGGGGATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGt                                                                                                                             >  6:103810/1‑147 (MQ=255)
                          aGCTGGGCAATAACGCACAGCTGGCGGCGTTACTGGGGATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGcc                                                                                                                        <  5:307589/148‑1 (MQ=255)
                            cTGGGCAATAACGCACAGCTGGCGGCGTTACTGGGGATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGtt                                                                                                                     <  2:377435/149‑1 (MQ=255)
                               ggCAATAACGCACAGCTGGCGGCGTTACTGGGGATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATg                                                                                                                  <  5:103810/149‑1 (MQ=255)
                               ggCAATAACGCACAGCTGGCGGCGTTACTGGGGATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATg                                                                                                                  <  5:252600/149‑1 (MQ=255)
                               ggCAATAACGCACAGCTGGCGGCGTTACTGGGGATCAAGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACAATGCCTGTGCAGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATg                                                                                                                  >  5:263786/1‑149 (MQ=255)
                                        gCACAGCTGGCGGCGTTACTGGGGATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTgg                                                                                                                                  <  8:244445/124‑1 (MQ=255)
                                        gCACAGCTGGCGGCGTTACTGGGGATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTgg                                                                                                                                  >  7:244445/1‑124 (MQ=255)
                                            aGCTGGCGGCGTTACTGGGGATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATAcc                                                                                                     >  7:104094/1‑149 (MQ=255)
                                              cTGGCGGCGTTACTGGGGATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCg                                                                                                    >  1:379749/1‑148 (MQ=255)
                                               tGGCGGCGTTACTGGGGATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGa                                                                                                  >  2:318711/1‑149 (MQ=255)
                                               tGGCGGCGTTACTGGGGATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGa                                                                                                  >  2:190102/1‑149 (MQ=255)
                                                 gcggcgTTACTGGGGATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGAc                                                                                                 <  4:39477/148‑1 (MQ=255)
                                                      gTTACTGGGGATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAgg                                                                                           <  5:144403/149‑1 (MQ=255)
                                                           tGGGGATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTc                                                                                      >  2:391252/1‑149 (MQ=255)
                                                             gggATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTccc                                                                                    >  4:243916/1‑148 (MQ=255)
                                                             gggATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCa                                                                                    >  7:363854/1‑149 (MQ=255)
                                                             gggATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCAGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCa                                                                                    >  8:34311/1‑149 (MQ=255)
                                                              ggATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGgt                                                                                                                >  5:243822/1‑120 (MQ=255)
                                                              ggATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGgt                                                                                                                <  6:243822/120‑1 (MQ=255)
                                                               gATCACGGACATGGGCGAAATTGAGCCGTTGGTGCCGTGGGCCCAACTGCCCAAGCCTGAGCCGGTACTGGCGCTGGCATCCTGGATTTAAGCGAGTCAGGGACGTAACCCGTAATAGTGTGGCGATACCGCACCGGAGGTGGTCCa                                                                                    >  6:269623/1‑147 (MQ=255)
                                                                    cGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTc                                                                             <  1:99620/149‑1 (MQ=255)
                                                                        cATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCt                                                                         >  3:345954/1‑149 (MQ=255)
                                                                         aTGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTg                                                                        <  4:345954/149‑1 (MQ=255)
                                                                          tGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTgg                                                                       >  6:378104/1‑149 (MQ=255)
                                                                            ggCGAAATTGAGCCGTTGGTGCCGTGGGGCTAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTGGt                                                                      <  6:84323/148‑1 (MQ=255)
                                                                            ggCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTgg                                                                                                            <  2:255398/110‑1 (MQ=255)
                                                                            ggCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTgg                                                                                                            >  1:255398/1‑110 (MQ=255)
                                                                            ggCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTg                                                                        >  3:257791/1‑145 (MQ=255)
                                                                            ggCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTg                                                                        <  4:262729/145‑1 (MQ=255)
                                                                            ggCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTg                                                                        >  3:262729/1‑145 (MQ=255)
                                                                            ggCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTg                                                                        <  4:257791/145‑1 (MQ=255)
                                                                               gAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTGGTGCAc                                                                  <  1:190102/149‑1 (MQ=255)
                                                                                   ttGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTGGTGCACGggc                                                              >  7:362667/1‑149 (MQ=255)
                                                                                        ccGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGGCTGGGACGTAAGCCGTTATGGGGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTGGTGCACGGGcggcgg                                                         <  7:96101/149‑1 (MQ=255)
                                                                                        ccGTTGCTGCGGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTGTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTGGTGCACGGGcggcgg                                                         <  4:46679/149‑1 (MQ=255)
                                                                                                gCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTGGTGCACGGGCGGCGGGCAAAGtt                                                 >  6:220880/1‑149 (MQ=255)
                                                                                                    tGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTGGTGCACGGGCGGCGGGCAAAGTTTTAt                                             <  2:379749/149‑1 (MQ=255)
                                                                                                        gCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTGGTGCACGGGCGGCGGGCAAAGTTTTATCGa                                          >  2:118023/1‑148 (MQ=255)
                                                                                                            aCTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTGGTGCACGGGCGGCGGGCAAAGTTTTATCGATAgcg                                     >  7:153654/1‑149 (MQ=255)
                                                                                                             cTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTGGTGCACGGGCGGCGGGCAAAGTTTTATCGATAgcgc                                    <  5:378104/149‑1 (MQ=255)
                                                                                                               gACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGGCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTGGTGCACGGGCGGCGGGCAAAGTTTTATCGATAGcgccg                                  <  8:362667/149‑1 (MQ=255)
                                                                                                                aCCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTGGTGCACGGGCGGCGGGCAAAGTTTTATCGATAGcgccgc                                 >  8:383817/1‑149 (MQ=255)
                                                                                                                aCCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGCCGATACCGGACCTGAGGTGGTCCGGCGCGTCGCCTGGTGCACGGGCGGCGGGCAAAGTTTTATCGATAGcgccgc                                 >  3:7685/1‑149 (MQ=255)
                                                                                                                  cATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTGGTGCACGGGCGGCGGGCAAAGTTTTATCGATAGCGCcgcgc                               >  8:367599/1‑149 (MQ=255)
                                                                                                                          tgCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTGGTGCACGGGCGGCGGGCAAAGTTTTATCGATAGCGCCGCGCGTTTTg                         <  8:153654/147‑1 (MQ=255)
                                                                                                                             cGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTGGTGCACGGGCGGCGGGCAAAGTTTTATCGATAGCGCCGCGCGTTTTGGCGTg                    <  7:383817/149‑1 (MQ=255)
                                                                                                                             cGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTGGTGCACGGGCGGCGGGCAAAGTTTTATCGATAGCGCCGCGCGTTTTGGCGTg                    >  7:64477/1‑149 (MQ=255)
                                                                                                                              gggACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTGGTGCACGGGCGGCGGGCAAAGTTTTATCGATAGCGCCGCGCGTTTTGGCGTg                    >  2:98047/1‑148 (MQ=255)
                                                                                                                                   tGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTGGTGCACGGGCGGCGGGCAAAGTTTTATCGATAGCGCCGCGCGTTTTGGCGt                     >  5:36744/1‑142 (MQ=255)
                                                                                                                                   tGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTGGTGCACGGGCGGCGGGCAAAGTTTTATCGATAGCGCCGCGCGTTTTGGCGt                     <  6:36744/142‑1 (MQ=255)
                                                                                                                                      aGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTGGTGCACGGGCGGCGGGCAAAGTTTTATCGATAGCGCCGCGCGTTTTGGCGTGGATGCttt            >  1:47097/1‑148 (MQ=255)
                                                                                                                                        cTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTGGTGCACGGGCGGCGGGCAAAGTTTTATCGATAGCGCCGCGCGTTTTGGCGTGGATGCTTttat         <  1:98047/149‑1 (MQ=255)
                                                                                                                                         tGGCTTCCTGGATTGAAGCGCGTCTGGGAAGTAAGCAGTTATGGTGTGGCGATACCGGAGGTGAGGTGGTCCAGCGCGTCGCCTGGTGCACGGTCGGCGGGCAAAGTTTTATCGATAGCGCCGACAGTTTTGGCGTGGATGCtttttt         >  5:273103/1‑146 (MQ=255)
                                                                                                                                            cTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTGGTGCACGGGCGGCGGGCAAAGTTTTATCGATAGCGCCGCGCGTTTTGGCGTGGATGCTTTTATTACt     <  7:34311/149‑1 (MQ=255)
                                                                                                                                               ccTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTGGTGCACGGGCGGCGGGCAAAGTTTTATCGATAGCGCCGCGCGTTTTGGCGTGGATGCTTTTATTACTg    <  7:367599/147‑1 (MQ=255)
                                                                                                                                               ccTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGGGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTGGTGCACGGGCGGCGGGCAAAGTTTTATCGATAGCGCCGCGCGTTTTGGCGTGGATGCTTTTATTACTgg   <  8:209616/148‑1 (MQ=255)
                                                                                                                                                cTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCATCGCCTGGTGCACGGGCGGCGGGCAAAGTTTTATCGATAGCGCCGCGCGTTTTGGCGTGGATGCTTTTATTACTGGc  >  5:318422/1‑148 (MQ=255)
                                                                                                                                                 tGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTGGTGCACGGGCGGCGGGCAAAGTTTTATCGATAGCGCCGCGCGTTTTGGCGTGGATGCtt             >  7:48964/1‑136 (MQ=255)
                                                                                                                                                 tGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTGGTGCACGGGCGGCGGGCAAAGTTTTATCGATAGCGCCGCGCGTTTTGGCGTGGATGCtt             <  8:48964/136‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCATTTGCCGCTTGACGCACATCCTGAGCTGGGCAATAACGCACAGCTGGCGGCGTTACTGGGGATCACGGTCATGGGCGAAATTGAGCCGTTGGTGCCGTGGGGCGAACTGACCATGCCTGTGCCGGGACTGGAGCTGGCTTCCTGGATTGAAGCGCGTCTGGGACGTAAGCCGTTATGGTGTGGCGATACCGGACCTGAGGTGGTCCAGCGCGTCGCCTGGTGCACGGGCGGCGGGCAAAGTTTTATCGATAGCGCCGCGCGTTTTGGCGTGGATGCTTTTATTACTGGC  >  NC_000913/743114‑743405

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: