Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,542,503 G→A 100% W699* (TGG→TGA feoB → ferrous iron transporter protein B and GTP‑binding protein; membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,542,5030GA97.7% 594.7 / ‑0.4 186W699* (TGG→TGAfeoBferrous iron transporter protein B and GTP‑binding protein; membrane protein
Reads supporting (aligned to +/- strand):  ref base G (3/3);  new base A (99/80);  total (103/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 190 total aligned reads displayed.

ACCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTT  >  NC_000913/3542406‑3542599
                                                                                                 |                                                                                                
aCCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTgagg                                                                                                <  1:3006163/100‑1 (MQ=255)
 ccTGGTTTTCGTCCTGCTGTATGCACTATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTgagg                                                                                                <  2:896157/99‑1 (MQ=255)
  cTGATTTTCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGc                                                                                              >  1:2185582/1‑100 (MQ=255)
   tGATTTTCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCt                                                                                             >  1:3488422/1‑100 (MQ=255)
   tGATTTTCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCt                                                                                             >  1:1647781/1‑100 (MQ=255)
   tGATTTTCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCAGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCt                                                                                             >  2:1735295/1‑100 (MQ=255)
    gATTTTCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTg                                                                                            <  1:381722/100‑1 (MQ=255)
       tttCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAat                                                                                         <  1:868488/100‑1 (MQ=255)
        ttCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAata                                                                                        <  1:665134/100‑1 (MQ=255)
        ttCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAata                                                                                        <  2:1917090/100‑1 (MQ=255)
        ntCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGCGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAata                                                                                        >  2:1129052/2‑100 (MQ=255)
         tCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAatat                                                                                       >  2:1563077/1‑100 (MQ=255)
          cGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATc                                                                                      >  2:655491/1‑100 (MQ=255)
             ccTGCTGTATGTACCATGTATCTCGGTGAGCGGGGCAATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGc                                                                                    <  2:2567082/99‑1 (MQ=255)
              ctgctgTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCtt                                                                                  >  1:155393/1‑100 (MQ=255)
              ctgctgTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCt                                                                                   <  2:415929/99‑1 (MQ=255)
                gctgTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTAc                                                                                <  1:1991431/100‑1 (MQ=255)
                gctgTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTAc                                                                                >  1:102657/1‑100 (MQ=255)
                  tgtttgtaCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTAc                                                                                <  2:536451/98‑1 (MQ=255)
                  tgtatgtaCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTc                                                                              >  1:3361574/1‑100 (MQ=255)
                   gtatgtaCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCa                                                                             >  1:85920/1‑100 (MQ=255)
                     atgtaCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACt                                                                           >  1:1847914/1‑100 (MQ=255)
                      tgtaCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTg                                                                          <  2:2972270/100‑1 (MQ=255)
                       gtgccAGGTATCTGGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTgg                                                                         <  2:1255515/97‑1 (MQ=255)
                        taCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGc                                                                        <  1:2541561/100‑1 (MQ=255)
                        taCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGc                                                                        >  2:3704649/1‑100 (MQ=255)
                          ccATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCaa                                                                      >  2:335882/1‑100 (MQ=255)
                          ccATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCaa                                                                      >  1:3345372/1‑100 (MQ=255)
                           cATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAAc                                                                     <  1:2264415/100‑1 (MQ=255)
                           cATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAAc                                                                     >  1:3590298/1‑100 (MQ=255)
                           cATGTATCTCGGTGATGGGCGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGGGAGAGCTGAATATCGCTTCCTCACTGGc                                                                        >  2:735108/1‑97 (MQ=255)
                            aTGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACa                                                                    <  1:3326292/100‑1 (MQ=255)
                            aTGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACa                                                                    <  1:818588/100‑1 (MQ=255)
                             ngTATCTCGGTGATGGGGGCTATAGCCCGTGAATCAAGCCGAGGCTGGATGAGCTTCTCAAACATGTCAGGGCTGAATCTCGCTTCCTCACTCGc                                                                        >  2:2080088/2‑95 (MQ=255)
                                 tctcGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGtt                                                               >  2:2060652/1‑100 (MQ=255)
                                  ctcGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGATGGCTGAATATCGCTTACTCACTGGCAACATTGTTc                                                              <  2:3139287/100‑1 (MQ=255)
                                  ctcGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTc                                                              <  1:3123689/100‑1 (MQ=255)
                                  ctcGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTc                                                              <  1:981951/100‑1 (MQ=255)
                                  ctcGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTc                                                              <  1:2808643/100‑1 (MQ=255)
                                  ctcGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTc                                                              >  1:3643686/1‑100 (MQ=255)
                                   tcGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCt                                                             >  1:3183127/1‑100 (MQ=255)
                                   tcGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCt                                                             >  1:1354671/1‑100 (MQ=255)
                                    cGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTa                                                            >  1:2243520/1‑100 (MQ=255)
                                    cGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTa                                                            >  1:100621/1‑100 (MQ=255)
                                    cGGTGATGGGGGCTATCGCCCGTGAATCAAACCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTa                                                            >  1:2366612/1‑100 (MQ=255)
                                    cGGTGACGGGGTCTATCGCCCGTCGATCAAGCCGTGGCTGGATAGGCTCCGCCATCCTGTGAGGTCTGAATATCGCTTACTCACTGGCAACCTTGTTCt                                                             <  2:2621757/99‑1 (MQ=255)
                                      gTGATGGGGGCTATCGCCCGTGAATGAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGCGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATc                                                          <  2:2221751/100‑1 (MQ=255)
                                       tGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCa                                                         <  1:2697323/100‑1 (MQ=255)
                                       tGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCa                                                         <  1:1449754/100‑1 (MQ=255)
                                       tGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCa                                                         >  1:3002140/1‑100 (MQ=255)
                                       tGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCa                                                         <  1:1806616/100‑1 (MQ=255)
                                       tGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCa                                                         <  2:3160324/100‑1 (MQ=255)
                                          tGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGt                                                      >  2:246312/1‑100 (MQ=255)
                                          cacgggCGATAGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAgg                                                       >  1:3540326/4‑99 (MQ=255)
                                            ggggCTATCGCCCATGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACCGGCCACATTGTTCTATCAGGTc                                                     >  2:2365788/1‑99 (MQ=255)
                                             gggCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGc                                                   >  1:2490714/1‑100 (MQ=255)
                                             gggCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGc                                                   >  2:2427805/1‑100 (MQ=255)
                                             gggCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGc                                                   >  2:624836/1‑100 (MQ=255)
                                             gggCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGc                                                   >  2:1835293/1‑100 (MQ=255)
                                             gggCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGc                                                   >  2:1990285/1‑100 (MQ=255)
                                              ggCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGcc                                                  >  1:164875/1‑100 (MQ=255)
                                                 tATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGc                                               <  1:827895/100‑1 (MQ=255)
                                                  aTCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCt                                              <  1:488750/100‑1 (MQ=255)
                                                     gCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACa                                           <  2:2791256/100‑1 (MQ=255)
                                                      cccGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAg                                          >  1:2157868/1‑100 (MQ=255)
                                                       ccGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGt                                         >  2:2926786/1‑100 (MQ=255)
                                                       ccGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGt                                         >  2:504463/1‑100 (MQ=255)
                                                         gTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTTTGAAGGCTGAATATCGCTTACGCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAg                                          <  2:3610871/97‑1 (MQ=255)
                                                         gTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCa                                       <  1:2129631/100‑1 (MQ=255)
                                                            aaTCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCa                                    <  1:69753/100‑1 (MQ=255)
                                                             aTCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGGCGCCAGCTACAGTCAGCAt                                   >  1:2293174/1‑100 (MQ=255)
                                                             aTCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCAt                                   >  2:715582/1‑100 (MQ=255)
                                                              tCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATc                                  >  1:928064/1‑100 (MQ=255)
                                                              tCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATc                                  >  1:3431906/1‑100 (MQ=255)
                                                                  gCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAAc                              >  1:1629090/1‑100 (MQ=255)
                                                                   ccGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACt                             <  1:100126/100‑1 (MQ=255)
                                                                     gTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTa                           >  1:611441/1‑100 (MQ=255)
                                                                       ggCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTAcn                         <  2:3361548/100‑2 (MQ=255)
                                                                        gCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACat                        <  2:345797/100‑2 (MQ=255)
                                                                          tGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGcc                      >  2:480058/1‑100 (MQ=255)
                                                                          tGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGc                       <  2:1181838/99‑1 (MQ=255)
                                                                             aTGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTgg                   >  1:1598533/1‑100 (MQ=255)
                                                                               gggCTTCTCCATCCTGTGACGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGgtg                 <  2:3237717/100‑1 (MQ=255)
                                                                               gatcTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGgtg                 <  2:813158/97‑1 (MQ=255)
                                                                               gatcTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGgtg                 <  2:2432555/97‑1 (MQ=255)
                                                                                ggCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGgtgt                >  1:2610552/1‑100 (MQ=255)
                                                                                ggCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGgtg                 >  2:2156939/1‑99 (MQ=255)
                                                                                 gCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGtgtg               >  2:3567323/1‑100 (MQ=255)
                                                                                 gCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGtgtg               >  1:1726445/1‑100 (MQ=255)
                                                                                    tctcCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCAt            <  1:1503770/100‑1 (MQ=255)
                                                                                    tctcCACCATGTGGGGGCTGGGTATCTCTTACTCACTGGCAACAATGTTCTTTCAGGTGGCCAGCTATAGTCAGCATCCAACTTACTGCCTGGTGTGCa             <  2:799071/99‑1 (MQ=255)
                                                                                      tcCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTc          <  2:3442851/100‑1 (MQ=255)
                                                                                       ccATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCt         >  2:1112734/1‑100 (MQ=255)
                                                                                       ccATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCt         >  1:2923098/1‑100 (MQ=255)
                                                                                         aTCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTgg       <  2:560708/100‑1 (MQ=255)
                                                                                         aTCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTgg       >  1:1571087/1‑100 (MQ=255)
                                                                                          tCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTGGc      <  2:3488395/100‑1 (MQ=255)
                                                                                           ccTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTGGCg     >  1:20629/1‑100 (MQ=255)
                                                                                            cTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTGGCgg    >  1:3618263/1‑100 (MQ=255)
                                                                                              gtgAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTGGCGGtt  >  2:1437873/1‑100 (MQ=255)
                                                                                                 |                                                                                                
ACCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTT  >  NC_000913/3542406‑3542599

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