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

GTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGC  >  NC_000913/2001773‑2002066
                                                                                                                                                 |                                                                                                                                                     
gTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCccg                                                                                                                                                 <  2:896791/151‑1 (MQ=255)
       gTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGAAGAAAAAAC‑CCCCGCCGGTg                                                                                                                                           <  1:443391/150‑1 (MQ=255)
       gTGAATTCACGAAAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGTTTTTTAATCGGACGATTAGGGGGTAAAAGGGGGGGTTTTTTGGGGGTTACAGGTAAATTCCAGGCAAAAAAAAC‑CCCCGCCGGTg                                                                                                                                           <  1:571600/150‑1 (MQ=255)
           aTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCggg                                                                                                                                      >  2:670387/1‑151 (MQ=255)
              cACGATAAACAGCCCTGCTTTATATGAGTTATGGGCATGATTATCCGTTTCTGCGGGGTTTTTAATCGGACGATTAGTGGGTAAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCgggg                                                                                                                                     <  1:521892/149‑1 (MQ=255)
              cACGATAAACAGCCCTGCGTTATATGAGTTATGGGCATGATTATCCGTTTCTGCGGGTTTTTTAATCGGACGATTGGGGGGTGAAAGGAGGGGTTATTTGGGGTTTACAGGTAAATTCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGaa                                                                                                                                   <  2:746976/151‑1 (MQ=255)
              cACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGa                                                                                                                                    >  2:1015595/1‑150 (MQ=255)
              aaCGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGa                                                                                                                                    >  1:338276/2‑150 (MQ=255)
                  aTAAACAGCCCTGCGTTATATGATTTATCGGCATGATTATCCTTTTCTGCGGGTTTTTTAACCGGACGATTAGGGGGTGAAATGGGGGTTTATTGGGGGTTTACAGGTAATTCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGc                                                                                                                                 <  2:464925/149‑1 (MQ=255)
                   tAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACg                                                                                                                              >  1:836900/1‑151 (MQ=255)
                    aaccAGCCCTGGGTTATATGAGTTATCGGCATGATTTCCCGTTTCTGCGGGTTTTTTAATCGGACGATTAGTGGGTGAAAGGGGGGTTTTTTTGGGGGTTACAGTTAAATTCCGGGCAGAAAAAAC‑CCCCGCCGGTGGCgggg                                                                                                                                     <  1:47788/140‑1 (MQ=255)
                    aaaCAGCCCTGCTTTATATGATTTTTCGGCATGATTATCCGTTTCTGCAGGTTTTTTAATCGGACGATTAGGGGGTAAAAGGAGGGGTTTTTGGGGGGTTACAGTTAAATTCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCAc                                                                                                                               <  2:821759/149‑1 (MQ=255)
                        aGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCt                                                                                                                         >  2:657181/1‑151 (MQ=255)
                         gCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTg                                                                                                                        >  1:433759/1‑151 (MQ=255)
                               ctttttatGATTTACCGGCATAATTCTCCGTTTCGGCGGGGTTTTTAATGGGAGGATTAGGGGGTGAAAGGGGGGGTTTTTTGGGGGTTAAAGGTAAATTCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaa                                                                                                                    <  2:555953/144‑1 (MQ=255)
                               cGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAt                                                                                                                  >  2:195416/1‑151 (MQ=255)
                                  tataTGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgc                                                                                                               >  2:440263/1‑151 (MQ=255)
                                     atGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCt                                                                                                            >  1:513952/1‑151 (MQ=255)
                                           tATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt                                                                                                      >  2:522698/1‑151 (MQ=255)
                                              cGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGc                                                                                                   >  1:244110/1‑151 (MQ=255)
                                               ggCATGATTATCCGTTTCTGCGGGGTTTTTAATCGGACGATTGGGGGGTAAAATGGGGGTTTATTGGGGGGTTACAGGTAAATTCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTg                                                                                                    <  2:1024146/149‑1 (MQ=255)
                                                   tttttttCCTTTTCTGCAGGGTTTTTAATCGTACTATTTTTGTGTAAAATGATTGTTTATTTGGGTTTTCCAGTTAAATTCCACCCAGAAAAAAC‑CCCCGCCGGTGGCGGTGAAGCACGTTGCTGACAAATTGCGCTTTATGTTTCCTGa                                                                                               <  1:668626/144‑1 (MQ=255)
                                                   tGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt                                                                                              >  1:260148/1‑151 (MQ=255)
                                                    gATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                             >  2:699704/1‑151 (MQ=255)
                                                      ttATCCGTTTCTGCGGGTTTTTTAATGGGACGATTAGTGGGTAAAATGAGGGTTTATTTGGGGTTTACGGGTAAATTCAAGGCAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt                                                                                              <  2:115982/148‑1 (MQ=255)
                                                          ccttttCTGCGGGGTTTTTAATCGGACGATTAGGGGGTAAAATGAGGGTTTATTGGGGGGTTACAGGTAAATTCCGGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggc                                                                                         <  2:814386/146‑1 (MQ=255)
                                                                 tGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGcc                                                                                >  1:697423/1‑151 (MQ=255)
                                                                 ggggggttttttAGCCGAACGATTGGGGGGTAAAAGGGGGGTTTATGGGGGGTTAAAAGATAAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGcc                                                                                <  2:853125/144‑1 (MQ=255)
                                                                   cAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                              >  2:475965/1‑151 (MQ=255)
                                                                     ggttttttAATCGGCCGTTTAGGGGGTAAAAGGAGGGTTTTTTGGGGGTTTAAGGGAAAATTCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGc                                                                                 <  2:957049/143‑1 (MQ=255)
                                                                              accGGACGATTAGGGGGTAAAATGAGGGTTTTTTTGGGGGTTACAGGTAAATCCCAGGCAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                   <  1:522602/149‑1 (MQ=38)
                                                                                   aCGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                               >  1:697719/1‑150 (MQ=255)
                                                                                    cGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAt                                                             >  1:244927/1‑151 (MQ=255)
                                                                                     gtttAGGGGGTAAAATGGGGGGTTATTGGGGGGTTACGGGTAAATTCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAt                                                             <  2:106551/148‑1 (MQ=38)
                                                                                          gTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCa                                                        >  2:683269/1‑150 (MQ=255)
                                                                                            gggTGAAATGGGGGGTTGTTTGGGGGTTACAGGTAAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaaa                                                      <  1:1028783/150‑1 (MQ=255)
                                                                                            gggTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAt                                                     >  2:150891/1‑151 (MQ=255)
                                                                                            gggTAAAATGAGGGTTTATTTGGGGGTTACAGTTAAATTCCAGGAAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                         <  1:191185/147‑1 (MQ=255)
                                                                                               tGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTc                                                   >  1:241319/1‑150 (MQ=255)
                                                                                                gAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGGCGGGGGCGGGGGAAGACGGTGCCGGCCAATTGCCGTTTATGTTGCCGGGTGCGGCGTAAACCCCTTATCCGGCCCACAAAAATGGGCCAAtt                                                    >  2:634899/1‑148 (MQ=255)
                                                                                                  aaTGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                                >  2:306163/1‑150 (MQ=255)
                                                                                                     gAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt                                            >  1:795579/1‑151 (MQ=255)
                                                                                                       ggggTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                          >  2:784428/1‑151 (MQ=255)
                                                                                                          tttATTTGGGGTTTAAAGTTAAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCa                                        <  1:699586/149‑1 (MQ=37)
                                                                                                           ttATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAt                                      >  2:733505/1‑151 (MQ=255)
                                                                                                            tATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                     >  1:746080/1‑151 (MQ=255)
                                                                                                              tttGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                                    <  2:387511/150‑1 (MQ=255)
                                                                                                              tttGGGGGTTACAGGTAAATCCCGGGCAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAt                                      <  1:311491/149‑1 (MQ=255)
                                                                                                              ttgggggtTTACAGGAAAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                     <  2:801534/146‑1 (MQ=37)
                                                                                                               ttGGGGGTTACAGGTAAATTCCGGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                                   <  2:554991/150‑1 (MQ=255)
                                                                                                                tGGGGGTTACAGTAAAATTCCAGGAAAAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCaa                                  <  1:886540/150‑1 (MQ=255)
                                                                                                                 ggggtTTACAGGAAAATTCCGGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCaa                                  <  2:396108/149‑1 (MQ=37)
                                                                                                                 gggggTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACt                                >  2:469503/1‑151 (MQ=255)
                                                                                                                 gggggTTACAGGTAAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACt                                <  1:664007/151‑1 (MQ=255)
                                                                                                                  gggtTTACAGGTAAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt                               <  2:55423/151‑1 (MQ=255)
                                                                                                                  ggggTTACAGGTAAATTCCAGGCAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt                               <  2:433845/151‑1 (MQ=255)
                                                                                                                    ggTTAAGGGTAATTCCCGGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                             <  1:383479/151‑1 (MQ=37)
                                                                                                                        aCAGTTAAATCCGAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGTCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                             <  2:38922/147‑1 (MQ=255)
                                                                                                                           ggaaaaTCCCAGGAAAAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg                          <  1:336530/144‑1 (MQ=255)
                                                                                                                           ggTAAATTCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCt                       <  1:182727/150‑1 (MQ=255)
                                                                                                                              aaatccgaggagaaaaaaac‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTAGCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCt                       <  1:119839/134‑1 (MQ=25)
                                                                                                                               aaTTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa                   >  2:388017/1‑150 (MQ=255)
                                                                                                                                 ttCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagc                >  2:1066469/1‑151 (MQ=255)
                                                                                                                                        cAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc         >  2:267099/1‑151 (MQ=255)
                                                                                                                                             aaaac‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAg    <  2:147303/151‑1 (MQ=255)
                                                                                                                                             aaaac‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAg    <  1:28232/151‑1 (MQ=255)
                                                                                                                                             aaaaa‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATc      >  2:989848/1‑149 (MQ=255)
                                                                                                                                                aa‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGc  >  2:191872/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGC  >  NC_000913/2001773‑2002066

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