Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,080,755 T→G 21.0% F231V (TTC→GTC)  metK → S‑adenosylmethionine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,080,7550TG21.0% 503.7 / 70.9 296F231V (TTC→GTC) metKS‑adenosylmethionine synthetase
Reads supporting (aligned to +/- strand):  ref base T (77/152);  new base G (64/3);  total (141/155)
Fisher's exact test for biased strand distribution p-value = 2.07e-21
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 314 total aligned reads displayed.

CGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACG  >  CP009273/3080607‑3080903
                                                                                                                                                    |                                                                                                                                                     
  aTGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttc                                                                                                                                                     <  2:3114881/149‑1 (MQ=255)
  aTGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttc                                                                                                                                                     <  2:514272/149‑1 (MQ=255)
    gCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttcgt                                                                                                                                                   <  2:2039042/149‑1 (MQ=255)
     cTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttcgtt                                                                                                                                                  <  1:904840/149‑1 (MQ=255)
      tGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTa                                                                                                                                                 <  2:2439333/149‑1 (MQ=255)
          gagCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCgg                                                                                                                                             >  2:3827464/3‑149 (MQ=255)
          gagCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCgg                                                                                                                                             <  1:3907742/147‑1 (MQ=255)
             cTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTgg                                                                                                                                          >  1:1583989/1‑149 (MQ=255)
               ttCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGGTATCGGTgggc                                                                                                                                        >  1:2004423/1‑147 (MQ=255)
                  cACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCaa                                                                                                                                    >  2:3370500/1‑149 (MQ=255)
                     tCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGGTATAGGTGGCCCAAAgg                                                                                                                                 >  1:2328287/1‑149 (MQ=255)
                          aCTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGAc                                                                                                                            <  2:1958166/149‑1 (MQ=255)
                           ctctGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCAATGGGTGACt                                                                                                                           >  2:2276430/1‑149 (MQ=255)
                              tGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTGGTTATCGGTGGCCCCATGGGTGACAGCg                                                                                                                        >  1:1415069/1‑149 (MQ=255)
                              tGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCg                                                                                                                        <  1:2825313/149‑1 (MQ=255)
                              tGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCg                                                                                                                        <  1:2598103/149‑1 (MQ=255)
                              tGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATAGGTGGCCCCAAGGGTGACAGCg                                                                                                                        >  1:4625049/1‑149 (MQ=255)
                              tGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTGGGTATCTGTGGCCCAATGGGTGACTGCg                                                                                                                        >  1:2182010/1‑149 (MQ=255)
                               gAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCgg                                                                                                                       <  2:2892322/149‑1 (MQ=255)
                                aaGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTGGGTATCGGGGGGCCCAAGGGTGGcggcggg                                                                                                                      >  1:3063078/1‑148 (MQ=255)
                                 agagATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGATTGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGGTTCGTTATCGGTGGCCCAATGGGTGACTGCggac                                                                                                                     <  2:2945055/149‑3 (MQ=255)
                                    gaTCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGa                                                                                                                  >  2:2745482/1‑149 (MQ=255)
                                       cGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTg                                                                                                               <  1:3627137/149‑1 (MQ=255)
                                       cGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTg                                                                                                               <  1:3032656/149‑1 (MQ=255)
                                        gACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTgg                                                                                                              >  2:3336482/1‑149 (MQ=255)
                                         aCCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGt                                                                                                             <  1:4253805/149‑1 (MQ=255)
                                         aCCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGt                                                                                                             >  2:589233/1‑149 (MQ=255)
                                          ccAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTc                                                                                                            <  1:2078916/149‑1 (MQ=255)
                                            aGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGt                                                                                                          <  2:922681/149‑1 (MQ=255)
                                               aaTCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCCAAGGGTGACTGCGGTCAGACTGGTCGGaaa                                                                                                       >  1:1675677/1‑149 (MQ=255)
                                                ttCGCTGCAAGTAGAGGTAATGGAAGAGATAATCAAGCCAATTCTGACCTCTGTATTGCTGACTTCTGCATCCTAATTCTTCATCAACCCGACCGGTTGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGTCTGGTCGTataa                                                                                                      <  2:1575595/148‑1 (MQ=255)
                                                aTCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTaaaa                                                                                                      >  1:4179442/1‑149 (MQ=255)
                                                aTCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGGTATCGGTGGGCCAAAGGGTGACTGCGGTCTGTCAGGTCGTaaaa                                                                                                      >  1:782968/1‑149 (MQ=255)
                                                aTCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGGGGCCCCCTGGGGGACTGCGGTCTGACTGGTCGTgaaa                                                                                                      >  1:3600548/1‑149 (MQ=255)
                                                 tCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAAt                                                                                                     <  2:711862/149‑1 (MQ=255)
                                                 tCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCCAAGGGTGACTGCGGTCTGACTGGTCGTAAAAt                                                                                                     >  1:3854355/1‑149 (MQ=255)
                                                  cGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAAtt                                                                                                    <  2:4001644/149‑1 (MQ=255)
                                                   gCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTa                                                                                                   <  1:4370215/149‑1 (MQ=255)
                                                    cTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTAt                                                                                                  <  1:2633171/149‑1 (MQ=255)
                                                    cTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGGTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTAt                                                                                                  <  2:4974998/149‑1 (MQ=255)
                                                     tGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTGGTTATAGGGGGGCCAATGGGTGACTGCGGTCTGACTGGTCGTTAAAatata                                                                                                 >  1:5051846/1‑148 (MQ=255)
                                                     tGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATc                                                                                                 <  1:1860484/149‑1 (MQ=255)
                                                        aagaagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGGTATGGGTGGCCCAAAGGGTGACAGCGGTCTTACTGGTGGTAAAAATATAGtt                                                                                              >  1:95480/1‑149 (MQ=255)
                                                         agaagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTg                                                                                             <  1:4956203/149‑1 (MQ=255)
                                                           aagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGAt                                                                                           <  2:4600565/149‑1 (MQ=255)
                                                           aagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGAt                                                                                           <  1:1024828/149‑1 (MQ=255)
                                                              cGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATAcc                                                                                        <  1:805441/149‑1 (MQ=255)
                                                               ggTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCt                                                                                       <  2:1721889/149‑1 (MQ=255)
                                                                 tAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTAc                                                                                     <  2:1558462/149‑1 (MQ=255)
                                                                     ggAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTTAAATTATCGTTGATACCTAcggcg                                                                                 >  2:2804907/1‑149 (MQ=255)
                                                                      gAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATAGGTGGCCCAATGGGTTACTGCGGTCTGACTGGTCGGAAAATTATCGTTGATACCAAcggcgg                                                                                >  2:479836/1‑149 (MQ=255)
                                                                       aaGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACggcggc                                                                               >  2:3324850/1‑149 (MQ=255)
                                                                       aaGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTAGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACggcggc                                                                               <  2:796527/149‑1 (MQ=255)
                                                                          agaTCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATg                                                                            <  2:5071565/149‑1 (MQ=255)
                                                                          agaTCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCCGTGGCCCCAAGGGTGTCTGCGGTCTGACGGGTCGTAAAATTATCGTTGATAACTACGGCGGCCTg                                                                            >  2:4449744/1‑149 (MQ=255)
                                                                          agaTCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGCTCGTGTAGTTATCCGAGGCCCCAAGGGTGACTGCGGTCAGACAGGTCGTAAAAATATCTATGATAACAACGGCGGCAtt                                                                            >  2:4656666/1‑148 (MQ=255)
                                                                           gaTCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATgg                                                                           <  2:3856766/149‑1 (MQ=255)
                                                                           gaTCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATgg                                                                           <  2:566239/149‑1 (MQ=255)
                                                                            atcatcAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGgc                                                                          <  1:4788689/149‑1 (MQ=255)
                                                                            atcatcAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGgc                                                                          >  2:4532693/1‑149 (MQ=255)
                                                                              catcaAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGGTATCGGTGGCCCCAAGGGTGACTGCGGTCTGACTGGTCGTAAAATTATAGTTGATAACAACCGCGGCATGgcgc                                                                        >  1:2922820/1‑149 (MQ=255)
                                                                               atcaAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATGGGTGGCCCAAATGGTTACTGCGGTCTGACGGGTCGTTAAATTATCGGTGATAACTAAGGCGGCATGGcgcg                                                                       >  2:21549/1‑149 (MQ=255)
                                                                                tcaAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCCAAGGGTGTCTGCGGTCTGACTGGTCGTAAAATTATAGTTGATAACTACGGCGGCATGGCGCGt                                                                      >  2:1426530/1‑149 (MQ=255)
                                                                                  aaGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCa                                                                    >  2:738985/1‑149 (MQ=255)
                                                                                   aGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCAc                                                                   <  1:1520164/149‑1 (MQ=255)
                                                                                    gCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGGTATCGGTGGGCCCAAGGGGGACTGCGGGATTTCTGGTGGGTAAAATATAGTTGATAAAAAAGGCGGCATGGGGCGTCCCg                                                                  >  1:1372767/1‑149 (MQ=255)
                                                                                     ccAATTCTGCCCGCTGAATGGCTGACTTCTGCCCCCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTTTCCGTGGCCCCAATGGGTAATGCGGTCTTACTTGTAGTAAAATTATCGTTTTTAACTAACGCGGCCTGGCGCGTCAccg                                                                 >  2:1861651/1‑147 (MQ=255)
                                                                                         ttCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGc                                                             <  1:4976369/149‑1 (MQ=255)
                                                                                         ttCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTAGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGc                                                             >  1:1965113/1‑149 (MQ=255)
                                                                                          tCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCg                                                            <  2:1473344/149‑1 (MQ=255)
                                                                                                cGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCAt                                                      <  1:4522279/149‑1 (MQ=255)
                                                                                                     aaTGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGGGGCCCAAGGGGTGACTGCGGTCTGACGGGTCGTAAAATTATCGTTGATACCTACGGCGGCAAGGCGCGTCACGGTGGCGGGGCATtctca                                                 >  2:2689621/1‑148 (MQ=255)
                                                                                                     aaTGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGGTATAGGTGGCCCCAAGGGTGACTGCGGTCGGACTGGTCGTAAAAATATAGTTGTTAACAACCGCGGCATGGCGCGTCCCCGTGGCGGGGCATTctct                                                 >  1:4235124/1‑149 (MQ=255)
                                                                                                         gCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTAAGGCGGCATTGCGCGTCACGGTGGCGGTGCATTCTCTGGTa                                             >  2:2514073/1‑149 (MQ=255)
                                                                                                            gACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGTCTGCGGTCTGACTGGTCGTAAAATTATAGTTTATTCCTACGGCGGCCTGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAg                                          >  2:4382266/1‑149 (MQ=255)
                                                                                                            gACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGGTAACGGTGGCCCAATGGGTGAAGGCGGGCTGACTGGTCGTAAAATTATCGTTGATACCAACGGCGGCAGTGCGCGGCACCGTGGCGGGGCAATCACAGGTaaat                                          >  1:1273359/1‑148 (MQ=255)
                                                                                                             aCTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCCAAGGGTGACAGCGGTCTGACAGGTCGGAAAATTATAGTTGATACCAACCGCGGCAATGCGCGGCACCGGGGCCGTGCAATCTCAGGGaaata                                         >  1:3352166/1‑147 (MQ=255)
                                                                                                                 cTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCa                                     <  1:3494645/149‑1 (MQ=255)
                                                                                                                  tGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCAt                                    <  1:2535192/149‑1 (MQ=255)
                                                                                                                   gCCACCAAATTCTTCCTCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATc                                   <  2:1183978/149‑1 (MQ=255)
                                                                                                                   gCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACGGCGGTCTGACTGGTCGTAAAAATATCGGTTATAACAAAGGCGGCATGGCGCGGCACGGTGGCGGTGCATTATCTGGTAAATATCCAta                                   >  2:2415420/1‑148 (MQ=255)
                                                                                                                      accaAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACGGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa                                >  2:4279550/1‑149 (MQ=255)
                                                                                                                      accaAATTCTTCATCAACCCGACCGGTCGTTTCGTTATAGGTGGCCCAAATGGTGACTGCGGTCTGACTGGTCGTAAAAATATAGTTGATACCAAAGGCGGCAAGGCGCGTCACCGTGGCGGTGCATTCTCAGGTAAAGATAAATCaaa                                >  1:1192066/1‑149 (MQ=255)
                                                                                                                      accaAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATAGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAATATACATCaaa                                >  2:1645354/1‑149 (MQ=255)
                                                                                                                          aaTTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTg                            <  1:5150335/149‑1 (MQ=255)
                                                                                                                          aaTTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTg                            <  1:1384452/149‑1 (MQ=255)
                                                                                                                          aaTTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCGGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCCCGGTGGCGGGGCATTCTCTGGTAAAGATCCAAAAAAAGTg                            >  2:2070394/1‑149 (MQ=255)
                                                                                                                            ttcttcATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGa                          <  1:4263626/149‑1 (MQ=255)
                                                                                                                             tcttcATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGAc                         <  1:1801944/149‑1 (MQ=255)
                                                                                                                             tcttcATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGAc                         <  2:3002764/149‑1 (MQ=255)
                                                                                                                             tcttcATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGAc                         <  1:2510115/149‑1 (MQ=255)
                                                                                                                                 catcaACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGtt                     <  1:856283/149‑1 (MQ=255)
                                                                                                                                 catcaACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCGGACGGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGGGGACCggt                     >  2:3263720/1‑147 (MQ=255)
                                                                                                                                  atcaACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                    <  2:4179442/149‑1 (MQ=255)
                                                                                                                                  atcaACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                    <  2:2701165/149‑1 (MQ=255)
                                                                                                                                     aaCCCGACCGGTCGTTTAGTTATAGGTGGCCCAAAGGGTGACTGCGGTCTGACTGGTCGTAAAATTATAGTTGATAACAACGGCGGCATGGCGCGTCACGGTGGCGGTGCCTTCTCTGGTAAACAAACATCAAAATTGGACCGTTcccc                 >  2:3603879/1‑147 (MQ=255)
                                                                                                                                      aCCCGACCGGTCGTTTAGGTATCGGGGGCCCACAGGGTGACAGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTAAGGCGGCATGGCGCGGCACGGTGGCGGTGCAATCTCTGGTAAAGATACATCAAAAGTGGACCGTTCCGCa                >  2:3094819/1‑149 (MQ=255)
                                                                                                                                       cccGACCGGTCGTGTCGATATCGGTGGCCCAAATGGTGAATGCGGTCAGAATGGTCGTAAAAATATAGTTTATACCTACGGCGGCATGGCGCGTCACGGTGGAGGTGCATTATCTTTTAAACATCCAAAAAAAAAGGACCGTTCCCCac               >  2:106306/1‑148 (MQ=255)
                                                                                                                                           aCCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTa           <  1:1257635/149‑1 (MQ=255)
                                                                                                                                                   tttcgttATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCAc   >  2:875988/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACG  >  CP009273/3080607‑3080903

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