Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 872,975 C→A 22.5% E173* (GAA→TAA)  rimO ← ribosomal protein S12 methylthiotransferase; radical SAM superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273872,9750CA22.5% 496.1 / 99.4 203E173* (GAA→TAA) rimOribosomal protein S12 methylthiotransferase; radical SAM superfamily
Reads supporting (aligned to +/- strand):  ref base C (69/88);  new base A (25/21);  total (94/109)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 203 total aligned reads displayed.

GCTGACCATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAGGAATG  >  CP009273/872829‑873121
                                                                                                                                                  |                                                                                                                                                  
gCTGACCATGCTTGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCg                                                                                                                                                   <  1:1334695/148‑1 (MQ=255)
gCTGACCATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATTGGCGTAAGTATCCTGCGAGATCACCAGAATCTATTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCg                                                                                                                                                   <  1:321244/148‑1 (MQ=255)
      cATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATc                                                                                                                                             >  1:236088/1‑148 (MQ=255)
       aTGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCg                                                                                                                                            >  2:3350838/1‑148 (MQ=255)
       aTGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCg                                                                                                                                            >  2:2834407/1‑148 (MQ=255)
       aTGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCg                                                                                                                                            >  2:1600673/1‑148 (MQ=255)
       aTGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCg                                                                                                                                            >  2:1800628/1‑148 (MQ=255)
         gCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGa                                                                                                                                          >  1:845824/1‑148 (MQ=255)
            ggTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACgg                                                                                                                                       >  1:3242292/1‑148 (MQ=255)
                 ttACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCAc                                                                                                                                  <  2:615534/148‑1 (MQ=255)
                    ccGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAg                                                                                                                               >  2:1186361/1‑148 (MQ=255)
                     cGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAgg                                                                                                                              <  1:598055/148‑1 (MQ=255)
                        cTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTcg                                                                                                                           >  1:3399802/1‑148 (MQ=255)
                        cTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTcg                                                                                                                           >  1:3399728/1‑148 (MQ=255)
                        cTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTcg                                                                                                                           >  1:3399710/1‑148 (MQ=255)
                         tCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTcgc                                                                                                                          <  2:802527/148‑1 (MQ=255)
                         tCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTcgc                                                                                                                          <  2:2720540/148‑1 (MQ=255)
                                tgtgGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCAt                                                                                                                   <  2:3063322/148‑1 (MQ=255)
                                  tgGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAg                                                                                                                 <  2:247092/148‑1 (MQ=255)
                                    gAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGAc                                                                                                               <  2:3376951/148‑1 (MQ=255)
                                     aaGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACg                                                                                                              >  1:2831496/1‑148 (MQ=255)
                                     aaGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACg                                                                                                              >  1:199598/1‑148 (MQ=255)
                                        ccAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGaa                                                                                                           >  2:1884902/1‑148 (MQ=255)
                                          aGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGaata                                                                                                         <  1:2810346/148‑1 (MQ=255)
                                           gTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGaataa                                                                                                        <  2:2700331/148‑1 (MQ=255)
                                            tACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGaataat                                                                                                       <  2:1833185/148‑1 (MQ=255)
                                               gATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCa                                                                                                    <  2:451795/148‑1 (MQ=255)
                                                  gTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGaa                                                                                                 >  2:1052763/1‑148 (MQ=255)
                                                  gTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGaa                                                                                                 >  2:1051274/1‑148 (MQ=255)
                                                  gTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGaa                                                                                                 <  2:231953/148‑1 (MQ=255)
                                                   tttAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAg                                                                                                >  2:117743/1‑148 (MQ=255)
                                                    ttAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAgg                                                                                               >  2:923520/1‑148 (MQ=255)
                                                       aCATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGc                                                                                            <  2:500718/148‑1 (MQ=255)
                                                       aCATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGc                                                                                            <  2:3290896/148‑1 (MQ=255)
                                                       aCATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGc                                                                                            >  2:303408/1‑148 (MQ=255)
                                                       aCATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGc                                                                                            <  2:2921107/148‑1 (MQ=255)
                                                       aCATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGc                                                                                            >  1:2077416/1‑148 (MQ=255)
                                                        cATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCa                                                                                           <  2:236315/148‑1 (MQ=255)
                                                        cATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCa                                                                                           >  2:2575726/1‑148 (MQ=255)
                                                            aaCGCCATCGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCgg                                                                                       <  1:46217/148‑1 (MQ=255)
                                                            aaCGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCgg                                                                                       <  2:3572502/148‑1 (MQ=255)
                                                             aCGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGt                                                                                      <  1:3016918/148‑1 (MQ=255)
                                                             aCGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGt                                                                                      <  2:3736795/148‑1 (MQ=255)
                                                              cGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTg                                                                                     <  1:2505603/148‑1 (MQ=255)
                                                               gCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGa                                                                                    <  2:3636445/148‑1 (MQ=255)
                                                                ccATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGAt                                                                                   <  2:1872712/148‑1 (MQ=255)
                                                                ccATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCACATAGACGGAATAATGCAGAAGGTGCAGCGGTGAt                                                                                   <  2:2074416/148‑1 (MQ=255)
                                                                 cATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGAtt                                                                                  <  1:3657246/148‑1 (MQ=255)
                                                                 cATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGAtt                                                                                  >  1:3754860/1‑148 (MQ=255)
                                                                  aTAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTa                                                                                 >  2:1999218/1‑148 (MQ=255)
                                                                    aGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGACTAATGCAGAAGGTGCAGCGGTGATTACa                                                                               >  2:3457598/1‑148 (MQ=255)
                                                                     ggCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAg                                                                              <  2:2835100/148‑1 (MQ=255)
                                                                      gCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGc                                                                             <  2:3180446/148‑1 (MQ=255)
                                                                      gCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGc                                                                             <  2:1817369/148‑1 (MQ=255)
                                                                      gCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGc                                                                             <  2:3704553/148‑1 (MQ=255)
                                                                       cGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGcc                                                                            <  1:2306352/148‑1 (MQ=255)
                                                                             tATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGa                                                                      >  2:2878190/1‑148 (MQ=255)
                                                                               tCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGaaa                                                                    >  2:2639671/1‑148 (MQ=255)
                                                                                ccTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAAt                                                                   >  2:3079025/1‑148 (MQ=255)
                                                                                      agaTCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAg                                                             <  2:1940282/148‑1 (MQ=255)
                                                                                      agaTCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAg                                                             <  2:1933674/148‑1 (MQ=255)
                                                                                       gaTCACCAGACTCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCCGGGCGCCGCGCATAGACGGAATAATGCAGAAGGTGCCGCGGTGATTACAGCCTTCAGAAATTTTCAGa                                                            >  1:3524200/1‑148 (MQ=255)
                                                                                       gaTCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGa                                                            >  2:2310680/1‑148 (MQ=255)
                                                                                       gaTCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGa                                                            >  2:3527914/1‑148 (MQ=255)
                                                                                            ccttaaTCTCTTTATCGCCCGCATCTACCAGACGTTTAGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGc                                                       <  1:3138776/144‑1 (MQ=255)
                                                                                                  tctcTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATg                                                 <  2:1757659/148‑1 (MQ=255)
                                                                                                  tctcTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATg                                                 <  2:1964257/148‑1 (MQ=255)
                                                                                                  tctcTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATg                                                 <  1:303408/148‑1 (MQ=255)
                                                                                                   ctctTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGa                                                <  1:2834407/148‑1 (MQ=255)
                                                                                                   ctctTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGa                                                <  1:127113/148‑1 (MQ=255)
                                                                                                   ctctTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGa                                                >  1:1040312/1‑148 (MQ=255)
                                                                                                    tctTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGAc                                               >  1:1158562/1‑148 (MQ=255)
                                                                                                          aCGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACgcggcg                                         >  2:1378422/1‑148 (MQ=255)
                                                                                                          aCGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACgcggcg                                         >  2:3283456/1‑148 (MQ=255)
                                                                                                          aCGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACgcggcg                                         >  2:781421/1‑148 (MQ=255)
                                                                                                          aCGCCCGCATCTACCAGACGTTTCGCGTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACgcggcg                                         >  1:2439479/1‑148 (MQ=255)
                                                                                                           cGCCCGCATCTACCAGACGTGTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGt                                        >  1:1829571/1‑148 (MQ=255)
                                                                                                               cGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGt                                    >  2:3120948/1‑148 (MQ=255)
                                                                                                               cGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGGCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGt                                    >  2:1587987/1‑148 (MQ=255)
                                                                                                                   tCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTca                                <  2:2974670/148‑1 (MQ=255)
                                                                                                                    cTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCTCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTcac                               >  1:3138059/1‑148 (MQ=255)
                                                                                                                      aCCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCacac                             >  1:2144115/1‑148 (MQ=255)
                                                                                                                       ccAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACAcc                            >  2:2506751/1‑148 (MQ=255)
                                                                                                                        cAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCt                           >  1:3488930/1‑148 (MQ=255)
                                                                                                                        cAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCt                           <  2:3236178/148‑1 (MQ=255)
                                                                                                                          gACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTg                         <  2:2267700/148‑1 (MQ=255)
                                                                                                                          gACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTg                         >  1:3186549/1‑148 (MQ=255)
                                                                                                                            cGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGtt                       <  1:1247069/148‑1 (MQ=255)
                                                                                                                            cGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGtt                       <  2:1194318/148‑1 (MQ=255)
                                                                                                                              tttCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCt                     <  2:3780840/148‑1 (MQ=255)
                                                                                                                              tttCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCt                     >  1:433552/1‑148 (MQ=255)
                                                                                                                                 cGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGc                  <  1:2724318/148‑1 (MQ=255)
                                                                                                                                   cTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCAc                <  2:1699031/148‑1 (MQ=255)
                                                                                                                                     tCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCa              <  2:3416529/148‑1 (MQ=255)
                                                                                                                                          tAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTc         >  2:1544485/1‑148 (MQ=255)
                                                                                                                                            acacTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAg       >  1:676474/1‑148 (MQ=255)
                                                                                                                                               cTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAGGaa    >  1:2954642/1‑148 (MQ=255)
                                                                                                                                                ttCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAGGAAt   <  1:1216751/148‑1 (MQ=255)
                                                                                                                                                ttCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAGGAAt   >  1:115553/1‑148 (MQ=255)
                                                                                                                                                 tCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAGGAATg  <  2:878640/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCTGACCATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAGGAATG  >  CP009273/872829‑873121

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: