Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 394,473 G→T 5.1% intergenic (+55/‑33) yaiT → / → yaiV pseudogene, autotransporter family;putative structure; Not classified; interrupted by IS3; putative flagellin structural protein/putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913394,4730GT 5.1% 409.4 / 7.2 135intergenic (+55/‑33)yaiT/yaiVpseudogene, autotransporter family;putative structure; Not classified; interrupted by IS3; putative flagellin structural protein/putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (74/53);  new base T (3/5);  total (77/58)
Fisher's exact test for biased strand distribution p-value = 2.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

Only 100 of 135 total aligned reads displayed.

ACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGATG  >  NC_000913/394327‑394620
                                                                                                                                                  |                                                                                                                                                   
aCTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGAtt                                                                                                                                                  <  1:1468125/150‑1 (MQ=255)
 cTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGAttt                                                                                                                                                 <  1:1005673/150‑1 (MQ=255)
 cTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGAttt                                                                                                                                                 <  1:822169/150‑1 (MQ=255)
     aGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTTGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACg                                                                                                                                             <  1:55512/150‑1 (MQ=255)
     aGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACg                                                                                                                                             >  1:1129038/1‑150 (MQ=255)
       cgcCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACggg                                                                                                                                           <  1:1579717/150‑1 (MQ=255)
        gacTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTTGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAATCCTTGACGGGc                                                                                                                                          <  2:691441/148‑1 (MQ=255)
         ccTATACCGATGCTAACTACCTCGGTGGGGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAGTAATGCCCGTCAAATCCTTGACGGGCa                                                                                                                                         <  2:319709/150‑1 (MQ=255)
          cTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCAt                                                                                                                                        <  1:427992/150‑1 (MQ=255)
          atataCGGATGCTAACTACCTCGGTTGTTGTGACGTAGATAAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTCGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAATCCTTGACGGGCAt                                                                                                                                        <  1:263328/149‑1 (MQ=255)
           tataCCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCAtt                                                                                                                                       >  1:1678397/1‑150 (MQ=255)
              aCCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCTAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGAAGGGCATTACt                                                                                                                                    >  2:1015544/1‑150 (MQ=255)
               ccGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTg                                                                                                                                   >  2:367219/1‑150 (MQ=255)
               ccGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTg                                                                                                                                   >  2:1111133/1‑150 (MQ=255)
                 gATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCa                                                                                                                                 >  2:1250236/1‑150 (MQ=255)
                  aTGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCaa                                                                                                                                >  1:1091975/1‑150 (MQ=255)
                  aTGCTAACTACCTCGGGGGGGGGGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCaa                                                                                                                                <  1:873177/150‑1 (MQ=255)
                    gCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAg                                                                                                                              >  2:77813/1‑150 (MQ=255)
                        aCTACCTCGTTGGTGGTGACGTAGAGCAAGACTGGTCCGCTAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGAcg                                                                                                                          <  1:1563276/150‑1 (MQ=255)
                          tACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGAcgcg                                                                                                                        <  2:1317374/150‑1 (MQ=255)
                           aCCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACgcgc                                                                                                                       >  1:378174/1‑150 (MQ=255)
                           aCCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTAGACGGGCATTACTGCAAAGGACgcgc                                                                                                                       <  2:128731/150‑1 (MQ=255)
                                ggtggtggtGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATg                                                                                                                  <  2:118668/150‑1 (MQ=255)
                                 gtggtggtgACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGt                                                                                                                 <  1:754241/150‑1 (MQ=255)
                                   ggtggtgACGTAGATCTAGACTTGTCTGCGAACGTGTGTGTTAAATATACCTTGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTTAATAAAGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTa                                                                                                               <  1:64616/150‑1 (MQ=255)
                                   ggtggtgACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTa                                                                                                               <  1:1292685/150‑1 (MQ=255)
                                   ggtggtgACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTa                                                                                                               >  2:855356/1‑150 (MQ=255)
                                      ggtgACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGTATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCt                                                                                                            <  1:1116978/150‑1 (MQ=255)
                                        tgACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGt                                                                                                          >  2:1597097/1‑150 (MQ=255)
                                         gACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACTAGCAAATGTTATCTGTa                                                                                                         >  1:111826/1‑150 (MQ=255)
                                          aCGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAg                                                                                                        <  1:855356/150‑1 (MQ=255)
                                           cGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGt                                                                                                       >  2:676556/1‑150 (MQ=255)
                                           cGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGt                                                                                                       >  2:835982/1‑150 (MQ=255)
                                              aGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTaa                                                                                                    >  1:1171885/1‑150 (MQ=255)
                                               taTCAATACTTGTCCGCGAACGTGGGTTTTAAATATACCTGGTAATATTCTTCACTCCGAATAAATACTGGTAATTTAATATAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTaaa                                                                                                   <  1:1340660/149‑1 (MQ=255)
                                               gATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTaaa                                                                                                   <  1:367219/150‑1 (MQ=255)
                                                   aaGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACcta                                                                                               >  2:1677787/1‑149 (MQ=255)
                                                         ggTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAgg                                                                                         >  1:1038277/1‑150 (MQ=255)
                                                         ggTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAgg                                                                                         <  1:714146/150‑1 (MQ=255)
                                                         ggTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAgg                                                                                         <  1:687668/150‑1 (MQ=255)
                                                         ggTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTAAgg                                                                                         >  1:1246748/1‑150 (MQ=255)
                                                          gTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGa                                                                                        >  2:655551/1‑150 (MQ=255)
                                                           tCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGaa                                                                                       >  1:1357612/1‑150 (MQ=255)
                                                              gcgAACGTGGATGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAAttt                                                                                    >  2:38369/1‑150 (MQ=255)
                                                                   cGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTaa                                                                               >  2:175492/1‑150 (MQ=255)
                                                                      gggTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGTTAATTTAATATAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCt                                                                            <  1:1277851/150‑1 (MQ=255)
                                                                          gtTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGAt                                                                        >  1:780168/1‑150 (MQ=255)
                                                                           ttAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATa                                                                       <  2:1044178/150‑1 (MQ=255)
                                                                            tAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATaa                                                                      >  2:821629/1‑150 (MQ=255)
                                                                            tAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATaa                                                                      >  1:1068063/1‑150 (MQ=255)
                                                                               atttaCCTGGTAATATTCTTCACTCGTAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATg                                                                   <  1:1091802/147‑1 (MQ=255)
                                                                               atatACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATg                                                                   <  1:582314/150‑1 (MQ=255)
                                                                                 ataCCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGtt                                                                 <  1:524701/150‑1 (MQ=255)
                                                                                  taCCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGttt                                                                >  2:644628/1‑150 (MQ=255)
                                                                                  taCCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGttt                                                                >  2:286751/1‑150 (MQ=255)
                                                                                      tGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTcc                                                            >  2:1603315/1‑150 (MQ=255)
                                                                                      tGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTcc                                                            <  2:222536/150‑1 (MQ=255)
                                                                                      tGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTcc                                                            >  2:1687725/1‑150 (MQ=255)
                                                                                      tGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGATCGATAAATGTTTGTcc                                                            >  2:52903/1‑150 (MQ=255)
                                                                                         aaaTATTCTTCACTCCGAAGAAATTCTGGTATTTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAAAGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGa                                                         <  1:457264/149‑1 (MQ=255)
                                                                                            tatTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATAc                                                      >  1:177627/1‑150 (MQ=255)
                                                                                              ttcttcACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACgg                                                    <  1:1186446/150‑1 (MQ=255)
                                                                                              ttcttcACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACgg                                                    >  2:1478402/1‑150 (MQ=255)
                                                                                                cttcACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGCCGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTa                                                  >  2:895990/1‑150 (MQ=255)
                                                                                                cttcACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTa                                                  >  2:114087/1‑150 (MQ=255)
                                                                                                   cACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTAcgc                                               >  1:783589/1‑150 (MQ=255)
                                                                                                   cACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTAcgc                                               >  1:790177/1‑150 (MQ=255)
                                                                                                     cTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTAcgcga                                             >  2:794950/1‑149 (MQ=255)
                                                                                                     cTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACgcgc                                             >  1:1644171/1‑150 (MQ=255)
                                                                                                     cTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACgcgc                                             >  2:1148829/1‑150 (MQ=255)
                                                                                                       ccGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGAAGCGCAAATGTTAACTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCtt                                           >  1:1498334/1‑150 (MQ=255)
                                                                                                            taaaTACTGTTAATTTAATCTAAATAATGCCCGTCAAGTATTTTACGGGCATTACTGCAAAGGACGCGCAAATGATATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGt                                      <  1:835982/149‑1 (MQ=255)
                                                                                                             aaaTACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGtt                                     <  2:1038277/150‑1 (MQ=255)
                                                                                                               aTACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTta                                   >  2:1489504/1‑150 (MQ=255)
                                                                                                               aTACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTta                                   >  2:107314/1‑150 (MQ=255)
                                                                                                                    ggTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCg                                         >  2:807159/1‑139 (MQ=255)
                                                                                                                    ggTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTaataat                              >  2:1611317/1‑150 (MQ=255)
                                                                                                                    ggTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTaataat                              >  2:443663/1‑150 (MQ=255)
                                                                                                                     gTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATg                             <  1:179372/150‑1 (MQ=255)
                                                                                                                     gTAATTTAATCTAAATAATGCCAGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATg                             <  2:1085954/150‑1 (MQ=255)
                                                                                                                            aaTCTATATTATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCa                      <  2:675715/150‑1 (MQ=255)
                                                                                                                             aTCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCaa                     >  1:1224472/1‑150 (MQ=255)
                                                                                                                              tCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAg                    <  1:644628/150‑1 (MQ=255)
                                                                                                                               cTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGt                   >  2:819993/1‑150 (MQ=255)
                                                                                                                                tAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGtt                  >  2:1048546/1‑150 (MQ=255)
                                                                                                                                tAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGt                   <  2:47994/149‑1 (MQ=255)
                                                                                                                                tAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGt                   >  1:47994/1‑149 (MQ=255)
                                                                                                                                   ataatGCCCGGCAAGGATTTTACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAAACTCTTCAGGAATTTGGTTAGCTCGATAAATGTTTGTCCAGTTACGGTACGCGCTTCGAGTTTAATTATGAAAAGCAAGTtata               <  1:1015544/150‑1 (MQ=255)
                                                                                                                                     aatGCCCGTCAAATCCTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATAtt             >  2:146407/1‑150 (MQ=255)
                                                                                                                                       tGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAg                    <  2:575534/141‑1 (MQ=255)
                                                                                                                                       tGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATAtttt           >  1:308827/1‑150 (MQ=255)
                                                                                                                                        gCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTc          >  1:1653272/1‑150 (MQ=255)
                                                                                                                                        gCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTc          >  1:1614890/1‑150 (MQ=255)
                                                                                                                                        gCCCGTCAAATCCTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTc          >  2:645434/1‑150 (MQ=255)
                                                                                                                                         cccGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTcc         >  1:479320/1‑150 (MQ=255)
                                                                                                                                         cccGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTcc         <  2:226946/150‑1 (MQ=255)
                                                                                                                                          ccGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCa        <  1:175492/150‑1 (MQ=255)
                                                                                                                                            gTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGt      >  1:284241/1‑150 (MQ=255)
                                                                                                                                             tCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTg     >  2:196460/1‑150 (MQ=255)
                                                                                                                                                aGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGATg  <  1:1148829/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGATG  >  NC_000913/394327‑394620

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