Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,512,306 G→C 10.5% I13M (ATC→ATG yhfZ ← putative DNA‑binding transcriptional regulator YhfZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,512,3060GC10.5% 335.7 / 27.2 114I13M (ATC→ATGyhfZputative DNA‑binding transcriptional regulator YhfZ
Reads supporting (aligned to +/- strand):  ref base G (52/50);  new base C (4/8);  total (56/58)
Fisher's exact test for biased strand distribution p-value = 3.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-01

Only 100 of 121 total aligned reads displayed.

CGTTCAATCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAA  >  NC_000913/3512157‑3512454
                                                                                                                                                     |                                                                                                                                                    
ggTTCAATCCGTATCGCTCCGCTTGATTCGAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGtg                                                                                                                                                      <  1:1485802/149‑1 (MQ=255)
 gTTCAATCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCCGCGTCGtga                                                                                                                                                     >  1:201858/1‑150 (MQ=255)
 gTTCAATCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGtga                                                                                                                                                     <  2:83549/150‑1 (MQ=255)
 gTTCAATCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGtga                                                                                                                                                     >  1:1106647/1‑150 (MQ=255)
      aTCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATgacg                                                                                                                                                >  1:1519234/1‑150 (MQ=255)
         cGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGacgacg                                                                                                                                             <  2:1508867/150‑1 (MQ=255)
          gTACCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGc                                                                                                                                            <  2:846658/150‑1 (MQ=255)
           tATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGcc                                                                                                                                           <  1:508222/150‑1 (MQ=255)
           tATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTCATCACGACGct                                                                                                                                           <  2:1113540/150‑2 (MQ=255)
              cGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAATGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTCATCACGACGctggc                                                                                                                                        <  1:377220/150‑5 (MQ=255)
              cGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTCATCACGACGctggc                                                                                                                                        <  1:306476/150‑5 (MQ=255)
              cGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTCATCACGACGctggc                                                                                                                                        <  2:712938/150‑5 (MQ=255)
               gCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACGGGACTTTTCACACAAGTTATTAAGGGCACGCGTGGTGATGACGACGCcttct                                                                                                                                       >  2:688288/1‑150 (MQ=255)
               gCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCcttct                                                                                                                                       >  1:1164607/1‑150 (MQ=255)
                cTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTCATCACGACGctggccc                                                                                                                                      <  1:421935/150‑7 (MQ=255)
                   cGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCttttt                                                                                                                                   >  1:901892/1‑150 (MQ=255)
                     cTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTg                                                                                                                                 <  2:350020/150‑1 (MQ=255)
                     cTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTg                                                                                                                                 <  2:626753/150‑1 (MQ=255)
                       tGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGAt                                                                                                                               <  1:887608/150‑1 (MQ=255)
                          ttCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATaaa                                                                                                                            >  1:1120340/1‑150 (MQ=255)
                            ccAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACg                                                                                                                          <  1:296287/150‑1 (MQ=255)
                             cAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGt                                                                                                                         >  2:1453938/1‑150 (MQ=255)
                                 gTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGa                                                                                                                     >  1:1611154/1‑150 (MQ=255)
                                 gTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGa                                                                                                                     <  1:263433/150‑1 (MQ=255)
                                  ttttCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGAc                                                                                                                    <  1:249602/150‑1 (MQ=255)
                                   tttCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGAATTCATTTGCCAGCTCATCTATGTTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACg                                                                                                                   >  2:1566992/1‑150 (MQ=255)
                                    ttCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGCCGCCTTCTTTTTTGATAAACGTTCGACGc                                                                                                                  >  1:1373449/1‑150 (MQ=255)
                                      cAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGCCGACGCCTTCTTTTTTGATAAACGTTCGAGGCat                                                                                                                >  2:1403793/1‑150 (MQ=255)
                                       aaCGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGGatt                                                                                                               >  2:1069598/1‑149 (MQ=255)
                                            ggCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATaaa                                                                                                          <  1:1586060/150‑1 (MQ=255)
                                                tgacgcAGGCCAACGGATGAACGGCATTAATTTGCCAGCTCATCTATGGTTTTCAATCCATTAACGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGCCGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAAtttt                                                                                                      >  2:816083/6‑150 (MQ=255)
                                                 gCGTCAGGCCAACGGATGAACTGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTc                                                                                                     <  2:853175/150‑1 (MQ=255)
                                                 gCGTCAGGCCAACGGATGAACGGCTTTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTc                                                                                                     <  1:633306/150‑1 (MQ=255)
                                                 gCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTc                                                                                                     <  1:633303/150‑1 (MQ=255)
                                                 gCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTc                                                                                                     <  2:901892/150‑1 (MQ=255)
                                                   gTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAg                                                                                                   <  1:126732/150‑1 (MQ=255)
                                                   gTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAg                                                                                                   <  2:1519234/150‑1 (MQ=255)
                                                   gTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAg                                                                                                   <  1:509917/150‑1 (MQ=255)
                                                     cAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTa                                                                                                 >  2:1394442/1‑150 (MQ=255)
                                                       ggCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTaaa                                                                                               >  1:1117414/1‑150 (MQ=255)
                                                        gCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAAt                                                                                              <  2:87504/150‑1 (MQ=255)
                                                         ccAACGGATTAACGGCATTCATTTGCCATCTCATCTATGGTTTTCAATCTATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAAtt                                                                                             <  1:688288/150‑1 (MQ=255)
                                                         ccAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAAtt                                                                                             >  1:559930/1‑150 (MQ=255)
                                                         ccAACCCATCAACGCCATTCATTTCCCACCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAAtt                                                                                             <  1:705774/150‑1 (MQ=255)
                                                          cAACGGATGAACGGCATTCATTTGTCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTg                                                                                            >  1:309727/1‑150 (MQ=255)
                                                             cGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAat                                                                                         <  1:379935/150‑1 (MQ=255)
                                                                 tGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATAttt                                                                                     >  1:197693/1‑150 (MQ=255)
                                                                 tGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATAttt                                                                                     <  2:452855/150‑1 (MQ=255)
                                                                   aaCGGCACTCATTTGCCAGCTCATCTATGGTTTTCAATCGACTACCGCACTTTTCACCCAACACATAACCGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGGTCGACGCATAATAAATTTTCAGTAAATTGAATATTtat                                                                                   <  2:1488181/150‑1 (MQ=255)
                                                                    aCGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTtata                                                                                  >  2:576955/1‑150 (MQ=255)
                                                                    aCGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTtata                                                                                  <  2:802484/150‑1 (MQ=255)
                                                                     cGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTatat                                                                                 >  1:1289883/1‑150 (MQ=255)
                                                                        cattcattTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCtt                                                                              <  2:520983/150‑1 (MQ=255)
                                                                         attcattTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTc                                                                             <  1:1512855/150‑1 (MQ=255)
                                                                         attcattTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTc                                                                             <  2:655087/150‑1 (MQ=255)
                                                                          ttcattTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCa                                                                            >  2:289128/1‑150 (MQ=255)
                                                                           tcattTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAg                                                                           <  2:1141363/150‑1 (MQ=255)
                                                                           tcattTGCCAGCTCATCTATGCTTTTCAATCGATTACCGCACTTTTCAACCAACAAATACCGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTTAATTGAATATTTATATCTTCAg                                                                           <  2:1599147/150‑1 (MQ=255)
                                                                                 gCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAAttt                                                                     >  2:913055/1‑150 (MQ=255)
                                                                                  ccAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTg                                                                    <  2:1530925/150‑1 (MQ=255)
                                                                                   cAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGa                                                                   >  1:1447087/1‑150 (MQ=255)
                                                                                   cAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGa                                                                   >  1:1307772/1‑150 (MQ=255)
                                                                                     gCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAg                                                                 <  1:1394442/150‑1 (MQ=255)
                                                                                      cTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGa                                                                <  1:1169105/150‑1 (MQ=255)
                                                                                           cTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAgg                                                           >  2:952622/1‑150 (MQ=255)
                                                                                            tATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGgc                                                          >  1:1499213/1‑150 (MQ=255)
                                                                                             aTGGTTTTCAATCGAGTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGATCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATATGAAGATAAGgcg                                                         >  2:229810/1‑150 (MQ=255)
                                                                                               ggTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGa                                                                         >  2:1631671/1‑132 (MQ=255)
                                                                                               ggTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGa                                                                         >  2:327029/1‑132 (MQ=255)
                                                                                               ggTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGa                                                                         <  1:327029/132‑1 (MQ=255)
                                                                                               ggTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCa                                                       <  2:1120340/150‑1 (MQ=255)
                                                                                                  tttCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACa                                                    >  1:52325/1‑150 (MQ=255)
                                                                                                     cAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAAc                                                     <  2:15485/146‑1 (MQ=255)
                                                                                                     cAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAAc                                                     >  1:15485/1‑146 (MQ=255)
                                                                                                     cAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCg                                                 >  2:1356163/1‑150 (MQ=255)
                                                                                                            ttACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGAtgtg                                          <  2:458622/150‑1 (MQ=255)
                                                                                                             tACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTgg                                         >  1:1453837/1‑150 (MQ=255)
                                                                                                             tACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTgg                                         >  1:187452/1‑150 (MQ=255)
                                                                                                               ccGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGaa                                       >  1:1276368/1‑150 (MQ=255)
                                                                                                                 gCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAAc                                     >  1:869833/1‑150 (MQ=255)
                                                                                                                     ttttCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCg                                 <  2:1424111/150‑1 (MQ=255)
                                                                                                                      tttCACCCAACAAATAACGGGCCCGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGa                                >  2:1180680/1‑150 (MQ=255)
                                                                                                                          aCCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCa                            >  1:825686/1‑150 (MQ=255)
                                                                                                                            ccAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACTGCGAGGCACt                          >  1:284579/1‑150 (MQ=255)
                                                                                                                              aacaaATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTc                        <  1:987579/150‑1 (MQ=255)
                                                                                                                                   aTAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACAtt                   <  2:309727/150‑1 (MQ=255)
                                                                                                                                   aTAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTAAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACAtt                   >  2:686403/1‑150 (MQ=255)
                                                                                                                                     aaCGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACAtttt                 <  1:133365/150‑1 (MQ=255)
                                                                                                                                      aCGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACAttttt                >  2:555495/1‑150 (MQ=255)
                                                                                                                                         agaaggCGTCGTCATCACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCg             >  1:98273/6‑150 (MQ=255)
                                                                                                                                            aggCGTCGTCATCACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGat          >  2:370634/3‑150 (MQ=255)
                                                                                                                                            aggCGTCGTCATCACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGat          >  1:983428/3‑150 (MQ=255)
                                                                                                                                             aGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGatt         <  2:1135754/150‑1 (MQ=255)
                                                                                                                                             aGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGatt         <  2:1499213/150‑1 (MQ=255)
                                                                                                                                               cgtcgtGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGattat       >  2:51982/1‑150 (MQ=255)
                                                                                                                                                gtcgtGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGattatt      <  1:1453938/150‑1 (MQ=255)
                                                                                                                                                  cttgatgaCTACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTTAATATTTATATCTTCAGGAATTTGAAGATAAGGCTCCACAGCGAGATGTGGCAACGGCGAGGCACTTCACATTTTTTCGGATTATTGc    <  2:706507/148‑1 (MQ=255)
                                                                                                                                                   gtGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCa   >  1:831274/1‑150 (MQ=255)
                                                                                                                                                   gtGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCa   >  1:1601273/1‑150 (MQ=255)
                                                                                                                                                    tgatgaCGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCaa  <  1:1566992/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGTTCAATCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAA  >  NC_000913/3512157‑3512454

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