Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,385,516 A→T 10.2% intergenic (‑129/+301) tldD ← / ← yhdP putative peptidase/DUF3971‑AsmA2 domains protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,385,5160AT10.2% 686.6 / 12.9 288intergenic (‑129/+301)tldD/yhdPputative peptidase/DUF3971‑AsmA2 domains protein
Reads supporting (aligned to +/- strand):  ref base A (161/76);  new base T (2/32);  total (163/125)
Fisher's exact test for biased strand distribution p-value = 1.77e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 289 total aligned reads displayed.

CTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC  >  CP009273/3385368‑3385662
                                                                                                                                                    |                                                                                                                                                  
 ttACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAg                                                                                                                                                   >  1:1396848/1‑149 (MQ=255)
  tACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAgg                                                                                                                                                  >  2:1272602/1‑149 (MQ=255)
   aCCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGc                                                                                                                                                 <  1:1797259/149‑1 (MQ=255)
     cAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCaa                                                                                                                                               <  1:4839770/149‑1 (MQ=255)
            aGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGt                                                                                                                                        >  1:863090/1‑149 (MQ=255)
             gACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCATCCTACGGTTATGTTCAGGATAGTAGGCCTGATGAGACGAACAGCGTCGCATAAGGCCACGTCTGTc                                                                                                                                       >  2:1926113/1‑149 (MQ=25)
             gACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTc                                                                                                                                       >  2:4400821/1‑149 (MQ=255)
              aCTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCg                                                                                                                                      >  1:1275774/1‑149 (MQ=255)
                tCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGa                                                                                                                                    >  1:4137448/1‑149 (MQ=255)
                 cATCTTTTTTGATACTCGTTAGTTACTGAAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGAt                                                                                                                                   <  2:2954774/149‑1 (MQ=255)
                  aTCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATg                                                                                                                                  >  2:2181640/1‑149 (MQ=255)
                   tCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATgc                                                                                                                                 >  1:3180681/1‑149 (MQ=255)
                      tttttGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATgcggc                                                                                                                              >  1:2222374/1‑149 (MQ=255)
                         ttGCTACTCGTTTGTTACTGCTGTCGTGGATTTTTTCCGGCTGCCGGAGGGGGCGTGACCGCCCTATCCAGCCTAGGTTTTTGTTCCGGTTTGTCGGCCTTATAAGACGCACAGCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTa                                                                                                                           <  1:3421562/149‑1 (MQ=18)
                         ttGCTACTCGTTAGTTACTGCAGTAGCGGTTTTTTTCCGGCTGCCGGTTGCGGCGTGACCGCCCTCTCCAGCCTACGGTTATGTTCCGGTTTGTCGGCCTGATACGCCGCACCGCGGCGCATCAGGCAACGGCTGTCGGATGCGGCGTa                                                                                                                           <  2:2355135/149‑1 (MQ=37)
                           gCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGATAGTAGGCCTTATAAAACGAAAATCGTCGAAACAAACAACGGCTGTCAGATGCGGCGTaaa                                                                                                                         >  2:4055180/1‑149 (MQ=25)
                           gCTACTCGTGAGTTCCTGCAGCCGAGGATTTTGTACGGCTGCCGGTTGCGGCGTGACCGCCCTATCCAGCCTACGGTTATGTTCCGGTTGGTAGGCCTGATAAGACGCACCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTaaa                                                                                                                         <  2:4472391/149‑1 (MQ=18)
                            cTACTCTTTTGTTACTGCATTTGTGTTTTTTTTCCGGCGGCCGGAGGCGGCGTGACCGCCCTATCCAGCCTACGGTTTTGTTCCGGTTTGTAGGCCTTATAAGACGCACCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAAc                                                                                                                        <  2:4960547/149‑1 (MQ=17)
                               cTCGTTAGTTACTGCAGTAGAGTATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGTTTATGTTCCGGTTTGTCGGCCTGATAAGACGCACCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGcc                                                                                                                     <  2:4604117/149‑1 (MQ=255)
                                tctttagttaCTGCAGTCGAGTATTTTTTCCGGCTGCCGGAGGCGGCGTGAACGCCCTACCCAGCATACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGAATCAGGCAACGGCTGTCGTATGCGGCGTAAACGCCt                                                                                                                    <  2:572704/146‑1 (MQ=18)
                                  gttagttaCTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTa                                                                                                                  <  2:814582/149‑1 (MQ=255)
                                   tttgttaCTGCCTTCGGGGTTTTTTTCGGGCTGCCGGTGGCGGCGTGACCGCCATACCCAGCCTACGTTTATGTTCGGTTTTGTCGGCCTGATAAGACGCACCGCGTCGCAACCGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTAt                                                                                                                 <  1:5177291/146‑1 (MQ=9)
                                   ttagttaCTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTAt                                                                                                                 <  2:3349783/149‑1 (MQ=255)
                                      gttaCTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATACTACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATACg                                                                                                              >  1:1131810/1‑149 (MQ=255)
                                      gttaCTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCg                                                                                                              >  1:734996/1‑149 (MQ=255)
                                         aCTGCAGTAGTGTTTTTTTTACGGCTGCCGGTTGCGGCGTGAACGCCCTAGCCAGCCTACGTTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGAcc                                                                                                           <  2:4397583/149‑1 (MQ=18)
                                         aCTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCATGAAAAGACGCACAGCGACGCAAAAGCCAACGGCTCTCGGATGCGGCGTAAACGCCTTATCCGAcc                                                                                                           >  1:3509088/1‑149 (MQ=18)
                                           tGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGACGGATGCGGCGTAAACGCCGTATCCGACCTa                                                                                                         >  1:2512739/1‑149 (MQ=255)
                                           tGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGCCTGACGGATGCGGCGTAAACGCCATATCCGACCTa                                                                                                         >  1:3641049/1‑149 (MQ=255)
                                           tGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCAGCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTTTCAGACCTa                                                                                                         >  2:2580417/1‑149 (MQ=255)
                                            gCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCGGGATAGTTGGCCTGTGCCGACGCACATCGTCGCATCAGGCAACGCCTCTCTGATCCGGAGTAAACGACTTAGCAGACCTAc                                                                                                        >  1:5125590/1‑149 (MQ=21)
                                            gCAGTAGAGGATTTTTTACGGCTGCCGGAAGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGATTGTAGGCCTGATAAGACGCAAATCGTCGCATAAGGCAACGGCTGTCGGATGCTGCGTAAACGCCTGATCCGACCTAc                                                                                                        >  1:185844/1‑149 (MQ=18)
                                              aGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGACGGATGCGGCGTAAACGCCTTATCCGACCTACgg                                                                                                      >  1:3720010/1‑149 (MQ=255)
                                                ttgtgtaTTTTTTACGGCTGCCGGATGCGGCGTGAACGCCATATCCAGCCTAGGGTTTTGTTCCGGTTTGTTGGCCTGATAAGACGCACCGCGTCGAATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGtt                                                                                                    <  2:553808/143‑1 (MQ=17)
                                                  gtggtttttttACGGCTGCCGGAGGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTTGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTAt                                                                                                  <  2:3100531/144‑1 (MQ=255)
                                                  gagGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGAAAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTAt                                                                                                  >  1:4380521/1‑149 (MQ=255)
                                                    ggATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGt                                                                                                >  1:2337046/1‑149 (MQ=255)
                                                     gtttttttACGGCTGCGGGTTGCGGCGTGAACGCCCTACCCAGCCTAGGTTTTTGTTCGGGTTTGTTGGCCTTATAAGACGCACCGCGTCGCATCCGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGtt                                                                                               <  2:476465/147‑1 (MQ=17)
                                                     gATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCATGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCTACCTACGGTTATGtt                                                                                               <  2:1688315/149‑1 (MQ=255)
                                                      tttttttCCGGCGGCCGGGTGCGGCGTGAACGCCACACCCCGCCTACGTTTTGTTTCCGGTTTGTCGCCCTGACAAGACCCACCGCGGCGCTCCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTc                                                                                              <  1:545494/148‑1 (MQ=9)
                                                       ttttttCCGGCTGCCGGATGCGGCGTGAACGCCATATCCAGCCTACGGTTTTGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTcc                                                                                             <  2:4898874/149‑1 (MQ=255)
                                                           ttACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGTCTGTCGGAGGCGGCGTAGGCGCCTTGTCAGACCTACGGTTATGttaagttg                                                                                         >  2:751935/1‑143 (MQ=17)
                                                           tcaCGGCGGCCGGTTGGGGCGTCAACGCCCTAACCCGACAACGGTTTTGTTCCGTTTTGTAGGCCTTATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGttt                                                                                         <  2:241015/147‑1 (MQ=14)
                                                            tccggctgccggctgcGGCTTTAACGCCCTATCATGCCTTGGTTTTTTTTCGTGTTTTTAGCCCTTAAAAGACGCAATGCGTCGCATCTGCCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTg                                                                                        <  1:1482449/147‑1 (MQ=9)
                                                            tACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTg                                                                                        >  1:4101557/1‑149 (MQ=255)
                                                            tACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTg                                                                                        >  1:4022700/1‑149 (MQ=255)
                                                            tACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTg                                                                                        >  1:3667211/1‑149 (MQ=255)
                                                              cGGCTGCCGGATGCGGCGTTAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTa                                                                                      <  2:3266866/149‑1 (MQ=255)
                                                               ggcggcCGGCGGCGGCGGGACCGCCCTAGCCAGCACCCGGTTTGGTTCGGGTTTGTTGGCCTTATACGACCCACCGCGTCGCATCCGCCCACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAg                                                                                     <  1:1602487/149‑1 (MQ=9)
                                                               ggCTGCCGGTTGCGGCTTGAACGCCCTAACCAGCCTACGTTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACCGCGTCGAATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAg                                                                                     <  1:1243730/149‑1 (MQ=17)
                                                               ggCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTAGGGTTTTGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAg                                                                                     <  2:3027681/149‑1 (MQ=255)
                                                               ggCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCACCGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGGTTGTAg                                                                                     >  2:3617895/1‑149 (MQ=255)
                                                                gCTGCCGGATGCGGCGTGAACGCCCTACCCCGCCTACGTTTATGTTCCGTTTTGTCGCCCTGATAAGACGCACCGCGTCGAAACAGCCCACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAgg                                                                                    <  1:2486501/149‑1 (MQ=14)
                                                                     cGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGACGGATGCGGCGTCAACGCCGTATCCGACCTACGGTTATGTTCGGGATCGAGGCCTGa                                                                               >  2:4795413/1‑149 (MQ=14)
                                                                        aTGGGGCTTGACCGCCCTACCCACCCTAGGGTATTGTTCCGTTTTTTCGGCCTTAAAAGACGCACAGCGTCGCATCAGCCCACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATaa                                                                            <  1:181357/149‑1 (MQ=14)
                                                                        aTGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTGTGTAGGCCTGATaa                                                                            >  1:290696/1‑149 (MQ=255)
                                                                          gcgggttGACCGCCCTATCCCGCCTACGGTTATGTGCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGa                                                                          <  2:1059150/143‑1 (MQ=14)
                                                                           cggcgTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGAc                                                                         <  2:1242349/149‑1 (MQ=255)
                                                                                tgtaCGCCCTAGCTAGCCTAGGTTTTTGTTCGGTTTTGTGGGCCTTATACGACGCACTGGGTCGCTTCAGGCCACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGcaca                                                                    <  1:3809580/146‑1 (MQ=14)
                                                                                 gAACGCCCTATCCAGCCTAGGGTTTTGTTCGGGTTTGTCGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAg                                                                   <  1:3381763/149‑1 (MQ=255)
                                                                                 gAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAg                                                                   >  1:4738060/1‑149 (MQ=255)
                                                                                     gCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGACGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTc                                                               >  2:3796909/1‑149 (MQ=255)
                                                                                         tATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCAt                                                           >  2:2007404/1‑149 (MQ=255)
                                                                                          aTCCGGCCTACGGTTATGTTCCGGTTTGTAGGCATGATAAGACGCACAGCGTCGCAAAAGCCAACGCCTGACGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATCAGACGCACAGCGTCGCATc                                                          >  1:2291995/1‑149 (MQ=11)
                                                                                          aTCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATc                                                          <  2:3437265/149‑1 (MQ=255)
                                                                                           tCCAGCCTCGGTTTTTGTTCGGGTTTGTAGGCCTCATAAGACGCACAGGGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCa                                                         <  1:791663/149‑1 (MQ=11)
                                                                                           tCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTTATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCa                                                         <  2:3517541/149‑1 (MQ=255)
                                                                                           tCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCa                                                         >  2:3379539/1‑149 (MQ=40)
                                                                                           tCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCa                                                         >  2:2706713/1‑149 (MQ=40)
                                                                                            ccAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGACGCATCAg                                                        >  2:3753947/1‑149 (MQ=255)
                                                                                                ccTAGGTTTATGGTCCGTTTTGTAGGCCTTATACGACGCACAGCGACGCATATGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCaa                                                    <  1:177855/149‑1 (MQ=11)
                                                                                                  tACGGTTTTGTTCCGGTTTGTAGGCCTGATAAGACGCACCGCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACg                                                  <  2:2512705/149‑1 (MQ=255)
                                                                                                       ttATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGc                                             >  1:3351180/1‑149 (MQ=37)
                                                                                                         aTGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCg                                           >  2:1065064/1‑149 (MQ=37)
                                                                                                           gTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGa                                         >  1:4251957/1‑149 (MQ=37)
                                                                                                            ttCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGAt                                        >  2:2892439/1‑149 (MQ=35)
                                                                                                            ttCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGACGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCAGTTTCTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGAt                                        >  1:2155769/1‑149 (MQ=16)
                                                                                                             tcgggTTTGTCGGCTTTTTACGCCGCCCGGCGTCGCATAGGGCACCGGGTGTGGGAGGCGGCTTAAACGCCTTATCCGACCTCCGGTTTTGTTCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATg                                       <  2:1570861/146‑1 (MQ=11)
                                                                                                               gtttttGTCGGCCTTATAATACCCCCCGGGTCGCAACTGGCAAGGGCTGTGGGTTGCGGCTTAACCGCCTTAACCGACATACGTTTTTTTTCCGTTTTTTGGCCTTATAAGACGCACCGCGTCGCATCAGGCAACGGCTGCCGGATgcg                                     <  2:1867836/146‑1 (MQ=2)
                                                                                                                   ttGTCGGCATTAAAAGACGCACCGCGTCGAATCTGCCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACATCCGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGt                                 <  2:431647/149‑1 (MQ=11)
                                                                                                                         gCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCAGGTTGTAGGCCTGATAAGACGCAAAGAGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGc                           >  1:1049716/1‑149 (MQ=255)
                                                                                                                          ctttaTACGACGAACTGGGTCGCATATGGAAACGGCTGTCGGTTGCGGGGTAAACGCCTTATCCGACCTACGGTTTTTTTCCTTTTGTTGGCCTTATAAGCCGCACCGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGcc                          <  2:5009464/145‑1 (MQ=11)
                                                                                                                            tGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGTTAAGACGCACAGCGTCGCATCAGGCACCGGCTGCCGGATGCGGCGTAAACGccat                        >  2:1036971/1‑147 (MQ=19)
                                                                                                                            tGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTCAACGccat                        >  1:2761136/1‑147 (MQ=17)
                                                                                                                              aTAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCatat                      >  1:4901899/1‑149 (MQ=21)
                                                                                                                              aTAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTAATACGACCAACGGTTATGATCCGGTTGTAGGCCTGATAATACGAACAGCGACGCATCAGGCAACGGCTGCCGGATGCGGCGGAAACGCCCTAt                      >  1:2264283/1‑149 (MQ=11)
                                                                                                                                aaGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACGGCGGCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTGATcc                    >  1:2323873/1‑149 (MQ=255)
                                                                                                                                aaGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATcc                    >  2:5152492/1‑149 (MQ=37)
                                                                                                                                 aGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCATATCCg                   >  1:1986860/1‑149 (MQ=21)
                                                                                                                                 aGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCATATCCg                   >  2:1961767/1‑149 (MQ=21)
                                                                                                                                  gACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGGTTGTTGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGAGGCGGCGAAAACCCCTTATCAGa                  >  1:1683862/1‑149 (MQ=11)
                                                                                                                                   aCGCACCGCGTCGCATCGGGCAACGGGTGTGGGTGGCGGGGTAAACGCCTTAGCCGCCCTCGGGTTTTGTTCGGTTTGTCGGCCTGATAAGACGCACCGCGTCGCATCTGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGAc                 <  2:3146435/149‑1 (MQ=11)
                                                                                                                                    cgccccgCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGTTTTTGTTCCGTTTGTAGGCCTGATAAGACGCACCGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGAcc                <  2:4689800/143‑1 (MQ=11)
                                                                                                                                    cGCACAGGGTCCCATCTGGCAACGGATGGCGGTTGCGGCGTAAACCCCTTATCCGACCTACGTTTTTTTTCGGTTTGTCGGCCTGATAAGACGCCCCGCGTCGCTTCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGAcc                <  2:602364/149‑1 (MQ=11)
                                                                                                                                    cGCACAGCGTCGCATCTGGAAACGGCTGGCGGATGGGGCTTAAACGCCTTATCCGACATACGTTTTTGTTCCGTTTGTCGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGAcc                <  2:2096845/149‑1 (MQ=11)
                                                                                                                                        cggCGTCGAATCTGGCAACGGCTGGCGGTTGCGGCGTAAACGCCTTAACCTACCCAGGTTTTTGTTCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATAAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACg            <  2:1978908/147‑1 (MQ=11)
                                                                                                                                          gCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGt          >  1:1295378/1‑149 (MQ=37)
                                                                                                                                             tCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTAt       >  1:2598644/1‑149 (MQ=37)
                                                                                                                                             tCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTAt       >  1:3030856/1‑149 (MQ=37)
                                                                                                                                                 aTCTGGCACCGGCTGTCGGTTGGGGCTTAACCGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCTGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTc   <  2:3456586/149‑1 (MQ=11)
                                                                                                                                                    |                                                                                                                                                  
CTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC  >  CP009273/3385368‑3385662

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: