Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,080,755 T→G 22.4% 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,7550TG22.4% 267.8 / 49.6 192F231V (TTC→GTC) metKS‑adenosylmethionine synthetase
Reads supporting (aligned to +/- strand):  ref base T (55/91);  new base G (44/2);  total (99/93)
Fisher's exact test for biased strand distribution p-value = 4.08e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 208 total aligned reads displayed.

ATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCAC  >  CP009273/3080609‑3080902
                                                                                                                                                  |                                                                                                                                                   
  gCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttcg                                                                                                                                                  <  1:1290352/148‑1 (MQ=255)
  gCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttcg                                                                                                                                                  <  1:476905/148‑1 (MQ=255)
  gCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttcg                                                                                                                                                  <  1:515776/148‑1 (MQ=255)
  gCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttcg                                                                                                                                                  <  2:2879958/148‑1 (MQ=255)
  gCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttcg                                                                                                                                                  <  1:2955504/148‑1 (MQ=255)
    tGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtgtcgtt                                                                                                                                                >  1:3067895/1‑148 (MQ=255)
     gTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTa                                                                                                                                               >  2:1096149/1‑148 (MQ=255)
      tCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTAt                                                                                                                                              >  2:2079132/1‑148 (MQ=255)
      tCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTAt                                                                                                                                              >  2:267866/1‑148 (MQ=255)
      tCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTAt                                                                                                                                              >  2:2079144/1‑148 (MQ=255)
            tttCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTgg                                                                                                                                        >  1:2240374/1‑148 (MQ=255)
            tttCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGGTATAGGttg                                                                                                                                        >  1:2490177/1‑146 (MQ=255)
             ttCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGc                                                                                                                                       >  2:1994362/1‑148 (MQ=255)
             ttCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGc                                                                                                                                       >  2:2337117/1‑148 (MQ=255)
              tCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGcc                                                                                                                                      <  1:1742126/148‑1 (MQ=255)
                 aCTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCaa                                                                                                                                   >  2:385636/1‑148 (MQ=255)
                 aCTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTGGTTATAGGTGGcccca                                                                                                                                   >  1:2280158/1‑146 (MQ=255)
                 aCTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTTTAGGTGGcccca                                                                                                                                   >  1:1888539/1‑146 (MQ=255)
                  cTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAAt                                                                                                                                  >  2:2222603/1‑148 (MQ=255)
                    cAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATgg                                                                                                                                <  1:483265/148‑1 (MQ=255)
                       cACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTGGTTATCGGTGGCCCCAGGGGTg                                                                                                                             >  1:660495/1‑148 (MQ=255)
                          tctGAAGAGATCGACCAGAACTCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGCCCGGTCGGTTCGTTATACGGGGCCCACAGGGGGACt                                                                                                                          >  2:1025045/1‑148 (MQ=255)
                           ctGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGGTATCGGTGGCCCAATGGGTGACTg                                                                                                                         >  2:774340/1‑148 (MQ=255)
                              aaGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCgg                                                                                                                      >  2:1131125/1‑148 (MQ=255)
                               agagATCGACCAGACATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCCCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTGATCGGTGGCCCAATGGGTGACTGCGGt                                                                                                                     <  2:13331/148‑1 (MQ=255)
                               agagATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGt                                                                                                                     <  1:1743/148‑1 (MQ=255)
                               agagATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGt                                                                                                                     <  1:1752662/148‑1 (MQ=255)
                                 agaTCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCt                                                                                                                   >  2:2113481/1‑148 (MQ=255)
                                 agaTCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCt                                                                                                                   <  1:933807/148‑1 (MQ=255)
                                 agaTCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCt                                                                                                                   >  2:3651096/1‑148 (MQ=255)
                                  gaTCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTg                                                                                                                  <  2:594297/148‑1 (MQ=255)
                                   aTCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGa                                                                                                                 <  1:2042693/148‑1 (MQ=255)
                                   aTCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGa                                                                                                                 <  2:3623985/148‑1 (MQ=255)
                                        ccAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTGGTTTTCGGTGGCCCCATGGGTGTCTGCGGTCTGACTGGt                                                                                                            >  1:1494557/1‑148 (MQ=255)
                                         cAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTc                                                                                                           <  1:977312/148‑1 (MQ=255)
                                          aGAAATCGGTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGGCTGACTGGTCg                                                                                                          >  1:3766838/1‑148 (MQ=255)
                                              aTCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGTCTGGTCGGaaa                                                                                                      >  2:804131/1‑148 (MQ=255)
                                              aTCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTaaa                                                                                                      >  2:1778034/1‑148 (MQ=255)
                                                cGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTGGGTATAGGTGGCCCCCTGGGTGTCGGCGGTCTGTCAGGTCgggaaaa                                                                                                    >  1:750206/1‑142 (MQ=255)
                                                      aagaagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGt                                                                                              <  1:2582875/148‑1 (MQ=255)
                                                      aagaagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGt                                                                                              <  1:666476/148‑1 (MQ=255)
                                                       agaagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCCAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGTCGGGGCGTAAAATTATCGtt                                                                                             >  2:3760717/1‑148 (MQ=255)
                                                       agaagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGtt                                                                                             <  1:2376505/148‑1 (MQ=255)
                                                       agaagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGtt                                                                                             <  1:3488142/148‑1 (MQ=255)
                                                        gaagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTg                                                                                            <  2:1680813/148‑1 (MQ=255)
                                                        gaagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATGGGTGGCCCCATGGGTGTCTGCGGTGTGTCTGGTCGGTAAATTATAGTTg                                                                                            >  1:1701307/1‑148 (MQ=255)
                                                         aagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGa                                                                                           <  2:35555/148‑1 (MQ=255)
                                                         aagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGGCTGACTGGTCGTAAAATTATCGTTGa                                                                                           >  2:1442057/1‑148 (MQ=255)
                                                         aagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTAACGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGa                                                                                           >  2:2300808/1‑148 (MQ=255)
                                                         aagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGGTATAGGTGGCCCAAAGGGTGACTGCGGTCTGACTGGTCGTAAAATTATAGTTGa                                                                                           >  1:3785466/1‑148 (MQ=255)
                                                            cGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATAc                                                                                        <  2:2602677/148‑1 (MQ=255)
                                                            cGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATAc                                                                                        <  1:1758718/148‑1 (MQ=255)
                                                              gTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGGTATAGGTGGCCCCAAGGGTGTCTGGGGTCTGTCTGGGCGGAAAATTATCGTTGATACCt                                                                                      >  2:339303/1‑148 (MQ=255)
                                                               tAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTa                                                                                     <  1:1383663/148‑1 (MQ=255)
                                                                aaTGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTGGTTATAGGTGGGCCCAAGGGTGACGGCGGGCTGACTGGGCGTAAAATTATAGTTTATAACTAc                                                                                    >  1:3261470/1‑148 (MQ=255)
                                                                     aaGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTAcggcgg                                                                               <  1:3615675/148‑1 (MQ=255)
                                                                      agagATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACggcggc                                                                              <  2:311973/148‑1 (MQ=255)
                                                                      agagATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACggcggc                                                                              <  2:2867741/148‑1 (MQ=255)
                                                                      agagATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCGGACTGGGCGTAAAAATATCGTTGATAACTAAcgcggc                                                                              >  1:2168939/1‑148 (MQ=255)
                                                                      agagATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTGGTTATAGGTGGGCCCATGGGTGTCTGCGGTCTGCCTGGTGGTAAAATTATAGTTGATAACTAAggccgc                                                                              >  1:1079051/1‑148 (MQ=255)
                                                                       gagaTCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGGTATCGGTGGCCCAATGGGTGTCTGCGGGCTGACGGGTCGTAAAAATATCGTTGATACCTACGGCGGCa                                                                             >  1:2558376/1‑148 (MQ=255)
                                                                        agaTCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATGGGTGGCCCCATGGGTGACTGCGGTCTGACGGGTCGTAAAATTATAGTTTATACCTACggcggcct                                                                            >  2:884282/1‑146 (MQ=255)
                                                                        agaTCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCAt                                                                            >  2:2419920/1‑148 (MQ=255)
                                                                         gaTCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATg                                                                           <  1:2833311/148‑1 (MQ=255)
                                                                         gaTCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATg                                                                           <  2:1230037/148‑1 (MQ=255)
                                                                         gaTCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATg                                                                           <  1:629045/148‑1 (MQ=255)
                                                                               caAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTGGTTTTCGGTGGGCCAATGGGTGACGGCGGTCTGACTGGTCGGAAAATTATAGGTGTTAACGACCGCGGCATGGggcgg                                                                     >  1:1218618/1‑147 (MQ=255)
                                                                                aaGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATAGGTGGCCCAAAGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCTTTTATACATAAGGCGGCATGGCGCGTc                                                                    >  2:2926809/1‑148 (MQ=255)
                                                                                    cAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACgg                                                                <  1:2195649/148‑1 (MQ=255)
                                                                                       ttCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTGGTTTTGGGTGGGCCAAAGGGTGACAGCGGGGTGGCTGGTCGTTAAATTATAGTTTATAACGAAGGCGGCAATGCGCGGCACGGTgg                                                             >  1:751437/1‑148 (MQ=255)
                                                                                            cccGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATAGGTGGCCCAAAGGGTTACTGCGGTCTTACTGGTGGGAAAATTAACGTTGTTACCTACGGCGGCATGGCGCGGCCCGGTGGCGGTg                                                        >  1:2735267/1‑148 (MQ=255)
                                                                                            cccGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTGGTTTTGGGGGGGCCCAAGGGTGGCTGCGGTCTGACTGGTGGGTAAAATATAGTTTTTTACAACGGCGGCATGGCGCGGGACgggggggggg                                                        >  1:819803/1‑146 (MQ=255)
                                                                                              cGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCCATGGGTGACTGCGGTCTGACTGGTCGTCAAATTATCGTTGATAACTACGGCGGCCTGGCGCGGCACGGTGGCGGTGCa                                                      >  1:2403642/1‑148 (MQ=255)
                                                                                                cTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCAtt                                                    <  2:2848290/148‑1 (MQ=255)
                                                                                                cTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCAtt                                                    <  1:1994362/148‑1 (MQ=255)
                                                                                                     tGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTGGTTATCGGTGGGCCAATGGGTGACTGCGGTCTGACTGGTGGTAAAAATATCGTTTTTAACTACGGCGGCATGGCGCGTCAAGGGGGCGGTGCATTCTCTg                                               >  1:3444021/1‑148 (MQ=255)
                                                                                                     tGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTg                                               <  2:759864/148‑1 (MQ=255)
                                                                                                     tGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACAGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTg                                               >  2:2185626/1‑148 (MQ=255)
                                                                                                           aCTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAg                                         >  2:3459530/1‑148 (MQ=255)
                                                                                                            cTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGa                                        <  2:1890547/148‑1 (MQ=255)
                                                                                                              tCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATc                                      <  2:655595/148‑1 (MQ=255)
                                                                                                              tCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAAAGGGTGACTGCGGTCTGACTGGGCGTAAAATTATAGTTGATACCTAAGGCGGCATGGCGCGTCACCGTGGGGGTGGATTATCTGGTAAATATc                                      >  1:3182929/1‑148 (MQ=255)
                                                                                                               cTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATcc                                     <  2:1527556/148‑1 (MQ=255)
                                                                                                                 gCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCAt                                   <  2:3068079/148‑1 (MQ=255)
                                                                                                                   caccaAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATTGGTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCa                                 <  1:1467664/148‑1 (MQ=255)
                                                                                                                     ccaAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAAAGGGTGACTGCGGTGTGACTGGTCGTAAAATTATAGTTTATACCTAAGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa                               >  2:115518/1‑148 (MQ=255)
                                                                                                                      caAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaaa                              <  2:1625914/148‑1 (MQ=255)
                                                                                                                       aaaTTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAg                             <  1:2345793/148‑1 (MQ=255)
                                                                                                                       aaaTTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAg                             <  1:2344602/148‑1 (MQ=255)
                                                                                                                       aaaTTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGGGGCCCCATGGGTGTCTGCGGTCTGACTGGTGGGAAAATTATAGGTTATAACTACGGCGGCATTGGGCGTCACGGTGGCGGTGCATTCTATTGTAAATATACAaaaaaat                             >  1:2228110/1‑147 (MQ=255)
                                                                                                                        aaTTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGt                            <  2:3796737/148‑1 (MQ=255)
                                                                                                                        aaTTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGt                            <  2:2781089/148‑1 (MQ=255)
                                                                                                                          ttcttcATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGGCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTgg                          >  2:2877184/1‑148 (MQ=255)
                                                                                                                              tcatcaACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGTCTGCGGTCTGACTGGTCGTCAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGAccc                      >  2:2067690/1‑147 (MQ=255)
                                                                                                                                atcaACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGtt                    <  1:1546692/148‑1 (MQ=255)
                                                                                                                                    aCCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGGAAAAATATCGTTGATACCTAAGGCGGCATGGCGCGTCAAGGTGGCGGTGCATTCTCTGGTAAAGATCCAAAAAAAGTGGACCGTTCCGc                >  2:2700783/1‑148 (MQ=255)
                                                                                                                                       cGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGc             <  2:3162122/148‑1 (MQ=255)
                                                                                                                                          ccGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTTAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTa          <  2:3066605/148‑1 (MQ=255)
                                                                                                                                               cgtttcgttATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACgcag     >  1:3517108/1‑148 (MQ=255)
                                                                                                                                                gtttcgttATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACgcagc    <  1:2968735/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCAC  >  CP009273/3080609‑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: