Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,001,9180AC42.4% 116.5 / 16.5 66intergenic (‑129/+192)fliA/fliCRNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base A (35/3);  new base C (0/28);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 1.32e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCA  >  NC_000913/2001770‑2002067
                                                                                                                                                      |                                                                                                                                                      
gCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATAGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑cc                                                                                                                                                     >  1:1224441/1‑151 (MQ=255)
  ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑cccc                                                                                                                                                   <  1:223204/151‑1 (MQ=255)
       agagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCgg                                                                                                                                              >  2:1246961/1‑151 (MQ=255)
                 cACGAAAAACAGCCCTGCGTTATATGATTTATCGGCATGATTACCCTTTTCTGCAGGTTTTTTAATCGGACGATTGGGGGGGAAAATG‑GGGGTTTTTT‑GGGGGTTTACAGGTAAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGaa                                                                                                                                    <  2:1265977/151‑1 (MQ=255)
                   caaaaaaCAGCCTTGCGTTATATGAGTTATCGGCATGATTATCCTTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTAAAATGAGGGGTTTTTG‑GGGGG‑TTACAGTAAAATTCCAGGAAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAg                                                                                                                                   <  2:200477/146‑1 (MQ=255)
                     aTAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCAc                                                                                                                                >  1:826583/1‑151 (MQ=255)
                       aaaCAGCCCTGCTTTATATGATTTATCGGCATGATTATCCTTTTCTGCAGGGTTTTAAATCGGACGTTTGGGGGGTGAAATGGGGGTTTTTTT‑GGGGG‑TTACAGGTAAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAATCACg                                                                                                                               <  1:1226505/150‑1 (MQ=255)
                         aCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTg                                                                                                                            >  1:772462/1‑151 (MQ=255)
                              ccTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGAc                                                                                                                       >  2:1081651/1‑151 (MQ=255)
                               cTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGAc                                                                                                                       >  2:440201/1‑150 (MQ=255)
                                tGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                                      >  2:1018132/1‑150 (MQ=255)
                                     tataTGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgc                                                                                                                >  1:937358/1‑151 (MQ=255)
                                      atatGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcg                                                                                                               >  2:1244446/1‑151 (MQ=255)
                                         tGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCtt                                                                                                            >  1:257137/1‑151 (MQ=255)
                                                 cGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGc                                                                                                    >  2:847005/1‑151 (MQ=255)
                                                     aTGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                                 >  1:811698/1‑150 (MQ=255)
                                                      tGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGTTTAGGGGGTAAAATGAGGGGTTTTTG‑GGGGG‑TTACAGGTAAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGa                                                                                                <  2:193271/150‑1 (MQ=255)
                                                        aTTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgc                                                                                             >  2:1264453/1‑151 (MQ=255)
                                                         ttATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGAAAAAGGGGGGGTTATTG‑GGGGG‑TTACAGATAAATTCCAGGCAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt                                                                                               <  2:839779/148‑1 (MQ=255)
                                                           aTCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcgg                                                                                           >  2:225128/1‑150 (MQ=255)
                                                             ccGTTTCTGCAGGTTTTTTACGGGAATGTTTAGTGGGTGAAAAGAGGGTTTATTG‑GGGGG‑TTCCAGCTAACTTCGAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt                                                                                               <  2:1024297/144‑1 (MQ=255)
                                                             ccGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAAGGGGGGTTTATTT‑GGGGG‑TTACAGGTAAATTCCGGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggc                                                                                          <  1:1081631/149‑1 (MQ=255)
                                                                tttCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaa                                                                                      >  2:296526/1‑150 (MQ=255)
                                                                  tCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACg                                                                                   >  2:1104157/1‑151 (MQ=255)
                                                                    tGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGcc                                                                                 >  2:1075743/1‑151 (MQ=255)
                                                                     gcgggttttttAATCGGACGATTGGTGGTTAAAATGAGGGGTTTTTT‑GGGGG‑TAAAAGGAAAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGcc                                                                                 <  2:980227/144‑1 (MQ=255)
                                                                        gggTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                             >  1:737524/1‑151 (MQ=255)
                                                                           tttttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                          >  1:559252/1‑151 (MQ=255)
                                                                              ttAATGGAACGATTAG‑GGGGGAAAAGGGGGGTTTTTTGGGGGG‑TTAAAGGAAAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                        <  2:738254/150‑1 (MQ=25)
                                                                               tAATCGGACGATTAGTGGGTGAAAGGAGGGGTTTTTG‑GGGGG‑TAAAAGGTAAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                       <  1:796218/150‑1 (MQ=38)
                                                                                aaTGGGACGATTAGTGGGTAAAAGGAGGGTTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                      <  2:375143/150‑1 (MQ=38)
                                                                                aaTCGGACGATTGGGGGGTAAAAGGGGGGTTTATTG‑GGGGT‑TTACAGGTAAATCCCAGGAAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                       <  1:468875/149‑1 (MQ=37)
                                                                                     gACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaa                                                                 >  1:382738/1‑150 (MQ=255)
                                                                                        gATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                             >  2:859924/1‑151 (MQ=255)
                                                                                          ttggggggTAAAATGGGGGTTTATTT‑GGGGT‑TTACAGGTAAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAt                                                              <  2:591377/143‑1 (MQ=37)
                                                                                          ttggTGGGTGAAATGGGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgt                                                            <  1:353682/147‑1 (MQ=255)
                                                                                           tAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                          >  1:935472/1‑151 (MQ=255)
                                                                                           tAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                          >  2:227046/1‑151 (MQ=255)
                                                                                             gggggTGAAAGGGGGGTTTTTTT‑GGGGT‑TTACAGGAAAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                          <  1:489237/147‑1 (MQ=37)
                                                                                             gTGGGTAAAATGGGGGTTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                                           <  2:713480/148‑1 (MQ=255)
                                                                                             gTGGGTAAAATG‑GGGGTTTTTT‑GGGGGTTTACAGGTAAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa                                                        <  2:466786/151‑1 (MQ=37)
                                                                                             gTGGGTAAAAGGGGGGGTTTTTT‑GGGGG‑TTACAGGAAAATTCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                          <  1:1015593/149‑1 (MQ=38)
                                                                                              tGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa                                                        >  1:490067/1‑150 (MQ=255)
                                                                                                 gtaaaaTGAGGGTTTATTG‑GGGGT‑TTAAAGG‑AAACTTCAGGCAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAtt                                                     <  2:414943/147‑1 (MQ=25)
                                                                                                 ggaaaaatgggggTTTTTT‑GGGGGTTTAAAGGAAAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAtt                                                     <  1:895777/142‑1 (MQ=25)
                                                                                                    aaagggggggTTATTG‑GGGGG‑TTACAGTTAAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa                                                   <  2:356118/143‑1 (MQ=37)
                                                                                                    aaaTGGGGGTTTATTT‑GGGGG‑TTAGAGATAAATCCAAGCCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa                                                   <  1:38363/149‑1 (MQ=37)
                                                                                                    aaaTGGGGGTTTATTG‑GGGGT‑TTAAAGGAAAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa                                                   <  2:990745/149‑1 (MQ=37)
                                                                                                        gAGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt                                             >  1:461394/1‑151 (MQ=255)
                                                                                                          ggggTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                           >  1:415299/1‑151 (MQ=255)
                                                                                                              tttttt‑GGGGT‑TTAAAGGAAAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                                          <  1:991684/145‑1 (MQ=25)
                                                                                                                 ttt‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                                     >  1:125465/1‑150 (MQ=255)
                                                                                                                 ttt‑GGGGG‑TTAAAGGTAAATCCCAGAAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                      <  1:1198551/149‑1 (MQ=37)
                                                                                                                 ttt‑GGGGG‑TTAAAGGTAAATCCAAGAAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                      <  1:878705/149‑1 (MQ=25)
                                                                                                                  tt‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                                    >  1:623557/1‑150 (MQ=255)
                                                                                                                   t‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAAc                                  >  2:339704/1‑151 (MQ=255)
                                                                                                                      gggg‑TTACAGTAAATTTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt                                <  2:273406/151‑1 (MQ=37)
                                                                                                                            tACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg                           >  1:1010380/1‑151 (MQ=255)
                                                                                                                             aaaGGTAAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGc                          <  1:1075424/149‑1 (MQ=255)
                                                                                                                             aCAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGc                          >  2:767891/1‑151 (MQ=255)
                                                                                                                                gaaaaaTCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGcc                         <  1:1018113/146‑1 (MQ=255)
                                                                                                                                  tAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt                     >  1:409560/1‑151 (MQ=255)
                                                                                                                                   aaaTTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt                     >  2:893510/1‑150 (MQ=255)
                                                                                                                                   aaaTTCCAGGCAAAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa                    <  2:281450/151‑1 (MQ=255)
                                                                                                                                     aTTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa                    >  2:1232330/1‑149 (MQ=255)
                                                                                                                                             cAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTACGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc          >  1:672725/1‑151 (MQ=255)
                                                                                                                                               aaaaaaac‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAt        <  2:598360/150‑1 (MQ=255)
                                                                                                                                               aaaaaaac‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAt        <  1:1264427/150‑1 (MQ=255)
                                                                                                                                                aaaaaaa‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATc       >  1:861757/1‑151 (MQ=255)
                                                                                                                                                    aaa‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAgac   >  2:388237/1‑149 (MQ=255)
                                                                                                                                                      a‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCa  >  1:1183276/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
GCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCA  >  NC_000913/2001770‑2002067

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