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

AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGG  >  NC_000913/2001769‑2002065
                                                                                                                                                        |                                                                                                                                                   
aGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAc                                                                                                                                                    <  1:748044/151‑1 (MQ=255)
   ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAcccc                                                                                                                                                 >  2:728256/1‑151 (MQ=255)
    gTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccg                                                                                                                                                >  2:805973/1‑151 (MQ=255)
    gTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccg                                                                                                                                                >  1:325739/1‑151 (MQ=255)
      ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgc                                                                                                                                               <  2:1210950/150‑1 (MQ=255)
      ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGGTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgcc                                                                                                                                              >  1:1086145/1‑151 (MQ=255)
            tGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGc                                                                                                                                        >  1:242187/1‑151 (MQ=255)
             gAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCg                                                                                                                                       >  1:1039251/1‑151 (MQ=255)
                 tCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCgggg                                                                                                                                    >  1:133817/1‑150 (MQ=255)
                     gATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCGGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGc                                                                                                                                >  2:110905/1‑150 (MQ=255)
                       tAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACg                                                                                                                             >  2:176642/1‑151 (MQ=255)
                        aaaCAGCCCTGCGTTATATGATTTATCGGCATTATTATCCTTTTCTGCAGGTTTTTAAATCGGAGGATTAGGGGGTAAAAGGGGGGTTTATTGGGGGTTTACAGTTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGc                                                                                                                                <  2:1462874/147‑1 (MQ=255)
                        aaaCAGCCCTGCGTTATATGAGTGATCGGTATGATTATCTGTTTCTGGGGGTTTTTTAATGGGACGATTAGTGGGTGAAAGGAGGGTTTTTTTGGGGGTTAAAGAAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCa                                                                                                                               <  1:1328913/148‑1 (MQ=255)
                        aaaCAGCCCGGCTTTATATGAGTTATGGGCATGATTATCCTTTTCGGCGGGTTTTTTAATCGGACGATTGGTGGGTGAAATGGGGGTTTTTTTGGGGGTTACAGTTAAATCCAAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCAc                                                                                                                              <  1:452601/149‑1 (MQ=255)
                         aaCGGCCCTGCGTTATAGGATTTATCGGAATGATTATCCGTTTCTGCGGGTTTTTAAATCGGACGATTGGGGGGTAAAAGGGGGGTTTATTTGGGGGTTACAGGTAAATCCCGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGtt                                                                                                                           <  2:556864/151‑1 (MQ=255)
                         aaCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGtt                                                                                                                           >  2:1124697/1‑151 (MQ=255)
                               ccTGCTTTATATGAGTTATCGGCATGATTATCCGTTTCTGCGGGGTTTTTAATCGGACGATTAGGGGGTAAAATGAGGGTTTTTTGGGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGAc                                                                                                                     <  2:249641/151‑1 (MQ=255)
                                cTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTCGTGGGTGAACTGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAACACGTTGCTGAc                                                                                                                     >  1:898745/1‑150 (MQ=255)
                                  gCGTTATATGATTTATCGGCATGATTATCCGTTTCGGCGGGTTTTTAAATCGGACGATTGGTGGGTGAAATGGGGGTTTATTTGGGGGTTACAGGTAATTCCCGGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                                    <  1:985806/149‑1 (MQ=255)
                                  gCGTTATATGAGTTATGGGCATGATTATCCTTTTCTGCGGGTTTTTTAATCGGAGGATTGGGGGGTGAAATGAGGGTTTTTTTGGGGGTAAAAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGAc                                                                                                                     <  1:1247490/148‑1 (MQ=255)
                                   cGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAt                                                                                                                 >  1:552517/1‑151 (MQ=255)
                                     ttATATGAGTTATCGGCATGATTTTCGTTTTCTGCGGGTTTTTTAACGGGACGATGAGGGGGTGAAAGGAGGGTTTTTTTGGGGGTAAAAGTAAATTCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                                <  1:1278684/150‑1 (MQ=255)
                                     ttAAATGAGTTATCGGCATGATTGCCGTTTTCGGCAGGTTTTTTAACGGAACGATTGGTGGGGGAAAGGGGGGTTTATTTGGGGTTTGAGGATAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGc                                                                                                                         <  2:1481039/141‑1 (MQ=255)
                                         atGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCt                                                                                                           >  1:649192/1‑151 (MQ=255)
                                         atGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCt                                                                                                           >  2:970251/1‑151 (MQ=255)
                                           gatttATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAAGGGGGGTTTTTTTGGGGGTTACAGTAAAATCCGAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcg                                                                                                             <  1:1289552/144‑1 (MQ=255)
                                           gAGTTATCGGCATGATTATCCGTTTCTGCAGGTTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGTTTTTTTGGGGGTTACAGGTAAATTCCAGGAAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt                                                                                                         <  2:1208095/151‑1 (MQ=255)
                                           caGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt                                                                                                         >  2:477430/2‑151 (MQ=255)
                                              ttttgggaaTGATTCCCCGTTTCGCAAGGTTTTTAATTCGGAAGTTTAGGGGGTAAAAGGGGGGTTTTTTTGGGGGTTAAAGTAAATTCCCAGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTAt                                                                                                       <  1:550839/142‑1 (MQ=255)
                                               tATCGGCATGATTTTCCTTTTCTGCGGGTTTTTTAATGGGACGTTTAGTGGGTAAAAGGGGGGTTTTTTTGGGGTTTACAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATg                                                                                                      <  2:1026550/150‑1 (MQ=255)
                                                  cGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGc                                                                                                  >  1:973946/1‑151 (MQ=255)
                                                  cGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGc                                                                                                  >  2:1455531/1‑151 (MQ=255)
                                                  cGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGc                                                                                                  >  2:875477/1‑151 (MQ=255)
                                                   ggaaTGATTATCCTTTTCTGCGGGTTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGTTTATTGGGGGTTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt                                                                                                     <  2:801893/144‑1 (MQ=255)
                                                   ggaaTGATTATCCGTTTCTGCAGGTTTTTTAATCGGACGATTAGTGGGTGAAATGGGGGTTTTTTGGGGGTTTACAGGAAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGc                                                                                                  <  2:535321/147‑1 (MQ=255)
                                                   ggCATGATTATCCGTTTCTGGAGGGTTTTTAATGGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGtt                                                                                                    <  1:579605/148‑1 (MQ=255)
                                                   ggCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGcc                                                                                                 >  1:758062/1‑151 (MQ=255)
                                                    gCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg                                                                                                >  2:1209918/1‑151 (MQ=255)
                                                      aTGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGCGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGAGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                               >  2:381721/1‑150 (MQ=255)
                                                          ttATCCTTTTCTGCGGGTTTTTTAATCGGACGATTAGGGGGTAAAATGGGGGTTTTTTTGGGGTTTACAGGAAATTTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgc                                                                                           <  2:395586/150‑1 (MQ=255)
                                                          ttATCCTTTTCTGCAGGTTTTTTAACCGGACTATTAGTGGGTGAAATGAGGGTTTTTTGGGGGGTTACAGGAAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg                                                                                                <  2:1500376/145‑1 (MQ=255)
                                                                gTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGATATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAACACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggcg                                                                                       >  2:156588/1‑148 (MQ=255)
                                                                 tttCTGCAGGTTTTTTAATGGGACGTTTAGGGGGTGAAATGAGGGTTTATTGGGGGTTTACAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTa                                                                                     <  2:1190545/149‑1 (MQ=255)
                                                                 tttCTGCAGGTTTTTAAATCGGACGATTAGTGGGTAAAATGGGGGTTTTTTGGGGGTTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaa                                                                                    <  2:1020075/150‑1 (MQ=255)
                                                                 tttCTGCAGGGTTTTTAAGCGGACGTTTGGGGGGTAAAAGGGGGGTTTTTTTGGGGGTTACAGTTAAATTCCAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGt                                                                                      <  1:381385/148‑1 (MQ=255)
                                                                  ttttGCAGGTTTTTTAAAGGAAGAATTGGGGGGAGAAAGGGGGGTTTGTTGGGGGGTAAAAGTAAAATCCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTa                                                                                     <  2:1252700/145‑1 (MQ=255)
                                                                  ttCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                                  >  1:920840/1‑151 (MQ=255)
                                                                   ttcggaGGGTTTTTAATTCGGACTATTGGGGGGTAAAAGGAGGGTTTATTTGGGGTTAAGAGGTAATTTCCAGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTa                                                                                     <  1:220945/143‑1 (MQ=255)
                                                                   tCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                                  >  2:358255/1‑150 (MQ=255)
                                                                     tGCGGGTTTTTTAATCGGACGATTAGGGGGTAAAATGGGGGTTTTTTTGGGGTTAAAAGGTAATTCCAAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGc                                                                                <  1:410538/150‑1 (MQ=255)
                                                                      gcggggTTTTTAATCGGACGTTTAGGGGGTAAAATGGGGGGTTTTTTGGGGGTAACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACg                                                                                 <  2:1209930/145‑1 (MQ=255)
                                                                             ttttAATCGGACGATTAGTGGGTAAAAGGAGGGGTTTTTGGGGGGTTACAGGAAAATCCCGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                        <  2:1157974/150‑1 (MQ=38)
                                                                                  aTCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                   >  1:1490661/1‑150 (MQ=255)
                                                                                        aCGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                              >  1:1431061/1‑150 (MQ=255)
                                                                                        aCGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                              >  2:964014/1‑150 (MQ=255)
                                                                                          gATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                           >  2:1148458/1‑151 (MQ=255)
                                                                                               gTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCa                                                       >  2:278305/1‑150 (MQ=255)
                                                                                                tGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa                                                      >  2:1032570/1‑150 (MQ=255)
                                                                                                 gggTGAAATGAGGGTTTTTTTGGGGGTTACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAt                                                    <  1:559679/151‑1 (MQ=255)
                                                                                                 gggTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAt                                                    >  2:289736/1‑151 (MQ=255)
                                                                                                    tGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTc                                                  >  1:859637/1‑150 (MQ=255)
                                                                                                     gAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa                                                 >  2:1431539/1‑150 (MQ=255)
                                                                                                      aaaTGGGGGTTTTTTTGGGGGTTACAGAAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAtt                                                   <  1:818534/147‑1 (MQ=37)
                                                                                                      aaaTGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                               >  2:1121051/1‑151 (MQ=255)
                                                                                                      aaaGGAGGGTTTTTTTGGGGGTTAAAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCaa                                                <  1:291566/150‑1 (MQ=37)
                                                                                                       aaTGGGGGTTTTTTTGGGGGTTAAAGATAATTTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATa                                              <  1:385087/151‑1 (MQ=25)
                                                                                                       aaTGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                               >  1:176968/1‑150 (MQ=255)
                                                                                                        aTGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATa                                              >  2:700225/1‑150 (MQ=255)
                                                                                                        aTGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATa                                              >  1:163057/1‑150 (MQ=255)
                                                                                                         tGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                             >  1:1400982/1‑150 (MQ=255)
                                                                                                          gAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt                                           >  2:1110505/1‑151 (MQ=255)
                                                                                                            gggtTTATTGGGGGTTTACAGTTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                          <  2:168194/150‑1 (MQ=37)
                                                                                                             ggtttttttGGGGGTTACAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                                        <  1:408792/145‑1 (MQ=255)
                                                                                                                ttttttGGGGGTTACAGGTAAATTCCAGGCAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCa                                       <  2:558730/146‑1 (MQ=255)
                                                                                                                   tttGGGGGTTAAAGAAAATTCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                    <  1:960409/149‑1 (MQ=25)
                                                                                                                     ttggggttaaaagaaaaaTTCAAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                                  <  1:928032/135‑1 (MQ=25)
                                                                                                                      gggggTAACAGGAAAATTCAGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATTCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAAc                                <  2:1277386/150‑1 (MQ=37)
                                                                                                                       ggggTTACAGTTAAATTCAAGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt                              <  1:43375/151‑1 (MQ=37)
                                                                                                                         ggTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                            >  2:268490/1‑151 (MQ=255)
                                                                                                                         ggTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                            <  1:1515540/151‑1 (MQ=255)
                                                                                                                          gtttACGGTTAAATCCCAGACAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                            <  2:1552684/149‑1 (MQ=37)
                                                                                                                          gTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                            >  2:355667/1‑150 (MQ=255)
                                                                                                                                 ggaaatttCCAGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg                     <  1:805927/145‑1 (MQ=255)
                                                                                                                                    aaaTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt                   >  1:1174567/1‑150 (MQ=255)
                                                                                                                                      aTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa                  >  1:923646/1‑149 (MQ=255)
                                                                                                                                        tcccAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAag                <  1:277083/148‑1 (MQ=255)
                                                                                                                                          ccAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcg              <  2:1183284/150‑1 (MQ=255)
                                                                                                                                          caagaaaGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAag                <  1:87009/142‑1 (MQ=255)
                                                                                                                                           cAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc             >  1:479379/1‑150 (MQ=255)
                                                                                                                                           cAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcg              <  2:339917/149‑1 (MQ=255)
                                                                                                                                            aGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc             >  2:787362/1‑149 (MQ=255)
                                                                                                                                              gCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagc          >  1:1207959/1‑150 (MQ=255)
                                                                                                                                                aGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc        <  2:297699/150‑1 (MQ=255)
                                                                                                                                                aGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc        >  1:1160284/1‑150 (MQ=255)
                                                                                                                                                 gAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAt      >  1:482205/1‑151 (MQ=255)
                                                                                                                                                   aaaaaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATc     >  1:1166453/1‑150 (MQ=255)
                                                                                                                                                     aaaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCTGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAgg  >  2:944646/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                   
AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGG  >  NC_000913/2001769‑2002065

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: