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

AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCA  >  NC_000913/2001769‑2002063
                                                                                                                                                      |                                                                                                                                                 
aGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGTTTATCCGTTTCTGCAGGGGTTTTAATCGGACCATTAGGGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAc                                                                                                                                                  >  1:1049547/1‑151 (MQ=255)
   ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAcccc                                                                                                                                               <  2:926588/151‑1 (MQ=255)
   ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAcc                                                                                                                                                 <  2:577568/149‑1 (MQ=255)
          agTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTg                                                                                                                                        >  2:50529/1‑151 (MQ=255)
                ttCACGAAAAAAAGCCCTGCGTTATATGATTTATCGGCATGATTATCCGTTTCTGCGGGGTTTTAAACCGGACGATTAGGGGGTAAAAGGGGGGGTTTTTTGGGGGTTACAGGAAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCggg                                                                                                                                   <  1:1145055/150‑1 (MQ=255)
                  cACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGa                                                                                                                                 >  2:354527/1‑150 (MQ=255)
                     gATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGGAAATTCTAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGc                                                                                                                              >  2:767693/1‑150 (MQ=255)
                        aaaCAGCCCTGCTTTATATGATTTAGGGGCATGATTATCCGTTTCTGCAGGTTTTTTAATCGAACGTTTGGTGGGTAAAAGGGGGGTTTATTTGGGGTTTAGAGATAAATTCGAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCAc                                                                                                                            <  1:635156/149‑1 (MQ=255)
                        aaaCACCCCTGTTTTATATGAGTTATGGGAATGTTTATCCGTTTCGGCGGGGTTTTAAAGCGGAGGATTAGGGGGTAAAATGAGGGTTTTTTTGGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGt                                                                                                                          <  1:48854/151‑1 (MQ=255)
                                  gCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTg                                                                                                                     >  1:1265615/1‑146 (MQ=255)
                                     ttATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGTTTTTTAATCGGACGATTGGTGGGTAAAATGGGGGTTTTTTTGGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaaa                                                                                                                <  1:590506/148‑1 (MQ=255)
                                         atGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCt                                                                                                         >  1:127619/1‑151 (MQ=255)
                                          tGAGTTATCGGCATGCTTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGGtt                                                                                                        >  1:407465/1‑151 (MQ=255)
                                           gAGTTAGGGGCTTGATTATCTTTTTCTGCGGGGTTTTTAAGCGGACGATTAGGGGGGAAAATGGGGGTTTATTTGGGGTTTACAGGAAATTCCCAGGCAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCtt                                                                                                        <  2:347796/150‑1 (MQ=255)
                                              ttACCGGCATGATTATCCTTTTCTGCGGGGTTTTTAAGCGGATGATTGGTGGGGGAAAGGGGGGTTTTTTTGGGGGTAACAGGTAAATCCCAGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTa                                                                                                      <  1:1332957/149‑1 (MQ=255)
                                               tATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt                                                                                                   >  1:235003/1‑151 (MQ=255)
                                                aTCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGtt                                                                                                  >  2:183064/1‑151 (MQ=255)
                                                   ggCATGATTATCTTTTTCTGAGGGGTTTTTAATGGGACATTTAGTGGGTGAAAGGGGGTGTTTTTTGGGTGTTAGAGGAAAATTCCAGAAAAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTg                                                                                                 <  1:1006637/149‑1 (MQ=255)
                                                    gCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg                                                                                              >  1:975160/1‑151 (MQ=255)
                                                        gATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                          >  2:935530/1‑151 (MQ=255)
                                                        gATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                          >  2:842894/1‑151 (MQ=255)
                                                         aTTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgc                                                                                         >  2:197670/1‑151 (MQ=255)
                                                          ttctCCGTTTCGGCAGGTTTTTTAAGCGGACGTTTGGTGGGGAAAATGGGGGTTTTTTTGGGGGTTACAGGTAATTCCCAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt                                                                                           <  1:525877/145‑1 (MQ=255)
                                                          ttATCCTTTTCTGGAGGTTTTTTAAGCGGACGATTAGTGGGTAAAAGGGGGGTTTTTTTGGGGGTTACAGTAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                          <  2:1132986/149‑1 (MQ=255)
                                                                 tttgtggggggTTTTTACTGGGGGGAGTGGGGGGTGAAATGAGGGTTTTTTGGGGGGTAGCATTCAAGTCCCAGGAAGAAAAAACCACCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTa                                                                                   <  1:900346/141‑1 (MQ=255)
                                                                  ttCTGCAGGTTTTTAAATCGGAAGATTAGGGGGTAAAATGAGGGTTTTTTGGGGGGTTACAGTTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaaa                                                                                 <  2:1245691/150‑1 (MQ=255)
                                                                   tcgggaGGTTTTTTAATCGGACGTTGGGGGGGTAAAAGGGGGGTTTTTTTGGGGTTTACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                                <  2:846276/145‑1 (MQ=255)
                                                                    cTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGc                                                                              >  1:547063/1‑151 (MQ=255)
                                                                         ggttttttAATCGGACGTTTAGTGGGTAAAAGGGGGGTTTATTTGGGGGTTACAGGTAAATCCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                           <  2:1144589/146‑1 (MQ=255)
                                                                         gggTTTTTAATCGGACGATTGGGGGGTGAAAGGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTa                                                                          <  2:1092406/150‑1 (MQ=255)
                                                                         gggTTTTAAATCGAACGATTAGGGGGTAAAATGGGGGGTTATTGGGGGTTTACAGGTAATTCCCGGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                           <  2:215223/149‑1 (MQ=255)
                                                                            tttttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                      >  2:530454/1‑151 (MQ=255)
                                                                                taaggggACGATTAGGGGGTAAAATGAGGGTTTATTGGGGGGTTACAGGTAAATCCGAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                    <  2:1081709/144‑1 (MQ=38)
                                                                                 aaTCGGAGGATGGGTGGGGAAAATGAGGGTTTTTTTGGGGGTTAAAGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                 <  1:801546/151‑1 (MQ=37)
                                                                                 aaTCGAACGATTGGTGGGTGAAAGGGGGGTTTTTTGGGGGTTTAGAGGAAAATCCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                   <  2:1132829/149‑1 (MQ=37)
                                                                                   acGGACGGTTGGGGGGGAAAATGGGGGTTTTTTTGGGGGTTAAAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                               <  2:1130095/150‑1 (MQ=37)
                                                                                    cGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                               >  2:1038922/1‑150 (MQ=255)
                                                                                     ggACGTTTAGTGGGTAAAATGGGGGTTTTTTTGGGGTTTACAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                               <  2:1018579/149‑1 (MQ=37)
                                                                                     ggACGATTTGTGGGTGAAATGGGGGTTTTTTTGGGGGTTAAAGGTAACTGCCAGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaa                                                             <  1:1146491/151‑1 (MQ=37)
                                                                                      gACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaa                                                             >  2:737456/1‑150 (MQ=255)
                                                                                       acaataagtggtaaaaaTGAGTGTTTTTTTGGGGGTTACAGGAAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaaa                                                           <  1:1238177/137‑1 (MQ=37)
                                                                                       acaaTTAGTGGGTAAAAGGGGGGGTTTTTTGGGGGTTAGAGATAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                            <  1:675938/147‑1 (MQ=37)
                                                                                         gATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                         >  2:1200317/1‑151 (MQ=255)
                                                                                              gggggTGAAATGAGGGGTTATTGGGGGTTTACAGGTAATTTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa                                                    <  2:1181226/149‑1 (MQ=38)
                                                                                                gggTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaaa                                                   >  1:90494/1‑150 (MQ=255)
                                                                                                 ggTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAtt                                                 >  2:1242687/1‑151 (MQ=255)
                                                                                                     aaagggggggTTTTTTGGGGGTAAAAGTAAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAACTCAAt                                             <  1:185274/145‑1 (MQ=25)
                                                                                                     aaaTGGGGGTTTTTTTGGGGTTTACAGGTAAATCCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAtt                                                 <  2:1178969/147‑1 (MQ=37)
                                                                                                     aaaTGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                             >  1:201812/1‑151 (MQ=255)
                                                                                                           gggtTTTTTTGGGGTTAACAGAAAAATTCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                       <  1:511512/151‑1 (MQ=25)
                                                                                                           gggtTTTTTTGGGGGTTAAAGTTAAATTCCAGGCAAAAAAACACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                           <  2:155155/147‑1 (MQ=37)
                                                                                                               ttttttGGGGGTTAGAGGTAAATTCAAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCaa                                    <  1:255933/147‑1 (MQ=37)
                                                                                                                tATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                  >  2:1016302/1‑151 (MQ=255)
                                                                                                                tATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                  >  2:73210/1‑151 (MQ=255)
                                                                                                                  tttGGGGGTTACAGTTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                                 <  1:1302823/150‑1 (MQ=255)
                                                                                                                  tttGGGGGTTAAAGTTAAATTCGAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                  <  1:167318/149‑1 (MQ=255)
                                                                                                                     gggggTTACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCaa                               <  2:1176866/149‑1 (MQ=255)
                                                                                                                      gggtTTACAGTAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt                            <  2:201837/151‑1 (MQ=37)
                                                                                                                       ggggtaccagaaaaaTTCCGGCAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt                            <  1:62394/139‑1 (MQ=25)
                                                                                                                         gTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                          >  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)
                                                                                                                                                      |                                                                                                                                                 
AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCA  >  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: