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

AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCA  >  NC_000913/2001769‑2002063
                                                                                                                                                      |                                                                                                                                                 
aGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGTTTATCCGTTTCTGCAGGGGTTTTAATCGGACCATTAGGGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAAc                                                                                                                                                  >  1:1049547/1‑151 (MQ=255)
   ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAAcccc                                                                                                                                               <  2:926588/151‑1 (MQ=255)
   ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAAcc                                                                                                                                                 <  2:577568/149‑1 (MQ=255)
          agTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTg                                                                                                                                        >  2:50529/1‑151 (MQ=255)
                ttCACGAAAAAAAGCCCTGCGTTATATGATTTATCGGCATGATTATCCGTTTCTGCGGGGTTTTAAACCGGACGATTAGGGGGTAAAAGGGGGGGTTTTTTGGGGG‑TTACAGGAAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCggg                                                                                                                                   <  1:1145055/150‑1 (MQ=255)
                  cACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGa                                                                                                                                 >  2:354527/1‑150 (MQ=255)
                     gATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGGAAATTCTAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGc                                                                                                                              >  2:767693/1‑150 (MQ=255)
                        aaaCAGCCCTGCTTTATATGATTTAGGGGCATGATTATCCGTTTCTGCAGGTTTTTTAATCGAACGTTTGGTGGGTAAAAGGGGGGTTTATTTGGGGT‑TTAGAGATAAATTCGAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCAc                                                                                                                            <  1:635156/149‑1 (MQ=255)
                        aaaCACCCCTGTTTTATATGAGTTATGGGAATGTTTATCCGTTTCGGCGGGGTTTTAAAGCGGAGGATTAGGGGGTAAAATGAGGGTTTTTTTGGGGG‑TTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGt                                                                                                                          <  1:48854/151‑1 (MQ=255)
                                  gCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTg                                                                                                                     >  1:1265615/1‑146 (MQ=255)
                                     ttATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGTTTTTTAATCGGACGATTGGTGGGTAAAATGGGGGTTTTTTTGGGGG‑TTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaaa                                                                                                                <  1:590506/148‑1 (MQ=255)
                                         atGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCt                                                                                                         >  1:127619/1‑151 (MQ=255)
                                          tGAGTTATCGGCATGCTTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGGtt                                                                                                        >  1:407465/1‑151 (MQ=255)
                                           gAGTTAGGGGCTTGATTATCTTTTTCTGCGGGGTTTTTAAGCGGACGATTAGGGGGGAAAATGGGGGTTTATTTGGGGT‑TTACAGGAAATTCCCAGGCAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCtt                                                                                                        <  2:347796/150‑1 (MQ=255)
                                              ttACCGGCATGATTATCCTTTTCTGCGGGGTTTTTAAGCGGATGATTGGTGGGGGAAAGGGGGGTTTTTTTGGGGG‑TAACAGGTAAATCCCAGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTa                                                                                                      <  1:1332957/149‑1 (MQ=255)
                                               tATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt                                                                                                   >  1:235003/1‑151 (MQ=255)
                                                aTCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGtt                                                                                                  >  2:183064/1‑151 (MQ=255)
                                                   ggCATGATTATCTTTTTCTGAGGGGTTTTTAATGGGACATTTAGTGGGTGAAAGGGGGTGTTTTTTGGGTG‑TTAGAGGAAAATTCCAGAAAAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTg                                                                                                 <  1:1006637/149‑1 (MQ=255)
                                                    gCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg                                                                                              >  1:975160/1‑151 (MQ=255)
                                                        gATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                          >  2:935530/1‑151 (MQ=255)
                                                        gATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                          >  2:842894/1‑151 (MQ=255)
                                                         aTTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgc                                                                                         >  2:197670/1‑151 (MQ=255)
                                                          ttctCCGTTTCGGCAGGTTTTTTAAGCGGACGTTTGGTGGGGAAAATGGGGGTTTTTTTGGGGG‑TTACAGGTAATTCCCAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt                                                                                           <  1:525877/145‑1 (MQ=255)
                                                          ttATCCTTTTCTGGAGGTTTTTTAAGCGGACGATTAGTGGGTAAAAGGGGGGTTTTTTTGGGGG‑TTACAGTAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                          <  2:1132986/149‑1 (MQ=255)
                                                                 tttgtggggggTTTTTACTGGGGGGAGTGGGGGGTGAAATGAGGGTTTTTTGGGGGG‑TAGCATTCAAGTCCCAGGAAGAAAAAACCACCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTa                                                                                   <  1:900346/141‑1 (MQ=255)
                                                                  ttCTGCAGGTTTTTAAATCGGAAGATTAGGGGGTAAAATGAGGGTTTTTTGGGGGG‑TTACAGTTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaaa                                                                                 <  2:1245691/150‑1 (MQ=255)
                                                                   tcgggaGGTTTTTTAATCGGACGTTGGGGGGGTAAAAGGGGGGTTTTTTTGGGGT‑TTACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                                <  2:846276/145‑1 (MQ=255)
                                                                    cTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGc                                                                              >  1:547063/1‑151 (MQ=255)
                                                                         ggttttttAATCGGACGTTTAGTGGGTAAAAGGGGGGTTTATTTGGGGG‑TTACAGGTAAATCCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                           <  2:1144589/146‑1 (MQ=255)
                                                                         gggTTTTTAATCGGACGATTGGGGGGTGAAAGGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTa                                                                          <  2:1092406/150‑1 (MQ=255)
                                                                         gggTTTTAAATCGAACGATTAGGGGGTAAAATGGGGGGTTATTGGGGGT‑TTACAGGTAATTCCCGGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                           <  2:215223/149‑1 (MQ=255)
                                                                            tttttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                      >  2:530454/1‑151 (MQ=255)
                                                                                taaggggACGATTAGGGGGTAAAATGAGGGTTTATTGGGGGG‑TTACAGGTAAATCCGAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                    <  2:1081709/144‑1 (MQ=38)
                                                                                 aaTCGGAGGATGGGTGGGGAAAATGAGGGTTTTTTTGGGGG‑TTAAAGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                 <  1:801546/151‑1 (MQ=37)
                                                                                 aaTCGAACGATTGGTGGGTGAAA‑GGGGGGTTTTTTGGGGGTTTAGAGGAAAATCCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                   <  2:1132829/149‑1 (MQ=37)
                                                                                   acGGACGGTTGGGGGGGAAAATGGGGGTTTTTTTGGGGG‑TTAAAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                               <  2:1130095/150‑1 (MQ=37)
                                                                                    cGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                               >  2:1038922/1‑150 (MQ=255)
                                                                                     ggACGTTTAGTGGGTAAAATGGGGGTTTTTTTGGGGT‑TTACAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                               <  2:1018579/149‑1 (MQ=37)
                                                                                     ggACGATTTGTGGGTGAAATGGGGGTTTTTTTGGGGG‑TTAAAGGTAACTGCCAGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaa                                                             <  1:1146491/151‑1 (MQ=37)
                                                                                      gACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaa                                                             >  2:737456/1‑150 (MQ=255)
                                                                                       acaataagtggtaaaaaTGAGTGTTTTTTTGGGGG‑TTACAGGAAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaaa                                                           <  1:1238177/137‑1 (MQ=37)
                                                                                       acaaTTAGTGGGTAAAAGGGGGGGTTTTTTGGGGG‑TTAGAGATAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                            <  1:675938/147‑1 (MQ=37)
                                                                                         gATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                         >  2:1200317/1‑151 (MQ=255)
                                                                                              gggggTGAAATGAGGGGTTATTGGGGGT‑TTACAGGTAATTTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa                                                    <  2:1181226/149‑1 (MQ=38)
                                                                                                gggTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaaa                                                   >  1:90494/1‑150 (MQ=255)
                                                                                                 ggTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAtt                                                 >  2:1242687/1‑151 (MQ=255)
                                                                                                     aaagggggggTTTTTTGGGGG‑TAAAAGTAAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAACTCAAt                                             <  1:185274/145‑1 (MQ=25)
                                                                                                     aaaTGGGGGTTTTTTTGGGGT‑TTACAGGTAAATCCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAtt                                                 <  2:1178969/147‑1 (MQ=37)
                                                                                                     aaaTGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                             >  1:201812/1‑151 (MQ=255)
                                                                                                           gggtTTTTTTGGGGT‑TAACAGAAAAATTCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                       <  1:511512/151‑1 (MQ=25)
                                                                                                           gggtTTTTTTGGGGG‑TTAAAGTTAAATTCCAGGCAAAAAAACACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                           <  2:155155/147‑1 (MQ=37)
                                                                                                               ttttttGGGGG‑TTAGAGGTAAATTCAAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCaa                                    <  1:255933/147‑1 (MQ=37)
                                                                                                                tATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                  >  2:1016302/1‑151 (MQ=255)
                                                                                                                tATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                  >  2:73210/1‑151 (MQ=255)
                                                                                                                  tttGGGGG‑TTACAGTTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                                 <  1:1302823/150‑1 (MQ=255)
                                                                                                                  tttGGGGG‑TTAAAGTTAAATTCGAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                  <  1:167318/149‑1 (MQ=255)
                                                                                                                     ggggg‑TTACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCaa                               <  2:1176866/149‑1 (MQ=255)
                                                                                                                      gggt‑TTACAGTAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt                            <  2:201837/151‑1 (MQ=37)
                                                                                                                       ggggtaccagaaaaaTTCCGGCAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt                            <  1:62394/139‑1 (MQ=25)
                                                                                                                         g‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                          >  2:825492/1‑150 (MQ=255)
                                                                                                                           ttACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                          >  1:1227887/1‑149 (MQ=255)
                                                                                                                             aCAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGc                      >  2:57627/1‑151 (MQ=255)
                                                                                                                              cAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGcc                     >  2:39116/1‑151 (MQ=255)
                                                                                                                                  tAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt                 >  2:1297758/1‑151 (MQ=255)
                                                                                                                                    actcccgggCAGAAAAATCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCAGATGCAGCATAAACGCCTTATCCGGGCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGTCCTGAt                 <  1:310594/142‑1 (MQ=37)
                                                                                                                                     aTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa                >  1:882689/1‑149 (MQ=255)
                                                                                                                                        ccAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc           >  2:36862/1‑151 (MQ=255)
                                                                                                                                           ggCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagc        >  2:889273/1‑151 (MQ=255)
                                                                                                                                            gCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcg       >  1:778979/1‑151 (MQ=255)
                                                                                                                                             cAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc      >  2:757290/1‑151 (MQ=255)
                                                                                                                                              aGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc      >  2:863159/1‑150 (MQ=255)
                                                                                                                                                 aaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCa  <  2:647283/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCA  >  NC_000913/2001769‑2002063

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