Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,639,551 C→A 35.9% T269T (ACC→ACA rlmJ → 23S rRNA m(6)A2030 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,639,5510CA35.9% 473.7 / 299.9 300T269T (ACC→ACArlmJ23S rRNA m(6)A2030 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (95/97);  new base A (46/62);  total (141/159)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 302 total aligned reads displayed.

GCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATCGGCGGCGGCAGC  >  CP009273/3639404‑3639700
                                                                                                                                                   |                                                                                                                                                     
 cAAATTGCACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCggg                                                                                                                                                    <  2:1122975/150‑1 (MQ=255)
 cAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCggg                                                                                                                                                    <  2:3396791/150‑1 (MQ=255)
  aaaTTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGc                                                                                                                                                   >  2:3016010/1‑150 (MQ=255)
   aaTTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCa                                                                                                                                                  <  2:51872/150‑1 (MQ=255)
    aTTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCAc                                                                                                                                                 <  2:4166476/150‑1 (MQ=255)
    aTTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCAc                                                                                                                                                 <  2:609733/150‑1 (MQ=255)
    aTTGAACTGGCGGTACTGCCAGACAGCGATCGCCGAGCCATGACCGCTTCCGGCATGATTGAGATAAACCCGCCGTGGAAACTGGAACAACAGATGAAAAACGTGCTGCCGAGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGTGCAc                                                                                                                                                 >  2:359667/1‑150 (MQ=255)
       gAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGcc                                                                                                                                              >  2:4034953/1‑150 (MQ=255)
           tGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCg                                                                                                                                          >  2:1925850/1‑150 (MQ=255)
            ggCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGt                                                                                                                                         >  1:2711262/1‑150 (MQ=255)
            ggCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGt                                                                                                                                         <  2:3971129/150‑1 (MQ=255)
             tcGGTACTGCCAGAAAGCGATCGCCGTGGAATGACCGCTTCCGGCATGATTGTGATTAACCCTCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTTGCTGCACATCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCTTa                                                                                                                                        <  2:2305020/149‑1 (MQ=255)
             gCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTa                                                                                                                                        <  2:1523799/150‑1 (MQ=255)
               ggTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAg                                                                                                                                      <  1:2350554/150‑1 (MQ=255)
                 tACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCt                                                                                                                                    <  2:2400024/150‑1 (MQ=255)
                  aCTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCATGCACAGGGCACGCCACCGTAAGCTg                                                                                                                                   <  1:1077826/150‑1 (MQ=255)
                   cTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTgg                                                                                                                                  <  2:3949543/150‑1 (MQ=255)
                    tGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGa                                                                                                                                 <  2:3235653/150‑1 (MQ=255)
                    tGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGa                                                                                                                                 >  2:1319922/1‑150 (MQ=255)
                      ccAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATc                                                                                                                               <  2:2034336/150‑1 (MQ=255)
                       cAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCg                                                                                                                              >  2:1980031/1‑150 (MQ=255)
                         gACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTg                                                                                                                            >  2:2752243/1‑150 (MQ=255)
                         gACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTg                                                                                                                            >  2:1144097/1‑150 (MQ=255)
                           cAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGcc                                                                                                                          <  1:801318/150‑1 (MQ=255)
                            aGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGTCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCg                                                                                                                         <  1:1982859/150‑1 (MQ=255)
                                aTCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGt                                                                                                                     >  1:437055/1‑150 (MQ=255)
                                  cGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTaa                                                                                                                   <  1:2881312/150‑1 (MQ=255)
                                  cGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTaa                                                                                                                   >  2:2267046/1‑150 (MQ=255)
                                   gCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAAt                                                                                                                  >  2:160524/1‑150 (MQ=255)
                                    ccGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAAtt                                                                                                                 <  1:1002456/150‑1 (MQ=255)
                                     cGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTg                                                                                                                >  1:263243/1‑150 (MQ=255)
                                     cGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTg                                                                                                                >  2:1167965/1‑150 (MQ=255)
                                     cGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTg                                                                                                                >  2:2477723/1‑150 (MQ=255)
                                        ggCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAg                                                                                                             <  2:3966319/150‑1 (MQ=255)
                                          cATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGcc                                                                                                           <  2:3311889/150‑1 (MQ=255)
                                            tGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCAt                                                                                                         >  1:1329781/1‑150 (MQ=255)
                                            tGACCGCAACCGGCATGATTGTGATTAACCCGCCGTGGAAACAGGAACAACATATGAAAAACCTCCTCCCGTGGCTGCACACCAAACTCGTACCTGCATGCACCGGTCACGCCACAGTAAGCAGGATCGAGCCGGAGTAAATGCAGCCAt                                                                                                         >  1:540000/1‑150 (MQ=255)
                                                 gCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTg                                                                                                    <  1:3585813/150‑1 (MQ=255)
                                                  cTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTgg                                                                                                   >  2:1724207/1‑150 (MQ=255)
                                                   ttCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGc                                                                                                  >  2:2681495/1‑150 (MQ=255)
                                                     ccGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCAc                                                                                                >  2:3299823/1‑150 (MQ=255)
                                                        gCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTa                                                                                             <  1:2085593/150‑1 (MQ=255)
                                                          aTGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTAtt                                                                                           >  2:1559401/1‑150 (MQ=255)
                                                            gATTTTGATTAACCCGCCGTGGAAAATTGAACAACAGATGAATAAATTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCATAGTGCACGCCAGCGTAAGCTGGAACGTGCCTGAGTAATTGCAGACATTGCTGGCAACTATTAc                                                                                         <  2:3968471/150‑1 (MQ=255)
                                                             aTTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACg                                                                                        <  1:1833472/150‑1 (MQ=255)
                                                                 tgATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGtctc                                                                                    >  2:1661967/1‑150 (MQ=255)
                                                                  gATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTcg                                                                                   <  2:3469445/150‑1 (MQ=255)
                                                                    ttAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTcgcg                                                                                 <  1:1302936/150‑1 (MQ=255)
                                                                    ttAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTcgcg                                                                                 <  1:3235016/150‑1 (MQ=255)
                                                                      aaCCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCt                                                                               <  1:969313/150‑1 (MQ=255)
                                                                      aaCCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCt                                                                               >  1:2037967/1‑150 (MQ=255)
                                                                         ccgccgTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACa                                                                            <  1:2663308/150‑1 (MQ=255)
                                                                          cgccgTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACaa                                                                           <  1:1663577/150‑1 (MQ=255)
                                                                           gccgTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAAt                                                                          <  2:3509713/150‑1 (MQ=255)
                                                                           gccgTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAAt                                                                          >  1:3386609/1‑150 (MQ=255)
                                                                              gTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATcgc                                                                       <  1:3414986/150‑1 (MQ=255)
                                                                              gTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATcgc                                                                       >  1:2511759/1‑150 (MQ=255)
                                                                               tGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATcgcg                                                                      >  1:3791796/1‑150 (MQ=255)
                                                                                  aaaCTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTa                                                                   >  2:1907004/1‑150 (MQ=255)
                                                                                  aaaCTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTa                                                                   >  2:2956552/1‑150 (MQ=255)
                                                                                   aaCTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTaa                                                                  <  1:3685942/150‑1 (MQ=255)
                                                                                    aCTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAAt                                                                 <  2:1010968/150‑1 (MQ=255)
                                                                                    aCTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAAt                                                                 <  1:2423215/150‑1 (MQ=255)
                                                                                        gAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCAACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAAc                                                             >  1:355708/1‑150 (MQ=255)
                                                                                           caacaGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGAt                                                          >  2:2367730/1‑150 (MQ=255)
                                                                                            aacaGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATa                                                         >  1:2684817/1‑150 (MQ=255)
                                                                                             acaGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATaa                                                        <  2:1357917/150‑1 (MQ=255)
                                                                                               aGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAgg                                                      <  1:1479774/150‑1 (MQ=255)
                                                                                               aGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAgg                                                      <  1:2108761/150‑1 (MQ=255)
                                                                                                   gAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGacac                                                  >  1:2617301/1‑150 (MQ=255)
                                                                                                     aTAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACtt                                                >  1:2230096/1‑150 (MQ=255)
                                                                                                     aTAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACtt                                                <  2:2544492/150‑1 (MQ=255)
                                                                                                       aaCGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTg                                              <  2:2966408/150‑1 (MQ=255)
                                                                                                       aaCGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTg                                              >  1:1619849/1‑150 (MQ=255)
                                                                                                        aCGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGt                                             <  2:3700689/150‑1 (MQ=255)
                                                                                                            gctgcCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATg                                         <  2:1496288/150‑1 (MQ=255)
                                                                                                             ctgcCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGa                                        <  2:4111004/150‑1 (MQ=255)
                                                                                                               gcCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACt                                      <  1:3937651/150‑1 (MQ=255)
                                                                                                               gcCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCAGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACt                                      <  1:1911120/150‑1 (MQ=255)
                                                                                                                 cGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTaa                                    >  2:3859804/1‑150 (MQ=255)
                                                                                                                 cGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTaa                                    >  1:3082205/1‑150 (MQ=255)
                                                                                                                  gTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTaaa                                   <  2:2166964/150‑1 (MQ=255)
                                                                                                                  gTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTaaa                                   <  1:4046670/150‑1 (MQ=255)
                                                                                                                       tGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGTCACTTTGTCATGACTAAACACTa                              <  1:3027683/150‑1 (MQ=255)
                                                                                                                         cacaGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATg                            <  2:798055/150‑1 (MQ=255)
                                                                                                                         cacaGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATg                            >  1:3946111/1‑150 (MQ=255)
                                                                                                                           caGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGAt                          <  2:176404/150‑1 (MQ=255)
                                                                                                                               aaaCTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACa                      <  1:1521734/150‑1 (MQ=255)
                                                                                                                                aaCTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACAt                     <  2:703138/150‑1 (MQ=255)
                                                                                                                                aaCTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACAt                     <  1:2919254/150‑1 (MQ=255)
                                                                                                                                  cTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCg                   <  1:3137646/150‑1 (MQ=255)
                                                                                                                                  cTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCg                   >  1:3056996/1‑150 (MQ=255)
                                                                                                                                   tGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACGTCGc                  >  1:820866/1‑150 (MQ=255)
                                                                                                                                    ggTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGcc                 >  2:3195027/1‑150 (MQ=255)
                                                                                                                                      ttCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCAt               <  2:1872703/150‑1 (MQ=255)
                                                                                                                                      ttCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCAt               <  2:533288/150‑1 (MQ=255)
                                                                                                                                          ggcaggcaCCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATcggc           >  2:2931183/1‑150 (MQ=255)
                                                                                                                                          ggcaggcaCCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAAAGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATcggc           >  1:1248186/1‑150 (MQ=255)
                                                                                                                                          ggcaggcaCAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATcggc           <  1:3088479/150‑1 (MQ=255)
                                                                                                                                                caCCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATCggcggcggc     >  2:1407984/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATCGGCGGCGGCAGC  >  CP009273/3639404‑3639700

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: