Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,080,755 T→G 20.5% 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.5% 491.1 / 70.3 263F231V (TTC→GTC) metKS‑adenosylmethionine synthetase
Reads supporting (aligned to +/- strand):  ref base T (45/160);  new base G (56/2);  total (101/162)
Fisher's exact test for biased strand distribution p-value = 1.35e-26
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 288 total aligned reads displayed.

TCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACG  >  CP009273/3080606‑3080903
                                                                                                                                                     |                                                                                                                                                    
 cGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgttt                                                                                                                                                     <  2:2767717/150‑1 (MQ=255)
 cGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgttt                                                                                                                                                     <  1:1883557/150‑1 (MQ=255)
 cGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgttg                                                                                                                                                     >  1:413413/1‑149 (MQ=255)
  gATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttc                                                                                                                                                    <  1:176629/150‑1 (MQ=255)
  gATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtgtc                                                                                                                                                    >  1:4207979/1‑150 (MQ=255)
  gATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtgtc                                                                                                                                                    >  2:1152203/1‑150 (MQ=255)
    tGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttcgt                                                                                                                                                  <  2:3664639/150‑1 (MQ=255)
     gCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttcggt                                                                                                                                                 >  1:1523024/1‑148 (MQ=255)
      cTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGGTa                                                                                                                                                >  1:5234948/1‑150 (MQ=255)
       tGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTAt                                                                                                                                               <  1:2842211/150‑1 (MQ=255)
        gTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATc                                                                                                                                              <  2:4358284/150‑1 (MQ=255)
            tGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTg                                                                                                                                          <  2:4175476/150‑1 (MQ=255)
             gCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTgg                                                                                                                                         <  2:4507770/150‑1 (MQ=255)
                   cACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAAt                                                                                                                                   <  2:4262951/150‑1 (MQ=255)
                       cAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGt                                                                                                                               <  2:1774732/150‑1 (MQ=255)
                        aGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTg                                                                                                                              <  2:3750784/150‑1 (MQ=255)
                            ctctGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGAATg                                                                                                                          >  2:3899154/1‑150 (MQ=255)
                             tctGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGc                                                                                                                         <  1:3853621/150‑1 (MQ=255)
                                  agagATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCt                                                                                                                    <  2:1384884/150‑1 (MQ=255)
                                  agagATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCt                                                                                                                    <  2:37011/150‑1 (MQ=255)
                                  agagATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCt                                                                                                                    <  1:5116848/150‑1 (MQ=255)
                                  agagATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCt                                                                                                                    <  1:1094711/150‑1 (MQ=255)
                                    agaTCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGa                                                                                                                  <  1:3068879/150‑1 (MQ=255)
                                     gaTCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGAc                                                                                                                 <  1:4360834/150‑1 (MQ=255)
                                      aTCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACt                                                                                                                <  1:3741177/150‑1 (MQ=255)
                                         gACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGt                                                                                                             <  2:4309267/150‑1 (MQ=255)
                                            cAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGt                                                                                                          <  1:695892/150‑1 (MQ=255)
                                            cAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGt                                                                                                          <  2:4588972/150‑1 (MQ=255)
                                            cAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGt                                                                                                          <  2:4589856/150‑1 (MQ=255)
                                             aGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTGGTTATCGGTGGCCCAAGGGGGGGCTGCGGGCTGACTGGTCggt                                                                                                         >  1:1928914/1‑148 (MQ=255)
                                              gAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTaa                                                                                                        <  2:226358/150‑1 (MQ=255)
                                                aaTCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTaaaa                                                                                                      <  1:1661311/150‑1 (MQ=255)
                                                aaTCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCCATGGGTGACTGCGGTCTGACTGGTCGTaaaa                                                                                                      >  2:658156/1‑150 (MQ=255)
                                                 aTCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAAt                                                                                                     <  2:4932061/150‑1 (MQ=255)
                                                 aTCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAAt                                                                                                     <  2:3029211/150‑1 (MQ=255)
                                                      tGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGGGGGCCCATGGGTGTCGGCGGGGTGTCTGGTGGTAAAATTAtat                                                                                                >  1:992087/1‑148 (MQ=255)
                                                        cAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGtt                                                                                              <  2:3444734/150‑1 (MQ=255)
                                                        cAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATAGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGtt                                                                                              >  2:2148318/1‑150 (MQ=255)
                                                         aagaagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTAGTTATAGGTGGCCCAAAGGGTGAAAGCGGTCTGACTGGTCGTAAAATTATAGTTg                                                                                             >  2:64289/1‑150 (MQ=255)
                                                         aagaagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGGTATAGGGGGCCCCAACGGTGACTGCGGGCGGACTGGGCGTAAAATTATAGTTg                                                                                             >  1:4503806/1‑150 (MQ=255)
                                                          agaagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGa                                                                                            <  2:1457186/150‑1 (MQ=255)
                                                           gaagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGAt                                                                                           <  1:4606533/150‑1 (MQ=255)
                                                           gaagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGAt                                                                                           >  1:4805061/1‑150 (MQ=255)
                                                            aagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATa                                                                                          <  1:1987714/150‑1 (MQ=255)
                                                            aagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTGGTTATAGGTGGCCCACATGGTGTGTGGGGTCTGGCTGGGGGTTAAATTATAGtttttt                                                                                          >  1:3560234/1‑146 (MQ=255)
                                                              gCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATAcc                                                                                        <  1:3236879/150‑1 (MQ=255)
                                                               cGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTGGTTTTCGGTGGCCCAATTGGTGTCTGCGGTCTGAAGTGTCGTGAAATTATAGTTGTTaacg                                                                                       >  1:2470656/1‑147 (MQ=255)
                                                                ggTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATAGGTGGCCCCATGGGTGACTGCGGTCTGGCTGGTCGTAAAATTATCGTTGATaaata                                                                                      >  2:3606338/1‑146 (MQ=255)
                                                                 gTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTAc                                                                                     <  1:2773293/150‑1 (MQ=255)
                                                                  tAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTAcg                                                                                    <  1:2094356/150‑1 (MQ=255)
                                                                  tAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGGTATCGGTGGCCCAATGGGTGACTGCGGGCTGACTGGTCGTAAAATTATCGTTGATACCTAcg                                                                                    >  2:1959481/1‑150 (MQ=255)
                                                                   aaTGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGTCTGGTGGTGAAATTATCGTTGATACCTAcgg                                                                                   >  2:2954940/1‑150 (MQ=255)
                                                                        aaGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCa                                                                              <  1:4278435/150‑1 (MQ=255)
                                                                           agaTCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATgg                                                                           <  1:2155341/150‑1 (MQ=255)
                                                                             atcatcAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGgcg                                                                         <  2:2302873/150‑1 (MQ=255)
                                                                                  caAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCa                                                                    <  1:5007625/150‑1 (MQ=255)
                                                                                   aaGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTTTAGGTGGGCCAAAGGGTGTCTGCGGTCTGGATGGTCGGAAAAATATAGTTGGTAAAGACCGCGGCAAGGCGCGTCAc                                                                   >  1:4043125/1‑150 (MQ=255)
                                                                                   aaGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCAc                                                                   <  1:1116820/150‑1 (MQ=255)
                                                                                    aGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGGGGCCCCATGGGTGACTGCGGTCTGACGGGTCGTAAAATTATCGGTGATACCTACGGCGGCATGGCGCGTCCCg                                                                  >  2:1051894/1‑150 (MQ=255)
                                                                                     gCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTGGTTATAGGTGGCCCAAAGGGTGTCTGCGGTCTGACTGGTCGTAAAATTATAGTTGATAACAACGGCGGCATGGCGCGTCACgg                                                                 >  1:3041469/1‑150 (MQ=255)
                                                                                       cAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTg                                                               <  1:598902/150‑1 (MQ=255)
                                                                                              gCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGGAAAATTATCGTTGAAACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGc                                                        >  2:1696702/1‑150 (MQ=255)
                                                                                              gCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATAGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGc                                                        >  2:732499/1‑150 (MQ=255)
                                                                                                 cGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCAtt                                                     <  2:4988312/150‑1 (MQ=255)
                                                                                                  gCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATtc                                                    <  1:3354174/150‑1 (MQ=255)
                                                                                                   cTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTGGGTATGGGTGGGCCCATGGGTGTGTGCGGTCTTTCTGGTGGTTAAAATATAGTTGATAACAAAGGCGGCATGGCGCGTGACCGTGGCGGTGCAttat                                                   >  1:2371089/1‑148 (MQ=255)
                                                                                                   cTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTAAGGCGGCATGGCGCGTCACGGTGGCGGTGCAttat                                                   >  2:4285248/1‑148 (MQ=255)
                                                                                                     gAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTctct                                                 <  1:4006954/150‑1 (MQ=255)
                                                                                                           cTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTaaa                                           >  2:5197892/1‑150 (MQ=255)
                                                                                                            tGACTTCTGCCACCAAATTCGTCATCAACCCGACCGGTCGTTTCGGTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAg                                          <  2:1554320/150‑1 (MQ=255)
                                                                                                              aCTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGAt                                        <  2:2359190/150‑1 (MQ=255)
                                                                                                              aCTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGAt                                        <  1:747536/150‑1 (MQ=255)
                                                                                                              aCTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGAt                                        <  1:2667705/150‑1 (MQ=255)
                                                                                                               cTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATc                                       <  2:3737272/150‑1 (MQ=255)
                                                                                                                  cTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTGGTTATAGGTGGGCCAATGGGTTACTTCGGTGTGACTTGTCGTAAAAATATAGTTTTTAAATACGGCGGCATTGCGCGTCACCGTTGCGGTTCATTCTCTTGTAAATATACaa                                    >  1:5154474/1‑149 (MQ=255)
                                                                                                                      caccaAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa                                <  1:2740020/150‑1 (MQ=255)
                                                                                                                        ccaAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAg                              <  2:528697/150‑1 (MQ=255)
                                                                                                                        ccaAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAg                              <  2:4049012/150‑1 (MQ=255)
                                                                                                                        ccaAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAg                              <  1:2501953/150‑1 (MQ=255)
                                                                                                                           aaTTCTTCATCAACCCGACCGGTCGTGTCGTTATAGGTGGGCCCAAGGGTGTCTGCGGTCTGGCTGGTCGTAAAATTATAGTTTTTAACAAAGGCGGCCAGGCGCGGCACCGGGGCGGTGCAATCTCTGGTAAAGATACAAAAAAAGttg                           >  2:935932/1‑148 (MQ=255)
                                                                                                                            aTTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGa                          <  2:4104375/150‑1 (MQ=255)
                                                                                                                            aTTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGa                          <  1:1195935/150‑1 (MQ=255)
                                                                                                                             ttcttcATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGAc                         <  1:2815131/150‑1 (MQ=255)
                                                                                                                              tcttcATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGAcc                        <  1:836862/150‑1 (MQ=255)
                                                                                                                               cttcATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCg                       <  2:3319677/150‑1 (MQ=255)
                                                                                                                               cttcATCAACCCGACCGGTCGTTTCGTTATAGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTAAAGTTTATACCAAAGGCGGCATGGCGCGTCACCGTGGCGGGTCAATATATTGTAAAGATCCATAAAAATTTGACCg                       >  2:5019398/1‑150 (MQ=255)
                                                                                                                                  catcaACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                    <  2:476746/150‑1 (MQ=255)
                                                                                                                                  catcaACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                    <  1:3730854/150‑1 (MQ=255)
                                                                                                                                   atcaACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTcc                   <  2:4908300/150‑1 (MQ=255)
                                                                                                                                      aaCCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGCTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCa                <  2:1147686/150‑1 (MQ=255)
                                                                                                                                       aCCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACCGTGGCGGTGCATTCTCTGGTAAAGATCCAACAAAAGTGGACCGTTCCGCAg               >  2:857620/1‑150 (MQ=255)
                                                                                                                                         ccGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGcc             <  1:3674072/150‑1 (MQ=255)
                                                                                                                                          cGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCt            <  2:68229/150‑1 (MQ=255)
                                                                                                                                          cGACCGGTCGTGTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCt            <  2:3113036/150‑1 (MQ=255)
                                                                                                                                           gACCGGTCGTTTCGTTATAGGTGGCCCAAAGGGTGACAGCGGTCTGACTGGTCGTAAAAATAACGTTGATACCAAACGCCGCAAGGCGCGTCACCGTGGCCGTGCAATAACAGGTAAAAAAACAACAAAAGTTGACCCTTACGcccccaa           >  2:2606265/1‑144 (MQ=255)
                                                                                                                                             ccGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACg         <  1:472365/150‑1 (MQ=255)
                                                                                                                                             ccGGTCGTGTCGTTATAGGTGGGCCAATGGGTGACTGCCGTCTTAATGGTCCTAAAATTATAGTTTATAACAACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTAACAGGTAAAGATCCATAAAAAGTGGACCGTTCCGCAGCCAACg         >  2:747904/1‑150 (MQ=255)
                                                                                                                                               ggTCGTTTCCTTATCCGTGGACCAATGGGTGACTTCGGTCTGACAGGTCGTAAAATAATCGTTGATACCAACGGCGGCATGGCGCGTCACGGTGGCGGTGCAATCTCTGGTAAAGATCCATCAAAAGAGGACCGTTCCGCAGCCAAAgca       >  2:5110908/1‑150 (MQ=255)
                                                                                                                                                   gtttcgttATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCAc   >  2:1324801/1‑150 (MQ=255)
                                                                                                                                                   gtttcgttATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCAc   >  2:3161238/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACG  >  CP009273/3080606‑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: