Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,031,861 G→A 7.7% G18S (GGC→AGC)  hyaE → putative HyaA chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,031,8610GA 7.7% 676.2 / 26.7 194G18S (GGC→AGC) hyaEputative HyaA chaperone
Reads supporting (aligned to +/- strand):  ref base G (97/82);  new base A (10/5);  total (107/87)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 200 total aligned reads displayed.

TAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCC  >  CP009273/1031714‑1032006
                                                                                                                                                   |                                                                                                                                                 
 aaTGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCgg                                                                                                                                                  <  1:451356/148‑1 (MQ=255)
  aTGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGc                                                                                                                                                 >  2:1481944/1‑148 (MQ=255)
  aTGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGc                                                                                                                                                 <  1:1355060/148‑1 (MQ=255)
  aTGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCAGc                                                                                                                                                 >  1:3108059/1‑148 (MQ=255)
    gAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTg                                                                                                                                               >  1:1001390/1‑148 (MQ=255)
    gAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTg                                                                                                                                               >  2:2859278/1‑148 (MQ=255)
    gAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTg                                                                                                                                               >  2:3225298/1‑148 (MQ=255)
      aTCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGa                                                                                                                                             >  1:1681354/1‑148 (MQ=255)
       tCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGAc                                                                                                                                            >  1:2915287/1‑148 (MQ=255)
       tCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGAc                                                                                                                                            >  1:3204311/1‑148 (MQ=255)
         gcgcTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGc                                                                                                                                          <  2:3359688/148‑1 (MQ=255)
           gcTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCca                                                                                                                                        <  2:220667/148‑1 (MQ=255)
           gcTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGACCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCca                                                                                                                                        <  2:220619/148‑1 (MQ=255)
            cTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCcag                                                                                                                                       >  2:702191/1‑148 (MQ=255)
            cTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCcag                                                                                                                                       >  1:513007/1‑148 (MQ=255)
             tGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCcagt                                                                                                                                      <  1:1634798/148‑1 (MQ=255)
              gTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCcagtc                                                                                                                                     <  1:2009045/148‑1 (MQ=255)
              gTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCcagtc                                                                                                                                     >  2:81296/1‑148 (MQ=255)
               tctcAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCcagtca                                                                                                                                    >  2:1648957/1‑148 (MQ=255)
               tctcAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCcagtca                                                                                                                                    >  1:264509/1‑148 (MQ=255)
               tctcAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCcagtca                                                                                                                                    >  2:2065165/1‑148 (MQ=255)
                ctcAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCcagtcag                                                                                                                                   <  2:2964548/148‑1 (MQ=255)
                 tcAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCcagtcagt                                                                                                                                  <  1:112213/148‑1 (MQ=255)
                  cAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTg                                                                                                                                 >  2:1592012/1‑148 (MQ=255)
                    aTTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGaa                                                                                                                               <  2:272671/148‑1 (MQ=255)
                    aTTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGaa                                                                                                                               <  2:389338/148‑1 (MQ=255)
                     ttATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAAt                                                                                                                              >  2:2415361/1‑148 (MQ=255)
                      tATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATc                                                                                                                             >  1:1486961/1‑148 (MQ=255)
                         gACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCg                                                                                                                          >  1:2434702/1‑148 (MQ=255)
                          aCTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGt                                                                                                                         >  2:3422930/1‑148 (MQ=255)
                          aCTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGt                                                                                                                         >  1:108024/1‑148 (MQ=255)
                              tctgtcGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTg                                                                                                                     >  2:1482982/1‑148 (MQ=255)
                              tctgtcGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTg                                                                                                                     >  2:1823060/1‑148 (MQ=255)
                              tctgtcGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCAGTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTg                                                                                                                     >  2:3254117/1‑148 (MQ=255)
                                 gtcGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTgacg                                                                                                                  >  1:1658301/1‑148 (MQ=255)
                                  tcGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTgacga                                                                                                                 <  1:3480982/148‑1 (MQ=255)
                                  tcGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTgacga                                                                                                                 >  1:2411254/1‑148 (MQ=255)
                                  tcGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTgacga                                                                                                                 >  1:2499570/1‑148 (MQ=255)
                                   cGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTgacgac                                                                                                                <  1:2521935/148‑1 (MQ=255)
                                    gATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACt                                                                                                               <  1:2814046/148‑1 (MQ=255)
                                      tGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTgg                                                                                                             <  2:725510/148‑1 (MQ=255)
                                      aggAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGAGAATCCCGTCTTGACGACTgg                                                                                                             >  2:2290971/2‑148 (MQ=255)
                                            aTTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACg                                                                                                       <  2:1773018/148‑1 (MQ=255)
                                            aTTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACg                                                                                                       <  2:1772514/148‑1 (MQ=255)
                                            aTTACGAAGGCGGTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACg                                                                                                       <  2:1387035/148‑1 (MQ=255)
                                              tACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCa                                                                                                     >  2:672955/1‑148 (MQ=255)
                                              tACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCa                                                                                                     <  2:2123748/148‑1 (MQ=255)
                                              tACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCa                                                                                                     >  2:672910/1‑148 (MQ=255)
                                                cGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAg                                                                                                   >  2:2739065/1‑148 (MQ=255)
                                                  aaGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAgcg                                                                                                 >  1:2770417/1‑148 (MQ=255)
                                                    ggCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGcc                                                                                               >  1:1955600/1‑148 (MQ=255)
                                                    ggCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCAGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGcc                                                                                               >  1:1536972/1‑148 (MQ=255)
                                                        ttCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGAc                                                                                           >  1:369850/1‑148 (MQ=255)
                                                         tCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACg                                                                                          <  1:3139870/148‑1 (MQ=255)
                                                         tCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACg                                                                                          >  2:2855130/1‑148 (MQ=255)
                                                                cgcCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCgtggt                                                                                   >  2:567630/1‑148 (MQ=255)
                                                                  ccAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGgtgt                                                                                 <  2:92795/148‑1 (MQ=255)
                                                                     gTACCGGATGACACAGGAGGAGCAGGTATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGttat                                                                              <  1:2892890/148‑1 (MQ=255)
                                                                       aCCGGATGACACCGGCGGCGCTGGGATGAGCAACGACACGCAATTTGTTGCGTTGTGGCACCGAATGGTGGCGCGCAGCTGGCCGCCAGTCTGTGAAGCCCGTCTTGACGACTGGCTTCCGCAAGCGCCAGACGGCGTGGTGttatta                                                                            <  1:368443/148‑1 (MQ=255)
                                                                       aCCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGttatta                                                                            >  1:3353863/1‑148 (MQ=255)
                                                                             tGACACCGGAGGAGCTGGGATGAGCAACGACACGCCATTTGATGCGTTGGGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCACGCGCCAGACGGCGTGGTGTTATTAAGCAGt                                                                      <  1:3026474/148‑1 (MQ=255)
                                                                                cacaGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGAc                                                                   <  2:1252532/148‑1 (MQ=255)
                                                                                cacaGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGAc                                                                   >  2:1620947/1‑148 (MQ=255)
                                                                                cacaGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGAc                                                                   >  2:1620705/1‑148 (MQ=255)
                                                                                   aggaggAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCg                                                                >  1:836682/1‑148 (MQ=255)
                                                                                    ggaggaGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGa                                                               >  2:3262466/1‑148 (MQ=255)
                                                                                    ggaggaGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGa                                                               >  1:3279985/1‑148 (MQ=255)
                                                                                     gaggagCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCTTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGaa                                                              <  1:1046239/148‑1 (MQ=255)
                                                                                     gaggagCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGaa                                                              >  1:2073822/1‑148 (MQ=255)
                                                                                     gaggagCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCAGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGaa                                                              >  2:2891164/1‑148 (MQ=255)
                                                                                        gagCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCAGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAacg                                                           <  2:1536972/148‑1 (MQ=255)
                                                                                             gggATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAacgcacgc                                                      <  2:2610408/148‑1 (MQ=255)
                                                                                                 tGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCaga                                                  >  2:2636770/1‑148 (MQ=255)
                                                                                                   aGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAgg                                                <  1:888552/148‑1 (MQ=255)
                                                                                                    gCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGt                                               <  1:1475397/148‑1 (MQ=255)
                                                                                                    gCAACGACACGCCATTTGATGCGTTGTGGCAACGAAGGGTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGt                                               <  1:191439/148‑1 (MQ=255)
                                                                                                    gCAACGACACGCCATTTGATGCGTTGGGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGt                                               <  1:1991031/148‑1 (MQ=255)
                                                                                                     cAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTc                                              <  2:1751392/148‑1 (MQ=255)
                                                                                                     cAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTc                                              <  2:2634651/148‑1 (MQ=255)
                                                                                                       aCGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAg                                            >  2:875671/1‑148 (MQ=255)
                                                                                                        cGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGc                                           >  2:613324/1‑148 (MQ=255)
                                                                                                            acGCCATTTGATGCGTGGTGGCAACGAAGGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATa                                       <  1:40898/148‑1 (MQ=255)
                                                                                                             cGCCATTTGTTGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGTACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCACGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATaa                                      <  1:1253763/148‑1 (MQ=255)
                                                                                                              gCCATTTGTTGCGGTGTGGCAACGAATGCTGTCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAAt                                     <  2:3204311/148‑1 (MQ=255)
                                                                                                                cATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATcc                                   >  1:2617649/1‑148 (MQ=255)
                                                                                                                 aTTTGATGCGTTGTGGCAACGAATGCTGGCGCGCAGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCg                                  <  2:3108059/148‑1 (MQ=255)
                                                                                                                   ttgttgCGTTGTGGCAACGATTGCTGGCGCGCGGCTGGACGCCAGTCTGTTAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGt                                <  1:1823060/148‑1 (MQ=255)
                                                                                                                        gCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGa                           <  2:3557232/148‑1 (MQ=255)
                                                                                                                            tgtgGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTgg                       >  2:717437/1‑148 (MQ=255)
                                                                                                                                gCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGaa                   <  1:2197713/148‑1 (MQ=255)
                                                                                                                                  aaCGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAAtt                 >  1:2895900/1‑148 (MQ=255)
                                                                                                                                  aaCGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAAtt                 >  2:2487464/1‑148 (MQ=255)
                                                                                                                                   aCGATTGCTGGCGCGCGGCTGTACGCATGTCAGTGAATCCCTTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTa                <  2:2434702/148‑1 (MQ=255)
                                                                                                                                   aCGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTa                <  1:2583100/148‑1 (MQ=255)
                                                                                                                                     gAATGCTGGCGCGCAGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACt              >  1:1427208/1‑148 (MQ=255)
                                                                                                                                       aTGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTgc            >  1:1226283/1‑148 (MQ=255)
                                                                                                                                        tGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTgcg           <  1:2327658/148‑1 (MQ=255)
                                                                                                                                            ggCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAg       <  2:1023601/148‑1 (MQ=255)
                                                                                                                                                cgcgGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCACGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTc   <  2:369850/148‑1 (MQ=255)
                                                                                                                                                cgcgGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTc   <  1:2737670/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TAATGAATCGCGCTGTCTCAATTATGACTGTCTGTCGATGGAAAATTACGAAGGCGTTCGCTTGCGCCAGTACCGGATGACACAGGAGGAGCAGGGATGAGCAACGACACGCCATTTGATGCGTTGTGGCAACGAATGCTGGCGCGCGGCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCC  >  CP009273/1031714‑1032006

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