Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 872,975 C→A 6.4% 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,9750CA 6.4% 590.1 / 12.1 163E173* (GAA→TAA) rimOribosomal protein S12 methylthiotransferase; radical SAM superfamily
Reads supporting (aligned to +/- strand):  ref base C (74/78);  new base A (3/8);  total (77/86)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 163 total aligned reads displayed.

GACCATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAGGAATG  >  CP009273/872832‑873121
                                                                                                                                               |                                                                                                                                                  
   cATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATc                                                                                                                                             >  1:1318418/1‑148 (MQ=255)
   cATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATc                                                                                                                                             >  1:490947/1‑148 (MQ=255)
    aTGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCg                                                                                                                                            >  2:2606958/1‑148 (MQ=255)
     tGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCTTAGGCGGAAGTATCCTGCGATATAACCTGAATCTCTTTAACGCCCGCATCTACCTGACGTTTCGCTTCACTTAACACTTCGCCAATCgg                                                                                                                                           <  2:676134/148‑1 (MQ=255)
           tttttACCGGCTCGCCGTTGTGGAAGCCAGTACGATTTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCt                                                                                                                                     <  1:104020/148‑1 (MQ=255)
           tttttACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCt                                                                                                                                     <  1:104048/148‑1 (MQ=255)
             tttACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCa                                                                                                                                   >  1:974339/1‑148 (MQ=255)
             tttACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAAAACTTCGCCAATCGGACGGCTCa                                                                                                                                   >  1:2669613/1‑148 (MQ=255)
                aCCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCa                                                                                                                                <  2:563666/148‑1 (MQ=255)
                 ccGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAg                                                                                                                               <  2:815744/148‑1 (MQ=255)
                   ggCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGt                                                                                                                             <  1:645773/148‑1 (MQ=255)
                       cGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTcgcc                                                                                                                         >  1:2700046/1‑148 (MQ=255)
                         ccGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTcgccgc                                                                                                                       >  1:1456361/1‑148 (MQ=255)
                          cGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCcgcg                                                                                                                      >  2:2644213/1‑148 (MQ=255)
                            ttGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCa                                                                                                                    >  2:13675/1‑148 (MQ=255)
                             tgtgGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCAt                                                                                                                   <  1:2976788/148‑1 (MQ=255)
                             tgtgGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCAt                                                                                                                   <  1:2631539/148‑1 (MQ=255)
                             tgtgGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCAt                                                                                                                   <  1:2344732/148‑1 (MQ=255)
                               tgGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAg                                                                                                                 >  2:943828/1‑148 (MQ=255)
                               tgGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAg                                                                                                                 >  2:1991048/1‑148 (MQ=255)
                               tgGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAg                                                                                                                 >  2:664/1‑148 (MQ=255)
                                 gAAGCCAGTACTATGTTTAACATAAACGCCATAGGCGTATGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGAc                                                                                                               <  1:570550/148‑1 (MQ=255)
                                 gAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGAc                                                                                                               <  2:1074393/148‑1 (MQ=255)
                                 gAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGAc                                                                                                               >  2:2976218/1‑148 (MQ=255)
                                    gCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGa                                                                                                            <  1:1110339/148‑1 (MQ=255)
                                     ccAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGaa                                                                                                           <  1:2071277/148‑1 (MQ=255)
                                        gTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGaataa                                                                                                        >  2:888862/1‑148 (MQ=255)
                                        gTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGaataa                                                                                                        >  1:3304691/1‑148 (MQ=255)
                                        gTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGaataa                                                                                                        >  2:2577601/1‑148 (MQ=255)
                                          aCGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATg                                                                                                      <  1:1566207/148‑1 (MQ=255)
                                           cGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGc                                                                                                     >  1:2819729/1‑148 (MQ=255)
                                             aTGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAg                                                                                                   >  1:1999214/1‑148 (MQ=255)
                                              tGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGa                                                                                                  >  2:188749/1‑148 (MQ=255)
                                              tGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGa                                                                                                  >  2:3256779/1‑148 (MQ=255)
                                                  tAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGt                                                                                              <  1:165081/148‑1 (MQ=255)
                                                   aaCATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTg                                                                                             <  1:243065/148‑1 (MQ=255)
                                                       tCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGc                                                                                         >  1:379250/1‑148 (MQ=255)
                                                       tCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGc                                                                                         >  2:3199525/1‑148 (MQ=255)
                                                         aaCGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCgg                                                                                       <  2:127061/148‑1 (MQ=255)
                                                          aCGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGt                                                                                      <  1:961828/148‑1 (MQ=255)
                                                           cGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCtgtg                                                                                     <  2:2136953/148‑1 (MQ=255)
                                                             ccATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGAt                                                                                   <  2:1690370/148‑1 (MQ=255)
                                                               aTAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTa                                                                                 <  1:820109/148‑1 (MQ=255)
                                                                tAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGGGATTAc                                                                                >  1:1406390/1‑148 (MQ=255)
                                                                 aGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACa                                                                               >  2:1560248/1‑148 (MQ=255)
                                                                  ggCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAg                                                                              >  1:493630/1‑148 (MQ=255)
                                                                   gggtaattATCCTGCGAGATCACCAGAATCTATTTCACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGc                                                                             <  1:2321401/141‑1 (MQ=255)
                                                                    cGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGcc                                                                            <  1:2532926/148‑1 (MQ=255)
                                                                    cGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGcc                                                                            <  1:2992438/148‑1 (MQ=255)
                                                                     ggAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCt                                                                           <  1:3195947/148‑1 (MQ=255)
                                                                     ggAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACCTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCt                                                                           <  1:3195941/148‑1 (MQ=255)
                                                                      gAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCtt                                                                          <  1:719025/148‑1 (MQ=255)
                                                                      gAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCtt                                                                          >  2:603428/1‑148 (MQ=255)
                                                                      gAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCtt                                                                          <  1:1105059/148‑1 (MQ=255)
                                                                       aaGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTc                                                                         >  1:2286766/1‑148 (MQ=255)
                                                                          tATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGa                                                                      >  1:2012349/1‑148 (MQ=255)
                                                                             ccTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAAt                                                                   <  1:1089720/148‑1 (MQ=255)
                                                                             ccTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAAt                                                                   <  1:783934/148‑1 (MQ=255)
                                                                                gCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAAtttt                                                                <  1:1914759/148‑1 (MQ=255)
                                                                                gCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCGTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCAAAGACGGAAAAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAAtttt                                                                >  1:2946452/1‑148 (MQ=255)
                                                                                   agaTCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAg                                                             <  2:2700087/148‑1 (MQ=255)
                                                                                     aTCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGAt                                                           <  2:3134873/148‑1 (MQ=255)
                                                                                     aTCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGAt                                                           <  2:2961480/148‑1 (MQ=255)
                                                                                     aTCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTAAGAt                                                           >  1:1132289/1‑148 (MQ=255)
                                                                                      tCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATa                                                          <  1:1393156/148‑1 (MQ=255)
                                                                                      tCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATa                                                          <  1:1982300/148‑1 (MQ=255)
                                                                                       cACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAg                                                         >  2:304301/1‑148 (MQ=255)
                                                                                          cAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCa                                                      <  2:1325106/148‑1 (MQ=255)
                                                                                          cAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCa                                                      <  2:437903/148‑1 (MQ=255)
                                                                                          cAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCa                                                      <  1:1142479/148‑1 (MQ=255)
                                                                                          cAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCa                                                      >  2:2858700/1‑148 (MQ=255)
                                                                                            gactctcTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATa                                                    <  2:1965481/145‑1 (MQ=255)
                                                                                            gAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATa                                                    <  1:3080606/148‑1 (MQ=255)
                                                                                            gAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATa                                                    <  2:340991/148‑1 (MQ=255)
                                                                                              aTCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAAt                                                  >  2:3040997/1‑148 (MQ=255)
                                                                                              aTCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAAt                                                  <  2:3304691/148‑1 (MQ=255)
                                                                                               tctcTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATg                                                 <  2:338950/148‑1 (MQ=255)
                                                                                                ctctTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGGCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGa                                                <  2:2028830/148‑1 (MQ=255)
                                                                                                ctctTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGa                                                <  2:2028782/148‑1 (MQ=255)
                                                                                                ctctTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGa                                                >  2:391968/1‑148 (MQ=255)
                                                                                                  ctTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGAcg                                              >  2:773952/1‑148 (MQ=255)
                                                                                                  ctTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGAcg                                              <  2:867994/148‑1 (MQ=255)
                                                                                                                tCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTca                                >  1:697570/1‑148 (MQ=255)
                                                                                                                 cTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTcac                               >  1:2255429/1‑148 (MQ=255)
                                                                                                                   aCCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCacac                             <  2:2706113/148‑1 (MQ=255)
                                                                                                                   aCCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCacac                             >  1:1200397/1‑148 (MQ=255)
                                                                                                                   aCCAGACGTTTCGCTTCACTTAACACATCGCCCAACGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTAAAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCacac                             >  1:893225/1‑148 (MQ=255)
                                                                                                                           tttCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCt                     >  1:431571/1‑148 (MQ=255)
                                                                                                                                cTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCAc                >  1:790172/1‑148 (MQ=255)
                                                                                                                                cTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCAc                >  2:2107243/1‑148 (MQ=255)
                                                                                                                                  tCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCa              <  1:1688964/148‑1 (MQ=255)
                                                                                                                                  tCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCa              <  2:840786/148‑1 (MQ=255)
                                                                                                                                    aCTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAgg            >  2:867790/1‑148 (MQ=255)
                                                                                                                                      ttAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCt          >  1:1992913/1‑148 (MQ=255)
                                                                                                                                        aaCACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCa        >  2:2897135/1‑148 (MQ=255)
                                                                                                                                          cacTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAgg      <  2:773601/148‑1 (MQ=255)
                                                                                                                                           acTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAGGa     <  2:459618/148‑1 (MQ=255)
                                                                                                                                            cTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAGGaa    <  1:664/148‑1 (MQ=255)
                                                                                                                                              tCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAGGAATg  >  2:838405/1‑148 (MQ=255)
                                                                                                                                              tCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAGGAATg  >  2:2495350/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GACCATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAGGAATG  >  CP009273/872832‑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: