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

AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGC  >  NC_000913/2001769‑2002066
                                                                                                                                                        |                                                                                                                                                    
aGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAAc                                                                                                                                                     <  1:1514555/151‑1 (MQ=255)
 gCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAAcc                                                                                                                                                    <  2:1026765/151‑1 (MQ=255)
 gCGGTATAGAGTCAATTCACTATAAACAGCCGTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCATGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTGCATGTAAAT‑TCCAG‑GCA‑GAAAAAAAc                                                                                                                                                     <  2:594663/150‑1 (MQ=255)
  cGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAAccc                                                                                                                                                   >  2:536926/1‑151 (MQ=255)
        agagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCgg                                                                                                                                             >  1:832121/1‑151 (MQ=255)
          agTGACTTCCCGATAAGCGGCCTTGTGTTATAGGGTTTTTCGGTATGTTTATCCGTTTCTGCAGGGTTTTTAAGCGGACGTTTAGTGGGTGAAAGGAGGGGTTATTTGGGGGTGACAGTTAATT‑TCCAG‑GAA‑AAAAAACACCCCGCCGGt                                                                                                                                            <  2:768853/150‑1 (MQ=255)
          agTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTg                                                                                                                                           >  1:96695/1‑151 (MQ=255)
            tGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGc                                                                                                                                         >  1:1550449/1‑151 (MQ=255)
             gAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCg                                                                                                                                        >  1:861267/1‑151 (MQ=255)
              atttCACAAAAAACAGCCCTGCGTTATATGGGTTATCGGCATGTTTACCCATTTCTGCAGGGTTTTTAATCGGGCGATTAGTGGGTGAAATGAGGGTTTATTTGGGGTTTACAGGTAAAT‑TCAAG‑GCA‑AAAAAAACCCCCGCCGGTGGc                                                                                                                                         <  2:789162/147‑1 (MQ=255)
                ttCAGGATAAACAGCCCGGGGTTATAGGATTTTTCGGCATGATTATCCTTTTCTGCGGGGTTTTTAATCGGACGATTAGGGGGTAAAAGGAGGGTTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAACCCCCGCCGGTGGCggg                                                                                                                                      <  1:1568890/150‑1 (MQ=255)
                ttCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCgggg                                                                                                                                     >  1:1269024/1‑151 (MQ=255)
                ttCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCgggg                                                                                                                                     >  2:1058796/1‑151 (MQ=255)
                     gATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGc                                                                                                                                 >  1:1338841/1‑150 (MQ=255)
                     gATAAACAGCCCTGCGTTATAGGAGTTATCGGAATGATTATCCGTTTCTGCGGGGTTTTAAAGGGGACGATTGGGGGGTAAAAGGGGGGTTTATTGGGGGGTTACAGGAAAAT‑CCAAG‑GCA‑GAAAAAACCCCCGCCGGTGGCGGGGaa                                                                                                                                   <  1:1469172/148‑1 (MQ=255)
                        aaaCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGt                                                                                                                             >  1:614267/1‑151 (MQ=255)
                           cAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCGGTGTTTTTAATCGGACGATTAGGGGGTAAAATGAGGGGTTTTTTGGGGTTTACAGGTAATT‑TCCAG‑GCA‑AAAAAAACCCCCGCCGGTGGCGTGGAAGCACGtt                                                                                                                            <  2:789170/149‑1 (MQ=255)
                            agcgatgagTTATATGAGTTATCGGCATGATTATTAGTTTCTGGAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTAAAGGTAAAT‑TCCAG‑GGA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTgat                                                                                                                         >  2:335391/9‑149 (MQ=255)
                               ccTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGAc                                                                                                                      >  1:356080/1‑151 (MQ=255)
                                 tGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                                     >  1:581387/1‑150 (MQ=255)
                                   cGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAt                                                                                                                  >  2:119428/1‑151 (MQ=255)
                                     ttgtatGATTTATGGGCATGATGTTCCTTTTGTGCAGGGTTTTTATTCGTTATTATAGGGATGGAATTGAGGGGTTATTGGGGGTTTAAGGATAATT‑CCCAG‑ACA‑GAAAAAAGCCCCGCCGGTGGCGGGGAAGCACGTTGCTg                                                                                                                        <  2:451040/140‑1 (MQ=255)
                                     ttattggtgTTATCGGAATGATTACCCGTTTCTGATGGGTTTTTAAGCGGACGATTGGTGGGTGAAAGGGGGGGTTTTTTGGGGGTTCAAGGAAAAT‑GCCGG‑GAA‑AAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaa                                                                                                                    <  2:1553490/139‑1 (MQ=255)
                                        tatGATTTCGCGGAATGATTTTCGGTTTCTGCAGGGTTTTAAATGGGACGATTAGTGGGTAAAATGAGGGTTTTTTTGGGGGTTACAGTTAAAT‑TCAAG‑GAA‑AAAAAAACCCCCGCCGGTGGCGGGGAAGCGCGTTGCTGACAAATTgc                                                                                                               <  2:905790/149‑1 (MQ=255)
                                         atGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCt                                                                                                            >  1:717437/1‑151 (MQ=255)
                                              ttgtggtaatatttATCCTTTTCGGCAGGTTTTTTAATCGGACTTGTGGGGGGTAAAATGAGGGTTTTTTGGGGGGTGAAAGTTACAT‑TCCAG‑AAA‑AAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTGAt                                                                                                        <  2:392902/138‑1 (MQ=255)
                                               tATCTGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt                                                                                                      >  2:1574689/1‑151 (MQ=255)
                                                    gCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg                                                                                                 >  2:994811/1‑151 (MQ=255)
                                                        gATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                             >  1:424360/1‑151 (MQ=255)
                                                          ttATCCGTTTCTGCGGGTTTTTTAACCGGAGAATTAGTGGGTGAAAGGGGGGTTTATTTGGGGGTAACAGTAATTT‑TCCAG‑GCA‑AAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGag                                                                                              <  2:232145/148‑2 (MQ=255)
                                                           tATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAACACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcgg                                                                                          >  2:658726/1‑151 (MQ=255)
                                                            aTCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGAGCTTTATGTTGCCGGATgcgg                                                                                          >  2:1126260/1‑150 (MQ=255)
                                                              ccGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGTTTATTTGGGGGTGACAGGTAAAT‑CCCGG‑GAA‑AAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTGTATGTTGCCGGATgcg                                                                                           <  2:72566/147‑1 (MQ=255)
                                                                 tttCTGCGGGGTTTTTAATCGGACGATTAGTGGGTGAAATGGGGGTTTATTTGGGGGTTACAGTAAAAT‑CCCGG‑GCA‑AAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaa                                                                                     <  1:396936/150‑1 (MQ=255)
                                                                    cTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGc                                                                                 >  1:1572777/1‑151 (MQ=255)
                                                                       cAGGGTTTTAAATCGAACGATGGGTGGGTAAAATGAGGGGTTTTTTGGGGGTTAGAGTTAAAT‑CCCAG‑GCA‑AAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaaa                                                                                    <  2:1180065/145‑1 (MQ=255)
                                                                        aGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTCTTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCAtt                                                                              >  2:1301062/1‑150 (MQ=255)
                                                                        aGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                              >  1:885908/1‑150 (MQ=255)
                                                                         ggttttttAATCGGACGATTAGGGGGTAAAATGGGGGTTTATTGGGGGGTTACAGGAAAAT‑CCCAG‑GCA‑AAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTa                                                                             <  1:36932/147‑1 (MQ=255)
                                                                         gggTTTTTAATCGGGCATTTGGTGGGTGAAAGGAGGGGTTTTTGGGGGGTTGAAGGAAATT‑TCCGG‑GAA‑AAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTGTATGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                            <  2:255684/151‑1 (MQ=255)
                                                                         gggTTTTCAATCGGACGGTTGGGGGGGAAAAGGGGGGGTTATTTGGGGGTAAGAGGAAATT‑TCCAG‑GCA‑GAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                              <  2:828702/149‑1 (MQ=255)
                                                                         gggTTTTAAACCGACCGTTGAGTGGGGGAAAGGAGGGGTTATTTGGGGGTTAAGGGTAATT‑CCCAG‑GCA‑GAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                              <  2:263294/149‑1 (MQ=255)
                                                                            ttttaAATCGGACGATTAGTGGGTAAAAGGAGGGTTTATTGGGGGGTTACAGATAAAT‑TCCAG‑GAA‑GAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                           <  1:21114/149‑1 (MQ=255)
                                                                             ttttAATCGGACGTGTAGTGGGTAAAATGGGGGGTTATTTGGGGGTTACAGTTAAAT‑TCCAG‑GCA‑GAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                         <  2:153283/150‑1 (MQ=39)
                                                                             ttttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                        >  2:117024/1‑151 (MQ=255)
                                                                              tttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                       >  2:1007437/1‑151 (MQ=255)
                                                                               ttAATCGGACGATTAGTGGGTAAAATGGGGGGTTATTTGGGGGTTAAAGGAAAAT‑CCCAG‑GCA‑AAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                       <  2:949144/150‑1 (MQ=38)
                                                                                tAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                     >  1:379533/1‑151 (MQ=255)
                                                                                 aaTCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                     >  1:1479165/1‑150 (MQ=255)
                                                                                 aaTCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GAA‑AAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                    <  1:201662/151‑1 (MQ=255)
                                                                                   tCGGACGATTAGGGGGTAAAATGGGGGGTTTTTGGGGGGTTACAGGAAAAT‑CCCAG‑CCA‑AAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                   <  1:595506/150‑1 (MQ=37)
                                                                                    cGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGCAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                                  >  1:473932/1‑150 (MQ=255)
                                                                                        ggATTAGTGGGTGAAA‑GAGGGGTTTTTTTGGGGGGTGAAGTTAAT‑TCCGGCGAA‑CAAAGATCCGCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTTCCGGATGCGTCGTAATCGCCTTATCCGGCCTACAAAAAt                                                             <  2:1633286/150‑1 (MQ=19)
                                                                                         gATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                            >  1:933191/1‑151 (MQ=255)
                                                                                          aTTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                            >  2:273131/1‑150 (MQ=255)
                                                                                           ttAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                                          >  2:530228/1‑151 (MQ=255)
                                                                                           tgaGTGGGTGAAGGGGGGGTTTATGTGGGGTTTACAGATAACT‑CCCAG‑CCA‑AAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAt                                                             <  1:204702/146‑1 (MQ=25)
                                                                                            tAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                         >  1:135743/1‑151 (MQ=255)
                                                                                             aGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                         >  2:34541/1‑150 (MQ=255)
                                                                                               tGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa                                                       >  1:1170167/1‑150 (MQ=255)
                                                                                                gggTAAAAGGGGGGTTTATTGGGGGGTTACAGGAAAAT‑CCCAG‑GAA‑AAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaaa                                                      <  1:1059640/150‑1 (MQ=37)
                                                                                                 gggaaaaggggggtTTATTTGGGGGTTACAGGTAAAT‑CCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAtt                                                    <  1:232689/141‑1 (MQ=255)
                                                                                                 ggTGAAATGGGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GAA‑AAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAtt                                                    <  1:546258/151‑1 (MQ=255)
                                                                                                 ggTAAAAGGGGGGGTTTTTTGGGGGTTAAAGGAAAAT‑TCCAG‑GAA‑AAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAtt                                                    <  2:503764/151‑1 (MQ=37)
                                                                                                   tGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTc                                                   >  2:1040645/1‑150 (MQ=255)
                                                                                                    gacaTGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGAACGTTGCTGACACATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAACTGTTCAAATTCa                                                  >  2:741014/4‑150 (MQ=255)
                                                                                                     aaaGGAGGGTTTATTGGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑AAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa                                                  <  1:1238308/149‑1 (MQ=255)
                                                                                                      aaTGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                                >  1:834029/1‑150 (MQ=255)
                                                                                                      aaTGAGGAGTTATTGGGGGGTTA‑AAGTAAATGTCCAG‑CAA‑AAAAAATTCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGTATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCaa                                                 <  2:1025150/149‑1 (MQ=25)
                                                                                                        tGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                              >  1:214194/1‑150 (MQ=255)
                                                                                                        tGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                              >  2:1032524/1‑150 (MQ=255)
                                                                                                           gggtTTATTTGGGGGTTCAAGAAAAAT‑CCCAG‑GAA‑AAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                          <  2:1092463/151‑1 (MQ=25)
                                                                                                            gggTTTTTTGGGGGTTAAAGGTAAAT‑ACCAG‑GAA‑AAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                                         <  2:1591631/151‑1 (MQ=37)
                                                                                                                   ttGGGGGTTACGGGTAATT‑TCCAG‑CAA‑AAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                                   <  2:1337888/150‑1 (MQ=37)
                                                                                                                     gggggTTGCAGGGAAGG‑GCGAG‑GAG‑GAAAAAAGCGCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGGCTTATCCGGCGTACAAAAATGTGCAAATGCAAGAAATTGCaggggaa                                  <  2:767499/149‑7 (MQ=21)
                                                                                                                     gggggTTACAGGT‑AAT‑TCCTG‑GGATGAAAAAATCCCCGCCGGTGGCGGGGAAGCACGTTCCTGACGATTCGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGAATACAAAAATGTGGAAATTAAATAAATTGCTATGCaa                                  <  2:657449/149‑1 (MQ=25)
                                                                                                                      ggggTTACAGGTAATT‑TCCAG‑GAA‑AAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt                               <  2:87866/151‑1 (MQ=255)
                                                                                                                            aCAGGAAAAT‑TCCAG‑GCA‑AAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAg                           <  2:581134/149‑1 (MQ=255)
                                                                                                                               ggTAAAT‑CCCGG‑GAA‑AAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg                          <  1:953453/147‑1 (MQ=255)
                                                                                                                                  aaaT‑CCCAG‑GCA‑AAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGa                     <  2:1066459/149‑1 (MQ=255)
                                                                                                                                    aT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGCTGCTGACAAATTCCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTCGTAGGCCTGATa                   >  2:1551141/1‑149 (MQ=255)
                                                                                                                                     aattcccagagaaaaaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg                          <  2:473912/131‑1 (MQ=37)
                                                                                                                                         cAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc              >  1:842116/1‑150 (MQ=255)
                                                                                                                                           g‑gCA‑AAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagc           <  1:991268/151‑1 (MQ=255)
                                                                                                                                             ggaaaaaaaacaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgca             <  2:918034/144‑1 (MQ=255)
                                                                                                                                             gacaaaaaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcag            <  1:1008249/145‑1 (MQ=255)
                                                                                                                                             gCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcg          >  1:743364/1‑151 (MQ=255)
                                                                                                                                                  aaaaaaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAgcgcagcgca        >  2:1412993/1‑149 (MQ=255)
                                                                                                                                                       aaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGCTGCGGCGTAAACTCCTTATCCAGCCTACAAAAATGTGCAACTTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGc  >  2:1349797/1‑150 (MQ=255)
                                                                                                                                                       aaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGc  >  2:1454611/1‑150 (MQ=255)
                                                                                                                                                       aaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGAATCAgg   >  2:762304/1‑149 (MQ=255)
                                                                                                                                                        |                                                                                                                                                    
AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAAT‑TCCAG‑GCA‑GAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGC  >  NC_000913/2001769‑2002066

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