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

AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGC  >  NC_000913/2001769‑2002066
                                                                                                                                                        |                                                                                                                                                    
aGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAc                                                                                                                                                     <  1:1514555/151‑1 (MQ=255)
 gCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAcc                                                                                                                                                    <  2:1026765/151‑1 (MQ=255)
 gCGGTATAGAGTCAATTCACTATAAACAGCCGTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCATGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTGCATGTAAATTCCAGGCAGAAAAAAAc                                                                                                                                                     <  2:594663/150‑1 (MQ=255)
  cGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAccc                                                                                                                                                   >  2:536926/1‑151 (MQ=255)
        agagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCgg                                                                                                                                             >  1:832121/1‑151 (MQ=255)
          agTGACTTCCCGATAAGCGGCCTTGTGTTATAGGGTTTTTCGGTATGTTTATCCGTTTCTGCAGGGTTTTTAAGCGGACGTTTAGTGGGTGAAAGGAGGGGTTATTTGGGGGTGACAGTTAATTTCCAGGAAAAAAAACACCCCGCCGGt                                                                                                                                            <  2:768853/150‑1 (MQ=255)
          agTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTg                                                                                                                                           >  1:96695/1‑151 (MQ=255)
            tGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGc                                                                                                                                         >  1:1550449/1‑151 (MQ=255)
             gAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCg                                                                                                                                        >  1:861267/1‑151 (MQ=255)
              atttCACAAAAAACAGCCCTGCGTTATATGGGTTATCGGCATGTTTACCCATTTCTGCAGGGTTTTTAATCGGGCGATTAGTGGGTGAAATGAGGGTTTATTTGGGGTTTACAGGTAAATTCAAGGCAAAAAAAACCCCCGCCGGTGGc                                                                                                                                         <  2:789162/147‑1 (MQ=255)
                ttCAGGATAAACAGCCCGGGGTTATAGGATTTTTCGGCATGATTATCCTTTTCTGCGGGGTTTTTAATCGGACGATTAGGGGGTAAAAGGAGGGTTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCggg                                                                                                                                      <  1:1568890/150‑1 (MQ=255)
                ttCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCgggg                                                                                                                                     >  1:1269024/1‑151 (MQ=255)
                ttCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCgggg                                                                                                                                     >  2:1058796/1‑151 (MQ=255)
                     gATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGc                                                                                                                                 >  1:1338841/1‑150 (MQ=255)
                     gATAAACAGCCCTGCGTTATAGGAGTTATCGGAATGATTATCCGTTTCTGCGGGGTTTTAAAGGGGACGATTGGGGGGTAAAAGGGGGGTTTATTGGGGGGTTACAGGAAAATCCAAGGCAGAAAAAACCCCCGCCGGTGGCGGGGaa                                                                                                                                   <  1:1469172/148‑1 (MQ=255)
                        aaaCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGt                                                                                                                             >  1:614267/1‑151 (MQ=255)
                           cAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCGGTGTTTTTAATCGGACGATTAGGGGGTAAAATGAGGGGTTTTTTGGGGTTTACAGGTAATTTCCAGGCAAAAAAAACCCCCGCCGGTGGCGTGGAAGCACGtt                                                                                                                            <  2:789170/149‑1 (MQ=255)
                            agcgatgagTTATATGAGTTATCGGCATGATTATTAGTTTCTGGAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTAAAGGTAAATTCCAGGGAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTgat                                                                                                                         >  2:335391/9‑149 (MQ=255)
                               ccTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGAc                                                                                                                      >  1:356080/1‑151 (MQ=255)
                                 tGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                                     >  1:581387/1‑150 (MQ=255)
                                   cGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAt                                                                                                                  >  2:119428/1‑151 (MQ=255)
                                     ttgtatGATTTATGGGCATGATGTTCCTTTTGTGCAGGGTTTTTATTCGTTATTATAGGGATGGAATTGAGGGGTTATTGGGGGTTTAAGGATAATTCCCAGACAGAAAAAAGCCCCGCCGGTGGCGGGGAAGCACGTTGCTg                                                                                                                        <  2:451040/140‑1 (MQ=255)
                                     ttattggtgTTATCGGAATGATTACCCGTTTCTGATGGGTTTTTAAGCGGACGATTGGTGGGTGAAAGGGGGGGTTTTTTGGGGGTTCAAGGAAAATGCCGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaa                                                                                                                    <  2:1553490/139‑1 (MQ=255)
                                        tatGATTTCGCGGAATGATTTTCGGTTTCTGCAGGGTTTTAAATGGGACGATTAGTGGGTAAAATGAGGGTTTTTTTGGGGGTTACAGTTAAATTCAAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCGCGTTGCTGACAAATTgc                                                                                                               <  2:905790/149‑1 (MQ=255)
                                         atGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCt                                                                                                            >  1:717437/1‑151 (MQ=255)
                                              ttgtggtaatatttATCCTTTTCGGCAGGTTTTTTAATCGGACTTGTGGGGGGTAAAATGAGGGTTTTTTGGGGGGTGAAAGTTACATTCCAGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTGAt                                                                                                        <  2:392902/138‑1 (MQ=255)
                                               tATCTGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt                                                                                                      >  2:1574689/1‑151 (MQ=255)
                                                    gCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg                                                                                                 >  2:994811/1‑151 (MQ=255)
                                                        gATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                             >  1:424360/1‑151 (MQ=255)
                                                          ttATCCGTTTCTGCGGGTTTTTTAACCGGAGAATTAGTGGGTGAAAGGGGGGTTTATTTGGGGGTAACAGTAATTTTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGag                                                                                              <  2:232145/148‑2 (MQ=255)
                                                           tATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAACACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcgg                                                                                          >  2:658726/1‑151 (MQ=255)
                                                            aTCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGAGCTTTATGTTGCCGGATgcgg                                                                                          >  2:1126260/1‑150 (MQ=255)
                                                              ccGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGTTTATTTGGGGGTGACAGGTAAATCCCGGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTGTATGTTGCCGGATgcg                                                                                           <  2:72566/147‑1 (MQ=255)
                                                                 tttCTGCGGGGTTTTTAATCGGACGATTAGTGGGTGAAATGGGGGTTTATTTGGGGGTTACAGTAAAATCCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaa                                                                                     <  1:396936/150‑1 (MQ=255)
                                                                    cTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGc                                                                                 >  1:1572777/1‑151 (MQ=255)
                                                                       cAGGGTTTTAAATCGAACGATGGGTGGGTAAAATGAGGGGTTTTTTGGGGGTTAGAGTTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaaa                                                                                    <  2:1180065/145‑1 (MQ=255)
                                                                        aGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTCTTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCAtt                                                                              >  2:1301062/1‑150 (MQ=255)
                                                                        aGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                              >  1:885908/1‑150 (MQ=255)
                                                                         ggttttttAATCGGACGATTAGGGGGTAAAATGGGGGTTTATTGGGGGGTTACAGGAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTa                                                                             <  1:36932/147‑1 (MQ=255)
                                                                         gggTTTTTAATCGGGCATTTGGTGGGTGAAAGGAGGGGTTTTTGGGGGGTTGAAGGAAATTTCCGGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTGTATGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                            <  2:255684/151‑1 (MQ=255)
                                                                         gggTTTTCAATCGGACGGTTGGGGGGGAAAAGGGGGGGTTATTTGGGGGTAAGAGGAAATTTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                              <  2:828702/149‑1 (MQ=255)
                                                                         gggTTTTAAACCGACCGTTGAGTGGGGGAAAGGAGGGGTTATTTGGGGGTTAAGGGTAATTCCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                              <  2:263294/149‑1 (MQ=255)
                                                                            ttttaAATCGGACGATTAGTGGGTAAAAGGAGGGTTTATTGGGGGGTTACAGATAAATTCCAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                           <  1:21114/149‑1 (MQ=255)
                                                                             ttttAATCGGACGTGTAGTGGGTAAAATGGGGGGTTATTTGGGGGTTACAGTTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                         <  2:153283/150‑1 (MQ=39)
                                                                             ttttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                        >  2:117024/1‑151 (MQ=255)
                                                                              tttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                       >  2:1007437/1‑151 (MQ=255)
                                                                               ttAATCGGACGATTAGTGGGTAAAATGGGGGGTTATTTGGGGGTTAAAGGAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                       <  2:949144/150‑1 (MQ=38)
                                                                                tAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                     >  1:379533/1‑151 (MQ=255)
                                                                                 aaTCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                     >  1:1479165/1‑150 (MQ=255)
                                                                                 aaTCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                    <  1:201662/151‑1 (MQ=255)
                                                                                   tCGGACGATTAGGGGGTAAAATGGGGGGTTTTTGGGGGGTTACAGGAAAATCCCAGCCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                   <  1:595506/150‑1 (MQ=37)
                                                                                    cGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGCAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                                  >  1:473932/1‑150 (MQ=255)
                                                                                        ggATTAGTGGGTGAAAGAGGGGTTTTTTTGGGGGGTGAAGTTAATTCCGGCGAACAAAGATCCGCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTTCCGGATGCGTCGTAATCGCCTTATCCGGCCTACAAAAAt                                                             <  2:1633286/150‑1 (MQ=19)
                                                                                         gATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                            >  1:933191/1‑151 (MQ=255)
                                                                                          aTTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                            >  2:273131/1‑150 (MQ=255)
                                                                                           ttAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                                          >  2:530228/1‑151 (MQ=255)
                                                                                           tgaGTGGGTGAAGGGGGGGTTTATGTGGGGTTTACAGATAACTCCCAGCCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAt                                                             <  1:204702/146‑1 (MQ=25)
                                                                                            tAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                         >  1:135743/1‑151 (MQ=255)
                                                                                             aGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                         >  2:34541/1‑150 (MQ=255)
                                                                                               tGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa                                                       >  1:1170167/1‑150 (MQ=255)
                                                                                                gggTAAAAGGGGGGTTTATTGGGGGGTTACAGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaaa                                                      <  1:1059640/150‑1 (MQ=37)
                                                                                                 gggaaaaggggggtTTATTTGGGGGTTACAGGTAAATCCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAtt                                                    <  1:232689/141‑1 (MQ=255)
                                                                                                 ggTGAAATGGGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAtt                                                    <  1:546258/151‑1 (MQ=255)
                                                                                                 ggTAAAAGGGGGGGTTTTTTGGGGGTTAAAGGAAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAtt                                                    <  2:503764/151‑1 (MQ=37)
                                                                                                   tGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTc                                                   >  2:1040645/1‑150 (MQ=255)
                                                                                                    gacaTGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGAACGTTGCTGACACATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAACTGTTCAAATTCa                                                  >  2:741014/4‑150 (MQ=255)
                                                                                                     aaaGGAGGGTTTATTGGGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa                                                  <  1:1238308/149‑1 (MQ=255)
                                                                                                      aaTGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                                >  1:834029/1‑150 (MQ=255)
                                                                                                      aaTGAGGAGTTATTGGGGGGTTAAAGTAAATGTCCAGCAAAAAAAATTCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGTATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCaa                                                 <  2:1025150/149‑1 (MQ=25)
                                                                                                        tGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                              >  1:214194/1‑150 (MQ=255)
                                                                                                        tGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                              >  2:1032524/1‑150 (MQ=255)
                                                                                                           gggtTTATTTGGGGGTTCAAGAAAAATCCCAGGAAAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                          <  2:1092463/151‑1 (MQ=25)
                                                                                                            gggTTTTTTGGGGGTTAAAGGTAAATACCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                                         <  2:1591631/151‑1 (MQ=37)
                                                                                                                   ttGGGGGTTACGGGTAATTTCCAGCAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                                   <  2:1337888/150‑1 (MQ=37)
                                                                                                                     gggggTTGCAGGGAAGGGCGAGGAGGAAAAAAGCGCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGGCTTATCCGGCGTACAAAAATGTGCAAATGCAAGAAATTGCaggggaa                                  <  2:767499/149‑7 (MQ=21)
                                                                                                                     gggggTTACAGGTAATTCCTGGGATGAAAAAATCCCCGCCGGTGGCGGGGAAGCACGTTCCTGACGATTCGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGAATACAAAAATGTGGAAATTAAATAAATTGCTATGCaa                                  <  2:657449/149‑1 (MQ=25)
                                                                                                                      ggggTTACAGGTAATTTCCAGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt                               <  2:87866/151‑1 (MQ=255)
                                                                                                                            aCAGGAAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAg                           <  2:581134/149‑1 (MQ=255)
                                                                                                                               ggTAAATCCCGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg                          <  1:953453/147‑1 (MQ=255)
                                                                                                                                  aaaTCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGa                     <  2:1066459/149‑1 (MQ=255)
                                                                                                                                    aTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGCTGCTGACAAATTCCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTCGTAGGCCTGATa                   >  2:1551141/1‑149 (MQ=255)
                                                                                                                                     aattcccagagaaaaaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg                          <  2:473912/131‑1 (MQ=37)
                                                                                                                                         cAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc              >  1:842116/1‑150 (MQ=255)
                                                                                                                                           ggCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagc           <  1:991268/151‑1 (MQ=255)
                                                                                                                                             ggaaaaaaaacaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgca             <  2:918034/144‑1 (MQ=255)
                                                                                                                                             gacaaaaaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcag            <  1:1008249/145‑1 (MQ=255)
                                                                                                                                             gCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcg          >  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)
                                                                                                                                                        |                                                                                                                                                    
AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGC  >  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: