Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,001,9180AC40.2% 172.0 / 20.9 93intergenic (‑129/+192)fliA/fliCRNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base A (47/8);  new base C (0/37);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 3.03e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGC  >  NC_000913/2001770‑2002066
                                                                                                                                                       |                                                                                                                                                     
gCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑cc                                                                                                                                                    <  1:554534/151‑1 (MQ=255)
  ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑cccc                                                                                                                                                  >  2:992582/1‑151 (MQ=255)
  ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑cccc                                                                                                                                                  <  2:772795/151‑1 (MQ=255)
  ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑cccc                                                                                                                                                  >  1:835369/1‑151 (MQ=255)
      taGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCccgccg                                                                                                                                              >  2:884192/1‑151 (MQ=255)
       agagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCgg                                                                                                                                             <  1:1457565/151‑1 (MQ=255)
         agTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTg                                                                                                                                           >  1:1386867/1‑151 (MQ=255)
             aaTTCACGATAAACAGCCCTGCTTTATATGATTTATGGGAATGATTATCCTTTTCTGCAGGTTTTTTAATCGGACGATTAGTGGGTAAAATGGGGGTTTATTTGGGGG‑TT‑AAAGGTAAATTCCAGGCAAAAAAAA‑C‑CCCCGCCGGTGGc                                                                                                                                         <  2:549481/149‑1 (MQ=255)
                 cACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGa                                                                                                                                    >  1:172310/1‑150 (MQ=255)
                   cGAAAAACACCCTTGTTTTATATGATTTACCGGCTTGTTTATCCTTTTCTGCGGGTTTTTTAATCGGACGATTAGTGGGTAAAATGGGGGTTTATTTGGGGT‑TT‑ACAGGAAAATTCCGGGCAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGc                                                                                                                                 <  1:443359/151‑1 (MQ=255)
                       aaaCAGCCCTGCGTTATAGGAGTTATCGGCTTTTTTTTCCGTTTATGCAGTTTTTTTACTGGTTCTTTTAGTGGTAAAATTGTGGGGTTTTTTGGGTG‑TG‑ACGGGAAATTTCCAGCCAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCa                                                                                                                                <  2:952907/148‑1 (MQ=255)
                       aaaCACCCCTGCGTTATATGATTTATCGGCTTGATTATCCGTTTCTGCGGGGTTTTTAATCGAACGATTAGGGGGTAAAAGGAGGGGTTATTTGGGGG‑TT‑AAAGTAAAATCCCAGGAAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCAc                                                                                                                               <  2:87245/149‑1 (MQ=255)
                           aGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCt                                                                                                                         >  1:1201443/1‑151 (MQ=255)
                            gCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTg                                                                                                                        >  1:1011225/1‑151 (MQ=255)
                              ccTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGAc                                                                                                                      >  1:821729/1‑151 (MQ=255)
                                tGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                                     >  1:300498/1‑150 (MQ=255)
                                   gTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                                 >  1:666014/1‑151 (MQ=255)
                                    ttATATGAGTTATCGGCATGTTTACCCGTTTCTGCGGGGTTTTTAATGGGACGATTGGTGGGTGAAATGAGGGTTTATTGGGGGG‑TT‑AAAGGAAAATCCCAGGCAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                                 <  1:582971/150‑1 (MQ=255)
                                    ttATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg                                                                                                                >  1:411609/1‑151 (MQ=255)
                                      atatGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcg                                                                                                              >  1:1555619/1‑151 (MQ=255)
                                             ttATGGGCATGATTATCCGTTTCTGCGGGTTTTTTAATCGGACGATTAGTGGGTAAAATGAGGGGTTTTTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt                                                                                                          <  1:1564750/148‑1 (MQ=255)
                                             ttATCGGCTTGATTATCCTTTTCTGCAGGTTTTTTAATCGGACGTTTAGTGGGTGAAAGGGGGGTTTATTGGGGGG‑TT‑ACAGGAAAATTCCAGGCAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcg                                                                                                              <  2:156213/144‑1 (MQ=255)
                                              tATCGGATTGTTTATCCGTTTCTGCGGGTTTTTTAATCGGACGATTAGGGGGTAAAATG‑GGGGTTTTTTGGGGGTTT‑ACAGGAAAATCCCAGGAAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATg                                                                                                       <  1:134192/150‑1 (MQ=255)
                                               atgggaaTGATTATCCGTTTCGGCAGGTTTTTTAATCGGAGGATTAGGGGGGAAAAGGGGGGTTTATTGGGGGT‑TT‑ACGGATAAATTCCAGCAAAAAAAAC‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt                                                                                                          <  2:78930/140‑1 (MQ=255)
                                                    cATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                                >  2:57812/1‑151 (MQ=255)
                                                        aTTATCTGTTTCTGCGGGGTTTTAAATCGGACGATTAGGGGGTAAAATGGGGGGTTATTGGGGGG‑TT‑AAAGGTAAATTCCAGGCAAAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                             <  1:1239211/150‑1 (MQ=255)
                                                         ttccccGTTTTTGAGGGTTTTTTAGTCGGTGATTTAGTGGGTCAAATGAGGGTTCAGTGGGGGC‑TA‑ACGGCTAATTTCCAGCAAGAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt                                                                                              <  1:509902/144‑1 (MQ=255)
                                                         ttATCCTTTTCGGCAGGTTTTTTAATCGGACGATTAGGGGGTAAAAGGGGGGTTTATTGGGGGT‑TT‑ACAGGTAAATTCCAGGAAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGa                                                                                               <  1:1130764/147‑1 (MQ=255)
                                                          tATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcgg                                                                                          >  1:755183/1‑151 (MQ=255)
                                                          tATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcgg                                                                                          >  2:715812/1‑151 (MQ=255)
                                                           aTCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggc                                                                                         >  1:1393449/1‑151 (MQ=255)
                                                            tCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggcg                                                                                        >  1:721190/1‑151 (MQ=255)
                                                              cGTTTCTGGAGGGTTTTTAAGCGGACGATTTGTGGGTGAATGGGGGGGTTTTTTGGGGT‑TT‑AAAGGTAAATACCAGGAAGAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                             <  2:635481/144‑1 (MQ=255)
                                                                tttCTGCAGGTTTTTTAATCGGACGATTAGTGGGTGAAAGGGGGGTTTTTTGGGGGG‑TT‑ACAGGAAAATCCCAGGAAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaa                                                                                     <  2:647643/150‑1 (MQ=255)
                                                                     gCAGGTTTTTTAATCGGACGATTGGGGGGTGAAAGGAGGGGTTATTGGGGGG‑TT‑ACAGTTAAATTCCAGGCAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCt                                                                               <  2:128214/151‑1 (MQ=255)
                                                                     gCAGGGTTTTAAATCGGACGATTGGGGGGTAAAATGAGGGTTTATTGGGGGG‑TT‑AAAGGTAAATCCCAGGAAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACg                                                                                  <  1:534550/148‑1 (MQ=255)
                                                                     gCAGGCATTTTAATCGACCGATTCGAGTGCGATATGTCGAGTTGTTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGCGACAGGACGTTGCTGACAAATTGCGCTTTATGTTTCCGGATGCGGCGTTAACg                                                                                  >  2:319163/1‑148 (MQ=255)
                                                                       aGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                              >  1:1240483/1‑150 (MQ=255)
                                                                        ggttttttATTGGGACGATTGGGGGGTAAAATGAGGGTTTATTTGGGTT‑TT‑ACAGAAAATTTCGAGGAAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGc                                                                                 <  1:1257184/143‑1 (MQ=255)
                                                                        ggttttttAATCGGACGATTAGTGGGTGAAATGAGGGTTTATTTGGGGG‑TT‑ACAGGAAAATTCCAGGCAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACg                                                                                  <  1:683654/142‑1 (MQ=255)
                                                                        gggTTTTTATCCGGATAATGAGTGGGGGAAAGGAGGGGTTATTTGGGGT‑TT‑ACGGGAAAATCCAAGGCAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                                   <  1:294418/144‑1 (MQ=255)
                                                                        gggTTTTTAATCGGACGATTAGTGGGTGAAATGCGGGGTTCTTTGGGGG‑TT‑ACAGGTAAATTACAGGCTGAAAAAA‑A‑CCCCGCCTTTGGCGGGGAAGCACGTTTCTGACAAATTACGCTTCATGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                            >  2:873117/1‑151 (MQ=255)
                                                                         ggTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                           >  1:75801/1‑151 (MQ=255)
                                                                         ggTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                           >  1:1104828/1‑151 (MQ=255)
                                                                           tttttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                         >  2:1416853/1‑151 (MQ=255)
                                                                           tttttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                         >  1:1099963/1‑151 (MQ=255)
                                                                            ttttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                        >  1:904394/1‑151 (MQ=255)
                                                                               tAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                     >  1:1538029/1‑151 (MQ=255)
                                                                               tAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                     >  2:509679/1‑151 (MQ=255)
                                                                                aaTCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                     >  1:312585/1‑150 (MQ=255)
                                                                                   cGGACGTTTATGGGGTAAAATGAGGGTTTTTTTGGGGG‑TT‑ACAGGTAAATCCCAGGCAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                   <  2:673153/149‑1 (MQ=38)
                                                                                    ggACGATTGGGGGGTAAAATG‑GGGGTTTTTTGGGGGTTT‑ACAGGAAAATCCCAGGCAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaa                                                                <  1:309434/151‑1 (MQ=37)
                                                                                    ggACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaa                                                                 >  1:646219/1‑150 (MQ=255)
                                                                                     gACGATTAGTGGGAGAAACGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAACA‑A‑CCCCGCCGGTGGCGGGGAAGCACGAAGCTGACAAATTGCGCTTTAAGTTGCCGGCGGCGGCGTACACACCCTATCTGGCCGACaaa                                                                >  2:109716/1‑150 (MQ=255)
                                                                                          ttggggggaaaaaTGGGGGGTTTTTTGGGGT‑TT‑ACAGGAAAATCCCAGGAAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaaa                                                              <  1:660530/137‑1 (MQ=25)
                                                                                          ttggggggTAAAAGGAGGGGTTATTGGGGGG‑TT‑ACAGGTAAATTCCAGGAAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                            <  2:985929/144‑1 (MQ=38)
                                                                                          ttAGTGGGTGAAAGGGGGGTTTATTTGGGGG‑TT‑ACAGGTAAATCCCAGGAAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgt                                                           <  1:1238797/150‑1 (MQ=38)
                                                                                             gTGGGAGAAATGGGGGGTTTTTTGGGGT‑TA‑ACGGGAAAATGCGAGGAAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgt                                                           <  1:178973/147‑1 (MQ=25)
                                                                                              tGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa                                                       >  2:945406/1‑150 (MQ=255)
                                                                                               gggTGAAAGGGGGGTTTTTTTGGGGT‑TT‑ACAGTTAAATTCCGGGAAAAAAAAA‑T‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaaa                                                      <  2:1234583/150‑1 (MQ=37)
                                                                                               ggaaaaaaTGAGGGGTGAATGGGGGG‑TT‑AAATAAAACTCCCAGTAAAAAAAAA‑A‑CACCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGGTGCGGCATAAACGCCTTATCCCGCCTACAAAAATGTGCaaa                                                      <  2:1204482/145‑1 (MQ=255)
                                                                                                      aTGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATa                                               >  1:801162/1‑150 (MQ=255)
                                                                                                      aTGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATa                                               >  2:677910/1‑150 (MQ=255)
                                                                                                        gggggtTTATTGGGGGT‑TT‑AAAGGTAAATTCCAGGCAGAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                              <  1:434839/147‑1 (MQ=37)
                                                                                                        gAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt                                            >  1:646871/1‑151 (MQ=255)
                                                                                                        gAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt                                            >  2:965496/1‑151 (MQ=255)
                                                                                                        gAGGGGTTATTGGGGGT‑TT‑ACAGGAAAATTCCAGGAAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaaa                                             <  1:1513624/150‑1 (MQ=37)
                                                                                                            gggtttttGGGGG‑TTAACAGGTAAATTCAAGGAAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                                         <  2:1201427/146‑1 (MQ=25)
                                                                                                             gTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCa                                        >  1:1372126/1‑150 (MQ=255)
                                                                                                              ttATTTGGGGG‑TT‑AGAGGTAAATTCCAGAAAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                              <  2:96459/143‑1 (MQ=255)
                                                                                                              ttATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAt                                      >  1:1539435/1‑151 (MQ=255)
                                                                                                               ttttttGGGG‑GT‑AAAGGTAAATTCCGGGGAAAAAAAACC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                                         <  2:858704/146‑1 (MQ=25)
                                                                                                               tATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                     >  2:612940/1‑151 (MQ=255)
                                                                                                                aTTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                                    >  2:106080/1‑151 (MQ=255)
                                                                                                                aTTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                                    >  2:217827/1‑151 (MQ=255)
                                                                                                                 ttgggggg‑TT‑AAAGTAAAATTCCAGGAAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                                    <  1:472559/147‑1 (MQ=37)
                                                                                                                    ggggg‑TT‑AAAGATAAATTCAAGCCAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                     <  1:558641/146‑1 (MQ=37)
                                                                                                                     gggg‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt                               >  1:373187/1‑151 (MQ=255)
                                                                                                                      ggggataacggttAAATTCAAGGAAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCaa                                  <  1:1108713/137‑1 (MQ=25)
                                                                                                                      ggg‑TT‑ACAGGTAAATTCCAGGCAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg                              <  1:1214048/151‑1 (MQ=255)
                                                                                                                       gg‑TT‑ACAGGAAAATCCCAGGAAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                             <  1:367868/151‑1 (MQ=255)
                                                                                                                             aCAGGTAAATCCCAGGCAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGc                         <  1:371481/151‑1 (MQ=255)
                                                                                                                              cAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGcc                        >  2:4494/1‑151 (MQ=255)
                                                                                                                              acggataaTTCCCAGCAAAAAAAA‑ACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGc                         <  1:1304215/145‑1 (MQ=25)
                                                                                                                                   aaaTTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt                    >  1:198176/1‑150 (MQ=255)
                                                                                                                                   aaaTCCCAGGCAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa                   <  2:789152/151‑1 (MQ=255)
                                                                                                                                   aaaTCCCAGGAAAAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCCACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGcc                        <  1:1088386/146‑1 (MQ=255)
                                                                                                                                    aaTTCCAGGCAGAAAAAA‑C‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATaa                  <  1:205458/151‑1 (MQ=255)
                                                                                                                                    aaTTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa                   >  1:941893/1‑150 (MQ=255)
                                                                                                                                        ccAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc              >  1:839549/1‑151 (MQ=255)
                                                                                                                                          aGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc              >  2:692611/1‑149 (MQ=255)
                                                                                                                                          aGGCAAAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgca             <  2:1172899/150‑1 (MQ=255)
                                                                                                                                           ggaaaaaaaaa‑c‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgca             <  1:438766/144‑1 (MQ=255)
                                                                                                                                               gAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg                          >  2:5929/1‑132 (MQ=255)
                                                                                                                                               aaaaaaa‑a‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAt       <  1:1282711/150‑1 (MQ=255)
                                                                                                                                                 aaaaa‑a‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATc      >  1:1157872/1‑150 (MQ=255)
                                                                                                                                                 aaaaa‑a‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATc      >  2:33554/1‑150 (MQ=255)
                                                                                                                                                   aaa‑a‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAgg   >  1:1087795/1‑151 (MQ=255)
                                                                                                                                                     a‑a‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGc  >  1:1283940/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
GCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGGTAAATTCCAGGCAGAAAAAA‑A‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGC  >  NC_000913/2001770‑2002066

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: ‑