Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,711,755 C→A 100% D45Y (GAC→TAC)  yfiF ← putative methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,711,7550CA98.7% 720.5 / 0.5 225D45Y (GAC→TAC) yfiFputative methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (3/0);  new base A (107/115);  total (110/115)
Fisher's exact test for biased strand distribution p-value = 1.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 226 total aligned reads displayed.

TTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAG  >  CP009273/2711612‑2711902
                                                                                                                                               |                                                                                                                                                   
tttACCCCTGATGCCCCGGTTATGGGCATTTTCCGGCGTCTCACCACCCGGCGCGCGGTAAACCGGGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCa                                                                                                                                                 <  2:1299589/148‑1 (MQ=255)
    ccACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACgag                                                                                                                                             >  1:78715/1‑148 (MQ=255)
      agtgcggCCACCGTGCGCGCCCTTTTCCGGCGTCCCAACCCCCGGCGCGCGGCAAACCGGGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGAt                                                                                                                                           <  1:2311741/142‑1 (MQ=255)
        tgttcccACGGTCAGCGGCCTTGTCGGGCGTCTCACCCCCCGGCGCGGGGCAACCCGGCCCCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGAttt                                                                                                                                         <  1:1624857/143‑1 (MQ=255)
           ggccccggTGAGCGGCCTTTGCCGGCGTCTCACCCCCCGGCGCGGGGTAACCGGGGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTcc                                                                                                                                      <  1:3523393/141‑1 (MQ=255)
             ccACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTg                                                                                                                                    >  2:2212734/1‑148 (MQ=255)
                ccGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGAc                                                                                                                                 >  1:1026083/1‑148 (MQ=255)
                  gTGATCGGCATTTGCCGGCTTCTCATCACCCGGCGCGCGGCAACCCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGAcgc                                                                                                                               <  2:1905740/148‑1 (MQ=255)
                  gTGAGCGGCCTTTGCCGGCGTCTCATCCCCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGAcgc                                                                                                                               <  2:102213/148‑1 (MQ=255)
                   tGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGAcgcc                                                                                                                              >  1:1383162/1‑148 (MQ=255)
                     aTCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGAcgcccg                                                                                                                            >  1:2514615/1‑148 (MQ=255)
                       cGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTGTCCTGGAcgcccgcc                                                                                                                          >  2:2079763/1‑148 (MQ=255)
                        ggCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACgcccgccc                                                                                                                         >  1:921168/1‑148 (MQ=255)
                        ggCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACgcccgccc                                                                                                                         >  1:1841214/1‑148 (MQ=255)
                        ggCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACgcccgccc                                                                                                                         >  1:1936981/1‑148 (MQ=255)
                        ggCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACgcccgccc                                                                                                                         >  1:2228572/1‑148 (MQ=255)
                         gCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGcccc                                                                                                                        >  2:1492368/1‑148 (MQ=255)
                          ccTTTTCCGGCGTCTCATCCCCCGGCGCGCGGCAACCCGGGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGccccc                                                                                                                       <  2:478764/148‑1 (MQ=255)
                            ttttCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCtt                                                                                                                     >  2:595359/1‑148 (MQ=255)
                             tttCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACTGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCttt                                                                                                                    >  2:2445426/1‑148 (MQ=255)
                                       ctcATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTAcgc                                                                                                          >  2:507517/1‑148 (MQ=255)
                                         catcaCCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCgg                                                                                                        >  2:3338968/1‑148 (MQ=255)
                                          atccCCGGGCGCGCGGTAACCGGGGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCggg                                                                                                       <  1:2191290/148‑1 (MQ=255)
                                          atcaCCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTCCGCggg                                                                                                       >  1:116067/1‑148 (MQ=255)
                                           tcaCCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGt                                                                                                      >  1:1895451/1‑148 (MQ=255)
                                             ccccGGCGCGCGGGAACCCGGGCGCCACGGCGAATCTTCCCACGTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTtg                                                                                                    <  2:2719792/147‑1 (MQ=255)
                                               ccGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTtgtg                                                                                                  >  2:2815988/1‑148 (MQ=255)
                                                gggCGCGCGGGACCCCGTGCCCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTgg                                                                                                 <  1:285333/147‑1 (MQ=255)
                                                 ggCGCGCGGGAACCGGGCCCCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTggg                                                                                                <  1:1822371/148‑1 (MQ=255)
                                                     cgcgGTAACCCGGGCGCCACGGCGAATCTTCCCACGTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACg                                                                                            <  2:2790159/148‑1 (MQ=255)
                                                          gAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTAt                                                                                       >  2:1559567/1‑148 (MQ=255)
                                                             aCCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAg                                                                                    >  1:1509325/1‑148 (MQ=255)
                                                              cggggccccACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGa                                                                                   <  2:330811/141‑1 (MQ=255)
                                                               cgggcgcCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGaa                                                                                  <  2:3098686/145‑1 (MQ=255)
                                                                gTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAat                                                                                 <  2:1707980/148‑1 (MQ=255)
                                                                gTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAat                                                                                 >  2:550984/1‑148 (MQ=255)
                                                                  gcgcCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAatca                                                                               >  2:870092/1‑148 (MQ=255)
                                                                   cgcCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAatcat                                                                              <  2:1355886/148‑1 (MQ=255)
                                                                               tCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGa                                                                  >  2:3285994/1‑148 (MQ=255)
                                                                               tCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGa                                                                  >  1:480445/1‑148 (MQ=255)
                                                                                cTTCCCACTCGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGac                                                                 <  2:3247413/148‑1 (MQ=255)
                                                                                 ttCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGaca                                                                >  1:1137331/1‑148 (MQ=255)
                                                                                 ttCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGaca                                                                >  1:50632/1‑148 (MQ=255)
                                                                                    ccACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGacatac                                                             <  2:142906/148‑1 (MQ=255)
                                                                                        ttGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCa                                                         >  2:2484827/1‑148 (MQ=255)
                                                                                        ttGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCa                                                         >  1:2978383/1‑148 (MQ=255)
                                                                                         tGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCAc                                                        >  2:715690/1‑148 (MQ=255)
                                                                                          gCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACt                                                       <  1:1513159/148‑1 (MQ=255)
                                                                                           cGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACtt                                                      <  2:2276832/148‑1 (MQ=255)
                                                                                            ggTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACttt                                                     <  2:1620718/148‑1 (MQ=255)
                                                                                               cACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGAc                                                  <  1:3115915/148‑1 (MQ=255)
                                                                                               cACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGAc                                                  <  1:848065/148‑1 (MQ=255)
                                                                                                aCGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCTTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGAcc                                                 <  2:3115681/148‑1 (MQ=255)
                                                                                                aCGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGAcc                                                 <  1:2875653/148‑1 (MQ=255)
                                                                                                aCGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGAcc                                                 <  2:1560973/148‑1 (MQ=255)
                                                                                                 cGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCt                                                >  1:1730891/1‑148 (MQ=255)
                                                                                                  ggggCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCtt                                               >  1:3036457/1‑148 (MQ=255)
                                                                                                  ggggCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCtt                                               >  1:3036080/1‑148 (MQ=255)
                                                                                                      cTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCg                                           <  1:2453904/148‑1 (MQ=255)
                                                                                                      cTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCg                                           <  1:3338968/148‑1 (MQ=255)
                                                                                                       tGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGc                                          <  2:1534077/148‑1 (MQ=255)
                                                                                                         aCTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCtt                                        <  2:2966245/148‑1 (MQ=255)
                                                                                                           tctGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCtttt                                      <  1:3231671/148‑1 (MQ=255)
                                                                                                            ctGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTa                                     >  2:1523789/1‑148 (MQ=255)
                                                                                                              gTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTAcc                                   >  2:426530/1‑148 (MQ=255)
                                                                                                               tttGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCt                                  <  2:1457295/148‑1 (MQ=255)
                                                                                                                ttGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCtt                                 <  1:656328/148‑1 (MQ=255)
                                                                                                                ttGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCtt                                 >  2:254665/1‑148 (MQ=255)
                                                                                                                 tGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCttt                                <  1:1377810/148‑1 (MQ=255)
                                                                                                                  gTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCtttc                               <  2:3507390/148‑1 (MQ=255)
                                                                                                                     aTCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCtttcatt                            <  1:1158168/148‑1 (MQ=255)
                                                                                                                     aTCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCtttcatt                            >  2:423215/1‑148 (MQ=255)
                                                                                                                     aTCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCtttcatt                            >  1:1089795/1‑148 (MQ=255)
                                                                                                                     aTCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCtttcatt                            <  1:3549154/148‑1 (MQ=255)
                                                                                                                     aTCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCtttcatt                            >  2:1939106/1‑148 (MQ=255)
                                                                                                                      tCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTGTTTCCtttcattt                           >  1:72163/1‑148 (MQ=255)
                                                                                                                      tCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCACGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGAACTTGCCGCTTTTACCtttcattt                           >  2:1893185/1‑148 (MQ=255)
                                                                                                                       cgcgGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCtttcatttc                          >  1:3174166/1‑148 (MQ=255)
                                                                                                                         cgggcgggGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTttcatttcat                        <  1:1248257/148‑1 (MQ=255)
                                                                                                                           ggcggggcgACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCg                      <  2:78715/148‑1 (MQ=255)
                                                                                                                           ggcggggcgACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCg                      <  2:1077307/148‑1 (MQ=255)
                                                                                                                            gcggggcgACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGt                     >  2:667038/1‑148 (MQ=255)
                                                                                                                             cggggcgACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGtt                    <  1:1932746/148‑1 (MQ=255)
                                                                                                                             cggggcgACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGtt                    <  1:870118/148‑1 (MQ=255)
                                                                                                                                 gcgACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATg                <  1:2687237/148‑1 (MQ=255)
                                                                                                                                 gcgACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATg                <  1:2661696/148‑1 (MQ=255)
                                                                                                                                 gcgACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATg                >  2:1978265/1‑148 (MQ=255)
                                                                                                                                 gcgACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATGTCATCGTTCAtt                >  2:2271684/1‑147 (MQ=255)
                                                                                                                                   gACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCt              <  1:504650/148‑1 (MQ=255)
                                                                                                                                    aCGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCtt             >  1:2264523/1‑148 (MQ=255)
                                                                                                                                    aCGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCtt             <  1:1426760/148‑1 (MQ=255)
                                                                                                                                      ggCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCtttt           <  2:2520428/148‑1 (MQ=255)
                                                                                                                                       gCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTtc          >  2:3057947/1‑148 (MQ=255)
                                                                                                                                        ccACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTtct         >  1:1492794/1‑148 (MQ=255)
                                                                                                                                        ccACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTtct         >  1:2063403/1‑148 (MQ=255)
                                                                                                                                         cACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTtctc        <  1:2528867/148‑1 (MQ=255)
                                                                                                                                         cACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTtctc        <  1:550984/148‑1 (MQ=255)
                                                                                                                                         cACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTtctc        <  2:26728/148‑1 (MQ=255)
                                                                                                                                           ccGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTcca      <  2:2430728/148‑1 (MQ=255)
                                                                                                                                               cAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCGTTTACCTTTCATTTAATCGTTCATGCGTTTATCCACCAg  >  1:2817968/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAG  >  CP009273/2711612‑2711902

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