Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,080,755 T→G 20.3% 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,7550TG20.3% 405.2 / 49.7 233F231V (TTC→GTC) metKS‑adenosylmethionine synthetase
Reads supporting (aligned to +/- strand):  ref base T (60/122);  new base G (49/2);  total (109/124)
Fisher's exact test for biased strand distribution p-value = 4.65e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 244 total aligned reads displayed.

CGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCAC  >  CP009273/3080607‑3080902
                                                                                                                                                    |                                                                                                                                                   
     cTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtgtggtt                                                                                                                                                >  1:4176127/1‑149 (MQ=255)
      tGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTa                                                                                                                                               <  1:989782/149‑1 (MQ=255)
       gTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTAt                                                                                                                                              <  2:2197171/149‑1 (MQ=255)
       gTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTAt                                                                                                                                              >  1:2799270/1‑149 (MQ=255)
       gTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGGTAt                                                                                                                                              >  1:4546394/1‑149 (MQ=255)
        tCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTGGGTATc                                                                                                                                             >  1:28883/1‑149 (MQ=255)
          gagCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCgg                                                                                                                                           <  1:4436043/147‑1 (MQ=255)
           agCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGGTATAggg                                                                                                                                          >  1:506347/2‑148 (MQ=255)
                tCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGccc                                                                                                                                     <  1:912656/149‑1 (MQ=255)
                  cACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATAGGTGGCCCaa                                                                                                                                   >  1:3019966/1‑149 (MQ=255)
                   aCTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAAt                                                                                                                                  <  2:670355/149‑1 (MQ=255)
                    cTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATg                                                                                                                                 <  1:1559931/149‑1 (MQ=255)
                    cTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAAGGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACGTCTGCCACCATATTCGTCATCAACCCGACCGGTCGGGTCGGTATGGGGGGcgcgctg                                                                                                                                 >  2:2909136/1‑143 (MQ=255)
                      cAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCCATggg                                                                                                                               >  1:1668027/1‑149 (MQ=255)
                       aGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGGTATCGGTGGGCCAATGGGt                                                                                                                              >  1:1502748/1‑149 (MQ=255)
                        gCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTg                                                                                                                             <  2:3064848/149‑1 (MQ=255)
                        gCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGGTATCGGTGGCCCAATGGGTg                                                                                                                             >  2:1534154/1‑149 (MQ=255)
                        gCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCAAAGGGTg                                                                                                                             >  2:648166/1‑149 (MQ=255)
                              tGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGCCTGCg                                                                                                                       >  2:3747702/1‑149 (MQ=255)
                                aaGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGt                                                                                                                     <  1:3977220/149‑1 (MQ=255)
                                aaGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGt                                                                                                                     <  2:819368/149‑1 (MQ=255)
                                   agaTCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTg                                                                                                                  <  2:4062732/149‑1 (MQ=255)
                                     aTCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGAc                                                                                                                >  2:2004220/1‑149 (MQ=255)
                                     aTCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATAGGTGGCCCAATGGGTGACTGCGGTCTGAc                                                                                                                >  1:2494896/1‑149 (MQ=255)
                                     aTCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTGGGTATCGGGGGCCCCAAGGGTGTCTGCGGTCTGAc                                                                                                                >  1:1005240/1‑149 (MQ=255)
                                       cGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTg                                                                                                              <  1:305674/149‑1 (MQ=255)
                                          ccAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGTCTGCGGTCTGACTgggc                                                                                                           >  2:3295569/1‑147 (MQ=255)
                                            aGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTGTGGCTGGTCGt                                                                                                         >  2:4591116/1‑149 (MQ=255)
                                             gAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTGGTTATAGGTGGCCCAATGGGTGTCTGCGGTCTGACTGGTCGtt                                                                                                        >  1:130136/1‑148 (MQ=255)
                                             gAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTa                                                                                                        >  2:4476268/1‑149 (MQ=255)
                                              aaaTCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGTCTGCGGTCTGACTGGTCGTaa                                                                                                       >  1:2942627/1‑149 (MQ=255)
                                                aTCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCCAAGGGTGACTGCGGTCTGACTGGTCGTaaaa                                                                                                     >  2:1013811/1‑149 (MQ=255)
                                                aTCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTaaaa                                                                                                     <  1:1594239/149‑1 (MQ=255)
                                                 tCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAAt                                                                                                    <  1:3341981/149‑1 (MQ=255)
                                                  cGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAAtt                                                                                                   <  1:3846283/149‑1 (MQ=255)
                                                     tGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATc                                                                                                >  1:2555706/1‑149 (MQ=255)
                                                        aagaagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCAAATGGTGACTGCGGTCTGACAGGTCGTAAAAATATAGtt                                                                                             >  1:1228274/1‑149 (MQ=255)
                                                         agaagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTg                                                                                            <  1:713932/149‑1 (MQ=255)
                                                         agaagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGGTATCGGGGGCCCCAAGGGTGTCTGCGGTCTGACGGGGCGTAAAATTATCGTTg                                                                                            >  2:2833/1‑149 (MQ=255)
                                                          gaagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGa                                                                                           >  2:3033011/1‑149 (MQ=255)
                                                           aagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGGCGTTATCGGTGGCCCAAAGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGGTGAt                                                                                          >  1:285819/1‑149 (MQ=255)
                                                            agCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATa                                                                                         <  2:3831467/149‑1 (MQ=255)
                                                                 tAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTAc                                                                                    <  2:4932901/149‑1 (MQ=255)
                                                                 tAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTAc                                                                                    <  1:913075/149‑1 (MQ=255)
                                                                 tAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTAc                                                                                    <  1:543677/149‑1 (MQ=255)
                                                                 tAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTAc                                                                                    <  1:543620/149‑1 (MQ=255)
                                                                 tAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTAc                                                                                    >  2:4502600/1‑149 (MQ=255)
                                                                  aaTGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTAcg                                                                                   <  2:1822809/149‑1 (MQ=255)
                                                                   aTGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAAAGGGTGTCTGCGGTCTGTCTGGTCGTTAAATTATCGTTGTTACCAAccg                                                                                  >  1:2698603/1‑147 (MQ=255)
                                                                   aTGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGGTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTAcgg                                                                                  <  2:1456785/149‑1 (MQ=255)
                                                                      gAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTAcggcgg                                                                               <  1:4163011/149‑1 (MQ=255)
                                                                       aaGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGGTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACggcggc                                                                              <  1:75967/149‑1 (MQ=255)
                                                                        agagATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCa                                                                             <  2:303525/149‑1 (MQ=255)
                                                                            atcatcAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGgc                                                                         <  2:2642010/149‑1 (MQ=255)
                                                                               atcaAGCCGATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGGTATCGGTGGCCCCATGGGTGACTGCGGTCTGACTGGTCGGAAAATTATAGTTGATACCAACGGCGGCATGGcgcg                                                                      >  2:5120334/1‑149 (MQ=255)
                                                                               atcaAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTGGGTATAGGGGGCCCCAAGGGGGGCGGCGGGCTGACGGGGGGGAAAATTATCGTGGTTACCGACGGCGGCCGGGcgcg                                                                      >  1:2268048/1‑149 (MQ=255)
                                                                                    gCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTGGTTATAGGTGGCCCAAATGGTGAATTCGGTCTGACTGGTCGTAAAAATATAGTTGATAAAAACCGCGGCATGGCGCGTCACg                                                                 >  1:4595532/1‑149 (MQ=255)
                                                                                     ccAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACgg                                                                <  2:3938887/149‑1 (MQ=255)
                                                                                       aaTTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGGTATCGGTGGGCCCATGGGGGTCGGCGGTCTGTCTGGTCGTTAAAATATCGGTGATAACGACGGCGGCCTGGCGCGGCACGGTg                                                              >  1:3295817/1‑149 (MQ=255)
                                                                                         ttCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGGTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGc                                                            <  2:4217896/149‑1 (MQ=255)
                                                                                         gtCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGc                                                            <  2:3551315/148‑1 (MQ=255)
                                                                                          tCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCg                                                           <  2:1075665/149‑1 (MQ=255)
                                                                                          tCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAAATATAGTTGATAACTACGGCGGCATGGCGCGTCACGGTGGCg                                                           >  2:5185954/1‑149 (MQ=255)
                                                                                           cTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCgg                                                          <  1:1816915/149‑1 (MQ=255)
                                                                                            tGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTGGGTATGGGTGGGCCCAAGGGTGTCTGCGGTCTGAATTGTGGTAAAAATATATTTGATAACAAAAGCGGCAATGCGCGTCACCGGGGCGGt                                                         >  1:1780229/1‑149 (MQ=255)
                                                                                            tGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACGGGTCGTAAAATTATCGGTGATACCTACGGCGGCATGGCGCGGCACGGTGGCggg                                                         >  2:4100755/1‑148 (MQ=255)
                                                                                               ccGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCa                                                      >  2:796487/1‑149 (MQ=255)
                                                                                                cGCTGAATGGCTGACTTCTGCCACCAAATTCGTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCAt                                                     <  2:311031/149‑1 (MQ=255)
                                                                                                 gCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCAtt                                                    <  2:4408711/149‑1 (MQ=255)
                                                                                                  cTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGTCTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATtc                                                   >  2:1646842/1‑149 (MQ=255)
                                                                                                    gAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATtctc                                                 <  2:1349666/149‑1 (MQ=255)
                                                                                                    gAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGAATGGTCGTAAAATTATAGTTGATAACTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATtctc                                                 >  2:3345354/1‑149 (MQ=255)
                                                                                                    gAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCCATGGGTGACGGCGGGCGGACTGGTGGTAAAATTATCGTTGATAACTACGGCGGCATGGCGCGTCACCGGGGCGGTGGATtctc                                                 >  2:5029880/1‑149 (MQ=255)
                                                                                                     aaTGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTctct                                                <  1:1027826/149‑1 (MQ=255)
                                                                                                     aaTGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCAATGGGTGACTGCGGGCTGACTGGTCGTAAAATTATCGTTGATACCAACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTctct                                                >  2:1196600/1‑149 (MQ=255)
                                                                                                       tGGCTGACTTCTGCCACCAAATTCTTCCTCACCCCGACCGGTCGTGTCGTGATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTgg                                              <  1:5141559/149‑1 (MQ=255)
                                                                                                         gCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTa                                            <  1:4325968/149‑1 (MQ=255)
                                                                                                           tGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTaaa                                          <  1:1355357/149‑1 (MQ=255)
                                                                                                                 cTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCa                                    <  2:1569904/149‑1 (MQ=255)
                                                                                                                   gCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATAGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATc                                  >  2:4676642/1‑149 (MQ=255)
                                                                                                                     caccaAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaa                                <  1:341626/149‑1 (MQ=255)
                                                                                                                      accaAATTCTTCATCAACCCGACCGGTCGTTTAGTTATAGGTGGCCCAAATGGTGACATCGGTCTTAAAGGTCGTAAAATTATAGTTGATAACAAACGCCGCAAGGCGCGTCACGGTGGCGGTGCATTCTCTTGTAAAGATACaaaaaa                               >  1:2085743/1‑144 (MQ=255)
                                                                                                                        caAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAg                             <  2:28883/149‑1 (MQ=255)
                                                                                                                        caAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAg                             <  1:1973545/149‑1 (MQ=255)
                                                                                                                         aaaTTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGt                            <  1:1144395/149‑1 (MQ=255)
                                                                                                                         aaaTTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGt                            >  2:3410768/1‑149 (MQ=255)
                                                                                                                          aaTTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAAATATCGTTGATACCTACGGCGGCAAGGCGCGTCACGGTGGCGGTGCAATCTCTGGTAAAGATCCATCAAAAGTg                           >  2:4143037/1‑149 (MQ=255)
                                                                                                                          aaTTCTTCATCAACCCGACCGGTCGTTTCGTTATAGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTTATAACTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAATATCAATAAAAAGtt                           >  2:1984761/1‑148 (MQ=255)
                                                                                                                             tcttcATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGAc                        <  1:4383148/149‑1 (MQ=255)
                                                                                                                              cttcATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCCATGGGTGTCTGCGGTCTGACTGGTCGGAAAATTATAGTTGTTACCTAACGCGGCCTGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAATATACAACACAAGTGGAcc                       >  2:2282351/1‑149 (MQ=255)
                                                                                                                                tcatcaACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGt                     <  1:3218425/149‑1 (MQ=255)
                                                                                                                                 catcaACCCGACCGGTCGTTTCGTTATAGGTGGCCCCATGGGTGACTGCGGTCTGACGGGTCGTAAAATTATAGTTGATACCTACGGCGGCATGGCGCGGCACGGTGGCGGTGCATTATCTGGGAAAGATCCAACCAAAGTGGACCGtt                    >  2:4451767/1‑149 (MQ=255)
                                                                                                                                   tcaACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTcc                  <  2:1111753/149‑1 (MQ=255)
                                                                                                                                      aCCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCa               <  2:2555706/149‑1 (MQ=255)
                                                                                                                                       cccGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGTCTGGTCGTAAAATTATCGTTGATACCAACGGCGGCATTGCGCGTCACGGTGGCGGTGCATTCTATGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAg              >  2:1238800/1‑149 (MQ=255)
                                                                                                                                        ccGACCGGTCGTGTCGTTATCGGTGGCCCAAAGGGTGACGGCGGTCTGACTGGTCGTGAAAATATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGc             >  2:5198988/1‑149 (MQ=255)
                                                                                                                                         cGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGcc            <  1:654863/149‑1 (MQ=255)
                                                                                                                                         cGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTAGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTACGCAGcc            >  2:3463980/1‑149 (MQ=255)
                                                                                                                                                   tttcgttATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTCAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCAc  >  1:2694515/1‑149 (MQ=255)
                                                                                                                                                   tttcgttATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCAc  <  1:3522670/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCAC  >  CP009273/3080607‑3080902

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: