Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,639,551 C→A 10.3% 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,5510CA10.3% 638.8 / 33.4 183T269T (ACC→ACArlmJ23S rRNA m(6)A2030 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (65/99);  new base A (11/8);  total (76/107)
Fisher's exact test for biased strand distribution p-value = 1.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 185 total aligned reads displayed.

GCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATCGGCGG  >  CP009273/3639404‑3639693
                                                                                                                                                   |                                                                                                                                              
gCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCAcc                                                                                                                                                >  1:1596130/1‑148 (MQ=255)
          cTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCa                                                                                                                                      <  2:93785/148‑1 (MQ=255)
             gCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCg                                                                                                                                   <  1:2919928/148‑1 (MQ=255)
                 tACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTTCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAg                                                                                                                               >  1:3036646/1‑148 (MQ=255)
                 tACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAg                                                                                                                               <  1:2277002/148‑1 (MQ=255)
                  aCTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGc                                                                                                                              <  1:508380/148‑1 (MQ=255)
                   cTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCt                                                                                                                             <  1:3023216/148‑1 (MQ=255)
                     gCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTgg                                                                                                                           <  1:2544490/148‑1 (MQ=255)
                     gCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTgg                                                                                                                           <  1:1474500/148‑1 (MQ=255)
                      ccAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGa                                                                                                                          <  1:594713/148‑1 (MQ=255)
                      ccAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGa                                                                                                                          <  1:594722/148‑1 (MQ=255)
                       cAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGAt                                                                                                                         <  1:91740/148‑1 (MQ=255)
                         gACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCg                                                                                                                       >  2:1224845/1‑148 (MQ=255)
                             gCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGcc                                                                                                                   >  1:1710066/1‑148 (MQ=255)
                               gATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATATCGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCgg                                                                                                                 <  1:2098144/148‑1 (MQ=255)
                                aTCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGa                                                                                                                <  1:2809114/148‑1 (MQ=255)
                                   gCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTa                                                                                                             <  2:2203059/148‑1 (MQ=255)
                                   gCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTa                                                                                                             <  2:2202574/148‑1 (MQ=255)
                                         gCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCa                                                                                                       >  2:2985194/1‑148 (MQ=255)
                                          cATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAg                                                                                                      <  1:662205/148‑1 (MQ=255)
                                           aTGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGc                                                                                                     >  1:2742542/1‑148 (MQ=255)
                                              aCCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCAt                                                                                                  >  2:1892886/1‑148 (MQ=255)
                                               ccGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCAtt                                                                                                 <  1:2399061/148‑1 (MQ=255)
                                               ccGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCAtt                                                                                                 <  1:1987600/148‑1 (MQ=255)
                                                 gCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGc                                                                                               <  1:1485020/148‑1 (MQ=255)
                                                  cTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCt                                                                                              <  1:2436997/148‑1 (MQ=255)
                                                   ttCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTg                                                                                             <  2:2911732/148‑1 (MQ=255)
                                                    tCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTgg                                                                                            >  1:1802592/1‑148 (MQ=255)
                                                      cGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCa                                                                                          <  1:815161/148‑1 (MQ=255)
                                                      cGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCa                                                                                          <  1:2511114/148‑1 (MQ=255)
                                                           tGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTAt                                                                                     >  1:1375027/1‑148 (MQ=255)
                                                             aTTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTa                                                                                   >  2:2248631/1‑148 (MQ=255)
                                                                gtgATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGt                                                                                <  2:2413285/148‑1 (MQ=255)
                                                                gtgATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGt                                                                                <  1:1577979/148‑1 (MQ=255)
                                                                  gATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGtct                                                                              >  1:1357987/1‑148 (MQ=255)
                                                                      aaCCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTcgcg                                                                          >  2:2440511/1‑148 (MQ=255)
                                                                        cccGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCt                                                                        >  2:1216625/1‑148 (MQ=255)
                                                                        cccGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCt                                                                        >  2:1307789/1‑148 (MQ=255)
                                                                         ccgccgTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTa                                                                       <  1:103859/148‑1 (MQ=255)
                                                                         ccgccgTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTa                                                                       <  1:2893670/148‑1 (MQ=255)
                                                                         ccgccgTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTa                                                                       <  1:103189/148‑1 (MQ=255)
                                                                            ccgTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACaa                                                                    <  2:2143527/148‑1 (MQ=255)
                                                                            ccgTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACaa                                                                    >  2:2556200/1‑148 (MQ=255)
                                                                             cgTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAAt                                                                   <  1:1437179/148‑1 (MQ=255)
                                                                             cgTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAAt                                                                   >  1:144255/1‑148 (MQ=255)
                                                                                   aaCTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGt                                                             <  2:343669/148‑1 (MQ=255)
                                                                                    aCTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTa                                                            >  2:1757745/1‑148 (MQ=255)
                                                                                      tGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAAt                                                          <  1:1093131/148‑1 (MQ=255)
                                                                                        gAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACACTCGCGGTAATCa                                                        <  1:3178348/148‑1 (MQ=255)
                                                                                        gAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCa                                                        >  1:3173964/1‑148 (MQ=255)
                                                                                        gAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCa                                                        >  1:3204489/1‑148 (MQ=255)
                                                                                        gAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCa                                                        <  2:2388589/148‑1 (MQ=255)
                                                                                        gAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCa                                                        <  1:661193/148‑1 (MQ=255)
                                                                                         aacaacAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCaa                                                       >  1:1813018/1‑148 (MQ=255)
                                                                                         aacaacAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCaa                                                       <  2:2706664/148‑1 (MQ=255)
                                                                                          acaacaGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAAc                                                      >  1:1937711/1‑148 (MQ=255)
                                                                                           caacaGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACg                                                     >  2:262410/1‑148 (MQ=255)
                                                                                            aacaGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGa                                                    <  1:2391086/148‑1 (MQ=255)
                                                                                            aacaGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGa                                                    <  1:3093882/148‑1 (MQ=255)
                                                                                             acaGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGAt                                                   >  2:1506411/1‑148 (MQ=255)
                                                                                              caGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATa                                                  >  1:395200/1‑148 (MQ=255)
                                                                                              caGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATa                                                  <  1:247363/148‑1 (MQ=255)
                                                                                               aGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATaa                                                 <  2:1231776/148‑1 (MQ=255)
                                                                                                 aTGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAgg                                               >  1:1981013/1‑148 (MQ=255)
                                                                                                   gAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGac                                             <  2:1710066/148‑1 (MQ=255)
                                                                                                      tAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACt                                          <  2:3285318/148‑1 (MQ=255)
                                                                                                       aaCGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACtt                                         >  1:2232022/1‑148 (MQ=255)
                                                                                                       aaCGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACtt                                         <  2:2967065/148‑1 (MQ=255)
                                                                                                       aaCGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACtt                                         <  2:668578/148‑1 (MQ=255)
                                                                                                        aCGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACttt                                        >  2:285303/1‑148 (MQ=255)
                                                                                                        aCGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACttt                                        >  2:740563/1‑148 (MQ=255)
                                                                                                         cGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTg                                       <  1:2939816/148‑1 (MQ=255)
                                                                                                           tgctgcCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTc                                     >  2:1589956/1‑148 (MQ=255)
                                                                                                              tgcCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATg                                  <  1:1272563/148‑1 (MQ=255)
                                                                                                              tgcCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATg                                  <  1:2921540/148‑1 (MQ=255)
                                                                                                                ccGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGAc                                >  2:1774017/1‑148 (MQ=255)
                                                                                                                ccGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACAGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGAc                                >  2:1773662/1‑148 (MQ=255)
                                                                                                                   tGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTaa                             >  1:270881/1‑148 (MQ=255)
                                                                                                                    ggCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTaaa                            <  1:1913215/148‑1 (MQ=255)
                                                                                                                     gCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAac                           >  1:343929/1‑148 (MQ=255)
                                                                                                                      cTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAaca                          <  1:1506411/148‑1 (MQ=255)
                                                                                                                         cacaGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTa                       <  1:1335004/148‑1 (MQ=255)
                                                                                                                         cacaGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTa                       <  1:1265328/148‑1 (MQ=255)
                                                                                                                         cacaGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGCCATGACTAAACACTa                       <  2:88295/148‑1 (MQ=255)
                                                                                                                           caGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATg                     <  1:326884/148‑1 (MQ=255)
                                                                                                                            aGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGa                    <  1:1217711/148‑1 (MQ=255)
                                                                                                                              cAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGAtt                  >  2:2280029/1‑148 (MQ=255)
                                                                                                                              cAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGAtt                  >  1:2478293/1‑148 (MQ=255)
                                                                                                                              cAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGAtt                  >  1:2254481/1‑148 (MQ=255)
                                                                                                                               aaaCTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTa                 <  1:2960196/148‑1 (MQ=255)
                                                                                                                                aaCTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTAc                <  2:563610/148‑1 (MQ=255)
                                                                                                                                aaCTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGAGTAc                <  2:563624/148‑1 (MQ=255)
                                                                                                                                    ggTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCg            >  2:401699/1‑148 (MQ=255)
                                                                                                                                     gTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGc           >  1:1456781/1‑148 (MQ=255)
                                                                                                                                      ttCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGcc          >  2:3008374/1‑148 (MQ=255)
                                                                                                                                      ttCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGcc          <  2:191595/148‑1 (MQ=255)
                                                                                                                                       tCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCa         <  2:734193/148‑1 (MQ=255)
                                                                                                                                        ccGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCAt        <  2:457916/148‑1 (MQ=255)
                                                                                                                                           gcaggcaCCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATcgg     <  2:2548150/148‑1 (MQ=255)
                                                                                                                                              ggcaCCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATcggcgg  >  2:1980657/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
GCAAATTGAACTGGCGGTACTGCCAGACAGCGATCGCCGTGGCATGACCGCTTCCGGCATGATTGTGATTAACCCGCCGTGGAAACTGGAACAACAGATGAATAACGTGCTGCCGTGGCTGCACAGCAAACTGGTTCCGGCAGGCACCGGGCACGCCACCGTAAGCTGGATCGTGCCGGAGTAATTGCAGCCATTGCTGGCACCTATTACGTCTCGCGCTACAATCGCGGTAATCAACGATAAGGACACTTTGTCATGACTAAACACTATGATTACATCGCCATCGGCGG  >  CP009273/3639404‑3639693

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: