Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 872,975 C→A 15.0% 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,9750CA15.0% 604.3 / 56.5 200E173* (GAA→TAA) rimOribosomal protein S12 methylthiotransferase; radical SAM superfamily
Reads supporting (aligned to +/- strand):  ref base C (84/85);  new base A (21/10);  total (105/95)
Fisher's exact test for biased strand distribution p-value = 7.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 206 total aligned reads displayed.

GGCTGACCATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAGGAATGG  >  CP009273/872828‑873122
                                                                                                                                                   |                                                                                                                                                   
ggCTGACCATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTc                                                                                                                                                     >  2:993473/1‑148 (MQ=255)
ggCTGACCATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTc                                                                                                                                                     >  2:183629/1‑148 (MQ=255)
ggCTGACCATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTc                                                                                                                                                     >  2:3510525/1‑148 (MQ=255)
      ccATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAAt                                                                                                                                               >  1:3421835/1‑148 (MQ=255)
      ccATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAAt                                                                                                                                               <  2:198756/148‑1 (MQ=255)
         tGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCgg                                                                                                                                            <  2:2307921/148‑1 (MQ=255)
         tGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCgg                                                                                                                                            <  2:2331407/148‑1 (MQ=255)
           cTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGAc                                                                                                                                          <  2:2422995/148‑1 (MQ=255)
              gTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGc                                                                                                                                       >  1:1435986/1‑148 (MQ=255)
                  ttACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCAc                                                                                                                                   >  1:177748/1‑148 (MQ=255)
                    aCCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCa                                                                                                                                 >  2:1011208/1‑148 (MQ=255)
                    aCCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCCTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCa                                                                                                                                 >  2:3400303/1‑148 (MQ=255)
                     ccGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAg                                                                                                                                >  2:1197096/1‑148 (MQ=255)
                     ccGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAg                                                                                                                                >  1:469828/1‑148 (MQ=255)
                       ggCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGt                                                                                                                              <  2:854560/148‑1 (MQ=255)
                         cTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTcg                                                                                                                            >  1:2786696/1‑148 (MQ=255)
                            gCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTcgccg                                                                                                                         >  2:916661/1‑148 (MQ=255)
                             ccGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTcgccgc                                                                                                                        >  2:1923318/1‑148 (MQ=255)
                                ttGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGACGCGCa                                                                                                                     >  2:2025557/1‑148 (MQ=255)
                                   tgGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAg                                                                                                                  >  2:1620489/1‑148 (MQ=255)
                                         ccAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGaa                                                                                                            >  1:1634421/1‑148 (MQ=255)
                                           aGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGaata                                                                                                          <  2:2737129/148‑1 (MQ=255)
                                              aCGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATg                                                                                                       <  1:673260/148‑1 (MQ=255)
                                                 aTGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAg                                                                                                    >  1:673627/1‑148 (MQ=255)
                                                     ttAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAgg                                                                                                <  1:2972324/148‑1 (MQ=255)
                                                      tAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACGTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGt                                                                                               <  1:2158406/148‑1 (MQ=255)
                                                         cATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCa                                                                                            >  2:772941/1‑148 (MQ=255)
                                                         cATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCa                                                                                            <  2:2307136/148‑1 (MQ=255)
                                                            cAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCg                                                                                         <  2:109252/148‑1 (MQ=255)
                                                               cGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGACTAATGCAGAAGGTGCAGCGGTg                                                                                      >  2:1453223/1‑148 (MQ=255)
                                                               cGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTg                                                                                      >  2:1452952/1‑148 (MQ=255)
                                                               cGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTg                                                                                      >  2:3489791/1‑148 (MQ=255)
                                                               cGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTg                                                                                      <  1:1044438/148‑1 (MQ=255)
                                                                   aTAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGACTAATGCAGAAGGTGCAGCGGTGATTa                                                                                  >  2:1863623/1‑148 (MQ=255)
                                                                   aTAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTa                                                                                  <  1:145791/148‑1 (MQ=255)
                                                                   aTAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTa                                                                                  >  2:61503/1‑148 (MQ=255)
                                                                    tAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTCACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTAc                                                                                 <  2:2769085/148‑1 (MQ=255)
                                                                    tAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTAc                                                                                 <  2:3057354/148‑1 (MQ=255)
                                                                     aGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACa                                                                                >  2:694315/1‑148 (MQ=255)
                                                                     aGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACa                                                                                <  1:834060/148‑1 (MQ=255)
                                                                       gCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGc                                                                              >  2:3364399/1‑148 (MQ=255)
                                                                        cGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGcc                                                                             <  1:1962960/148‑1 (MQ=255)
                                                                         ggAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCt                                                                            <  1:2734619/148‑1 (MQ=255)
                                                                             gTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAg                                                                        >  2:547034/1‑148 (MQ=255)
                                                                                tCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGaaa                                                                     <  1:2487002/148‑1 (MQ=255)
                                                                                tCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGaaa                                                                     <  1:2489133/148‑1 (MQ=255)
                                                                                   tGCGAGATCACCAGACTCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAAttt                                                                  <  1:1633289/148‑1 (MQ=255)
                                                                                   tGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAAttt                                                                  >  1:1861756/1‑148 (MQ=255)
                                                                                    gCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAAtttt                                                                 <  2:469828/148‑1 (MQ=255)
                                                                                       agaTCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAg                                                              <  1:2121482/148‑1 (MQ=255)
                                                                                         aTCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGAt                                                            <  2:2949459/148‑1 (MQ=255)
                                                                                         aTCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGAt                                                            >  2:2783219/1‑148 (MQ=255)
                                                                                         aTCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTAAGAAATTTTCAGAt                                                            <  1:2673113/148‑1 (MQ=255)
                                                                                             ccAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCCGCGGTGATTACAGCCTTCAGAAAATTTCAGATAGGc                                                        >  2:3097996/1‑148 (MQ=255)
                                                                                             ccAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGc                                                        >  2:1773165/1‑148 (MQ=255)
                                                                                              cAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCa                                                       <  2:1427572/148‑1 (MQ=255)
                                                                                                gAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATa                                                     <  2:1004532/148‑1 (MQ=255)
                                                                                                 aaTCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATaa                                                    >  2:1052656/1‑148 (MQ=255)
                                                                                                   tctcTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATg                                                  <  2:176210/148‑1 (MQ=255)
                                                                                                   tctcTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATg                                                  <  2:728191/148‑1 (MQ=255)
                                                                                                   tctcTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATg                                                  <  2:1120842/148‑1 (MQ=255)
                                                                                                   tctcTTTAACGCCCGCATCTACCAGACGTTTCGCTACACTTAACACTTAGCCAATCTTAAGGCTCACCACGTCGACGCGCATATACGGAATAAACCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTATCAGATAGGAATAATg                                                  >  1:2354054/1‑148 (MQ=255)
                                                                                                    ctctTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGa                                                 >  2:3293196/1‑148 (MQ=255)
                                                                                                      ctTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGAcg                                               <  2:2916332/148‑1 (MQ=255)
                                                                                                        ttAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGAcgcg                                             >  1:194172/1‑148 (MQ=255)
                                                                                                         tAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCCATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACgcgg                                            >  1:47577/1‑148 (MQ=255)
                                                                                                           aCGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACgcggcg                                          <  1:3066609/148‑1 (MQ=255)
                                                                                                            cGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGt                                         <  1:993473/148‑1 (MQ=255)
                                                                                                            cGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGt                                         <  1:1952450/148‑1 (MQ=255)
                                                                                                            cGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGt                                         >  2:475564/1‑148 (MQ=255)
                                                                                                             gCCCGCATCTACCAGACGTTTAGCTTCACTTAAAACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGAAGCGGCGTc                                        >  1:3055535/1‑148 (MQ=255)
                                                                                                              cccGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCa                                       >  1:3163431/1‑148 (MQ=255)
                                                                                                                 gCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGtt                                    <  1:1011208/148‑1 (MQ=255)
                                                                                                                  cATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGttt                                   <  2:956619/148‑1 (MQ=255)
                                                                                                                   aTCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTc                                  <  1:2810071/148‑1 (MQ=255)
                                                                                                                   aTCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTc                                  >  2:1730055/1‑148 (MQ=255)
                                                                                                                    tCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTca                                 <  2:150447/148‑1 (MQ=255)
                                                                                                                     cTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTcac                                >  2:2230586/1‑148 (MQ=255)
                                                                                                                     cTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATTCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTcac                                >  1:1336259/1‑148 (MQ=255)
                                                                                                                       aCCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCacac                              >  2:2097287/1‑148 (MQ=255)
                                                                                                                        ccAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACAcc                             >  2:837495/1‑148 (MQ=255)
                                                                                                                           gACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTg                          >  1:3307669/1‑148 (MQ=255)
                                                                                                                             cGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAGTAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGtt                        <  1:874389/148‑1 (MQ=255)
                                                                                                                             cGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGtt                        <  1:109577/148‑1 (MQ=255)
                                                                                                                             cGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGtt                        <  1:2765563/148‑1 (MQ=255)
                                                                                                                                ttCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTg                     >  1:508886/1‑148 (MQ=255)
                                                                                                                                 tCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTgg                    <  2:1634421/148‑1 (MQ=255)
                                                                                                                                 tCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTgg                    >  2:2023574/1‑148 (MQ=255)
                                                                                                                                    cTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCAc                 <  2:1569662/148‑1 (MQ=255)
                                                                                                                                    cTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCAc                 >  2:2726203/1‑148 (MQ=255)
                                                                                                                                    cTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCAc                 >  1:2291220/1‑148 (MQ=255)
                                                                                                                                     ttCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCAcc                >  2:248648/1‑148 (MQ=255)
                                                                                                                                        aCTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAgg             <  2:1037512/148‑1 (MQ=255)
                                                                                                                                          ttAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCt           <  1:1197096/148‑1 (MQ=255)
                                                                                                                                            aaCACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCa         >  1:2887023/1‑148 (MQ=255)
                                                                                                                                             acacTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAg        >  2:1975119/1‑148 (MQ=255)
                                                                                                                                               acTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAGGa      <  1:1900881/148‑1 (MQ=255)
                                                                                                                                                  tCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAGGAATg   <  2:2384721/148‑1 (MQ=255)
                                                                                                                                                   cGCCACTCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAGGAATgg  <  1:572833/148‑1 (MQ=255)
                                                                                                                                                   cGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAGGAATgg  <  1:2495754/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGCTGACCATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAGGAATGG  >  CP009273/872828‑873122

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: