Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,080,618 T→A 93.7% V185E (GTG→GAG)  metK → S‑adenosylmethionine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,080,6180TA93.7% 986.7 / 43.3 287V185E (GTG→GAG) metKS‑adenosylmethionine synthetase
Reads supporting (aligned to +/- strand):  ref base T (12/6);  new base A (128/141);  total (140/147)
Fisher's exact test for biased strand distribution p-value = 1.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

Only 100 of 287 total aligned reads displayed.

CAGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCG  >  CP009273/3080471‑3080764
                                                                                                                                                   |                                                                                                                                                  
cAGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGc                                                                                                                                                  >  2:3467656/1‑150 (MQ=255)
cAGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGc                                                                                                                                                  >  2:4671113/1‑150 (MQ=255)
 aGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCt                                                                                                                                                 >  1:2451803/1‑150 (MQ=255)
  gCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCtt                                                                                                                                                <  2:4669699/150‑1 (MQ=255)
   cacctatcacctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCttt                                                                                                                                               <  2:984736/150‑1 (MQ=255)
    acctatcacctatGCACCCCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTc                                                                                                                                              <  2:4703573/150‑1 (MQ=255)
    acctatcacctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTc                                                                                                                                              >  2:4213755/1‑150 (MQ=255)
      ctatcacctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCTGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCa                                                                                                                                            >  2:4461211/1‑150 (MQ=255)
       tatcacctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCAc                                                                                                                                           >  1:2844226/1‑150 (MQ=255)
        atcacctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACt                                                                                                                                          >  1:2010759/1‑150 (MQ=255)
        atcacctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACt                                                                                                                                          >  1:864314/1‑150 (MQ=255)
             ctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCa                                                                                                                                     <  2:3504624/150‑1 (MQ=255)
              tatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCCCTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCAc                                                                                                                                    >  1:3739898/1‑150 (MQ=255)
                 gCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCActct                                                                                                                                 <  1:4035643/150‑1 (MQ=255)
                  cacaCCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTg                                                                                                                                >  1:66157/1‑150 (MQ=255)
                      ccGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAaga                                                                                                                            >  1:2769239/1‑150 (MQ=255)
                        gTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTGTCCACTCAGCACTCTGTAgaga                                                                                                                          >  2:3871462/1‑150 (MQ=255)
                          cTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATc                                                                                                                        <  1:2194427/150‑1 (MQ=255)
                           tGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCg                                                                                                                       >  1:557521/1‑150 (MQ=255)
                            ggTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGa                                                                                                                      >  1:1687566/1‑150 (MQ=255)
                               aCAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCa                                                                                                                   <  2:219281/150‑1 (MQ=255)
                                 aGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGa                                                                                                                 >  2:3657829/1‑150 (MQ=255)
                                  gCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGaa                                                                                                                <  1:3644858/150‑1 (MQ=255)
                                        ggCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCt                                                                                                          <  2:2278855/150‑1 (MQ=255)
                                         gCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTg                                                                                                         <  2:3827857/150‑1 (MQ=255)
                                           tGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCa                                                                                                       <  1:4878938/150‑1 (MQ=255)
                                              aGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCaaga                                                                                                    <  2:141652/150‑1 (MQ=255)
                                              aGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCaaga                                                                                                    <  2:4127970/150‑1 (MQ=255)
                                               gTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCaagaa                                                                                                   <  1:5022620/150‑1 (MQ=255)
                                                tGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGAaagaag                                                                                                  <  1:1901598/150‑1 (MQ=255)
                                                 gCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGc                                                                                                 <  2:2677204/150‑1 (MQ=255)
                                                  cGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCg                                                                                                <  1:257108/150‑1 (MQ=255)
                                                   gTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCgg                                                                                               <  2:4681738/150‑1 (MQ=255)
                                                     aaaaaCGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTa                                                                                             >  2:3438414/1‑150 (MQ=255)
                                                     aaaaaCGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGGATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAACTCGCTGCACGACGCGGTa                                                                                             >  2:4592501/1‑150 (MQ=255)
                                                       aaaCGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAAt                                                                                           <  2:287307/150‑1 (MQ=255)
                                                       aaaCGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAAt                                                                                           <  1:1436404/150‑1 (MQ=255)
                                                          cGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGa                                                                                        >  2:1695335/1‑150 (MQ=255)
                                                              aCTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAagag                                                                                    >  1:3776110/1‑150 (MQ=255)
                                                              aCTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAagag                                                                                    >  2:842909/1‑150 (MQ=255)
                                                                tctGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGat                                                                                  <  1:4135625/150‑1 (MQ=255)
                                                                tctGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGat                                                                                  <  1:4646962/150‑1 (MQ=255)
                                                                  tGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGatca                                                                                >  2:304262/1‑150 (MQ=255)
                                                                  tGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGatca                                                                                >  1:1089761/1‑150 (MQ=255)
                                                                  tGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGatca                                                                                <  1:3295555/150‑1 (MQ=255)
                                                                   gCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGatcat                                                                               >  1:3931626/1‑150 (MQ=255)
                                                                      gTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCaa                                                                            <  2:1466921/150‑1 (MQ=255)
                                                                      gTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCaa                                                                            >  1:4957998/1‑150 (MQ=255)
                                                                      gTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCaa                                                                            <  1:1528413/150‑1 (MQ=255)
                                                                      gTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCaa                                                                            <  2:1184432/150‑1 (MQ=255)
                                                                        ggCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGc                                                                          >  1:4732065/1‑150 (MQ=255)
                                                                         gCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGcc                                                                         >  2:4274222/1‑150 (MQ=255)
                                                                              gcCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTc                                                                    <  2:1306917/150‑1 (MQ=255)
                                                                                 cGGACGCGAAAAGCCAGGTGACTTTTCAGTATGCCGACGGCAAAATCGTGGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGCTCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGc                                                                 <  2:743161/150‑1 (MQ=255)
                                                                                         aaaaGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGaa                                                         >  1:1021245/1‑150 (MQ=255)
                                                                                         aaaaGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGaa                                                         >  2:1238464/1‑150 (MQ=255)
                                                                                          aaaGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAAt                                                        <  1:3602243/150‑1 (MQ=255)
                                                                                           aaGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATg                                                       >  1:2266225/1‑150 (MQ=255)
                                                                                           aaGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATg                                                       <  2:3595508/150‑1 (MQ=255)
                                                                                            aGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATgg                                                      >  2:2242338/1‑150 (MQ=255)
                                                                                            aGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATgg                                                      <  2:538013/150‑1 (MQ=255)
                                                                                             gCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGc                                                     <  1:4221361/150‑1 (MQ=255)
                                                                                              ccAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCt                                                    <  2:3697310/150‑1 (MQ=255)
                                                                                                aGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGa                                                  >  2:4166919/1‑150 (MQ=255)
                                                                                                 ggTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGAc                                                 >  2:4845508/1‑150 (MQ=255)
                                                                                                   tGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACtt                                               <  1:1471805/150‑1 (MQ=255)
                                                                                                     aCTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAACTCGCTGCAAGAAGCGGTCATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCt                                             >  2:4528855/1‑150 (MQ=255)
                                                                                                     aCTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCACGAAGCGGTAAGGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCt                                             >  1:817399/1‑150 (MQ=255)
                                                                                                     aCTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCt                                             <  1:1719369/150‑1 (MQ=255)
                                                                                                      cTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTg                                            >  1:3849793/1‑150 (MQ=255)
                                                                                                        tttCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAAGTCTGCCCGCTGAATGGCTGACTTCTGcc                                          >  2:3109823/1‑150 (MQ=255)
                                                                                                           cAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGccacc                                       >  1:1248131/1‑150 (MQ=255)
                                                                                                              tATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCaaa                                    >  2:4322501/1‑150 (MQ=255)
                                                                                                                 gacgacGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAAttc                                 >  2:1378754/1‑150 (MQ=255)
                                                                                                                 gacgacGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAAttc                                 <  2:3658704/150‑1 (MQ=255)
                                                                                                                  acgacgGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAAttct                                <  1:4801180/150‑1 (MQ=255)
                                                                                                                  acgacgGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAAttct                                <  1:940444/150‑1 (MQ=255)
                                                                                                                    gacgGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAAttcttc                              >  2:1862855/1‑150 (MQ=255)
                                                                                                                    gacgGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAAttcttc                              <  1:4274222/150‑1 (MQ=255)
                                                                                                                        gCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTtcatca                          >  2:1099718/1‑150 (MQ=255)
                                                                                                                         cAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCaa                         <  2:1779429/150‑1 (MQ=255)
                                                                                                                         cAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCaa                         >  2:1234161/1‑150 (MQ=255)
                                                                                                                          aaaaTCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAAc                        <  2:2453921/150‑1 (MQ=255)
                                                                                                                           aaaTCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAAcc                       >  2:2623815/1‑150 (MQ=255)
                                                                                                                           aaaTCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAAcc                       <  1:4528070/150‑1 (MQ=255)
                                                                                                                             aTCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCg                     >  1:1346382/1‑150 (MQ=255)
                                                                                                                              tCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGa                    >  2:3313576/1‑150 (MQ=255)
                                                                                                                              tCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGa                    >  2:214675/1‑150 (MQ=255)
                                                                                                                               cGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGAc                   >  2:2726466/1‑150 (MQ=255)
                                                                                                                                gTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGAcc                  <  1:304262/150‑1 (MQ=255)
                                                                                                                                 ttGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGACGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCCCCAAATTCTTCATCAACCCGACCg                 >  1:3395179/1‑150 (MQ=255)
                                                                                                                                 ttGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCg                 <  1:3170855/150‑1 (MQ=255)
                                                                                                                                 ttGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCg                 <  1:1199977/150‑1 (MQ=255)
                                                                                                                                 ttGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCg                 <  1:3170222/150‑1 (MQ=255)
                                                                                                                                        cGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgttt          >  2:2945480/1‑150 (MQ=255)
                                                                                                                                         gATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttc         >  1:3232101/1‑150 (MQ=255)
                                                                                                                                         gATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttc         >  1:2292543/1‑150 (MQ=255)
                                                                                                                                           tGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttcgt       >  1:4623463/1‑150 (MQ=255)
                                                                                                                                            gCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttcgtt      <  1:2123402/150‑1 (MQ=255)
                                                                                                                                               gTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATc   >  1:2902136/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CAGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCG  >  CP009273/3080471‑3080764

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