Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 872,975 C→A 16.1% 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,9750CA16.1% 576.9 / 68.8 205E173* (GAA→TAA) rimOribosomal protein S12 methylthiotransferase; radical SAM superfamily
Reads supporting (aligned to +/- strand):  ref base C (98/74);  new base A (12/21);  total (110/95)
Fisher's exact test for biased strand distribution p-value = 3.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01

Only 100 of 205 total aligned reads displayed.

GGCTGACCATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAGGAAT  >  CP009273/872828‑873120
                                                                                                                                                   |                                                                                                                                                 
ggCTGACCATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCg                                                                                                                                                  >  1:2661489/1‑149 (MQ=255)
      ccATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATc                                                                                                                                            >  1:631435/1‑149 (MQ=255)
        aTGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCgg                                                                                                                                          <  2:1240943/149‑1 (MQ=255)
             ggTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGc                                                                                                                                     <  1:3742859/149‑1 (MQ=255)
             ggTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACACCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGc                                                                                                                                     <  1:2060691/149‑1 (MQ=255)
                    aCCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGATGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAg                                                                                                                              <  1:3927839/149‑1 (MQ=255)
                     ccGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAgg                                                                                                                             >  2:1835535/1‑149 (MQ=255)
                      cGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGt                                                                                                                            <  1:3002505/149‑1 (MQ=255)
                      cGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGt                                                                                                                            <  1:3002509/149‑1 (MQ=255)
                      cGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGt                                                                                                                            <  1:13019/149‑1 (MQ=255)
                        gCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCACTCGGACGGCTCACCAGGTcg                                                                                                                          >  2:144780/1‑149 (MQ=255)
                         cTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTcgc                                                                                                                         >  2:2072276/1‑149 (MQ=255)
                         cTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTcgc                                                                                                                         >  2:1674167/1‑149 (MQ=255)
                          tCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTcgcc                                                                                                                        >  2:2087109/1‑149 (MQ=255)
                          tCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTcgcc                                                                                                                        <  1:372978/149‑1 (MQ=255)
                               gTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCa                                                                                                                   <  2:2251902/149‑1 (MQ=255)
                                ttGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCAt                                                                                                                  >  1:2546572/1‑149 (MQ=255)
                                ttGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCAt                                                                                                                  >  2:1337938/1‑149 (MQ=255)
                                    ggAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGAc                                                                                                              >  2:370696/1‑149 (MQ=255)
                                       aGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGa                                                                                                           <  1:3402145/149‑1 (MQ=255)
                                        gCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGaa                                                                                                          <  1:2026160/149‑1 (MQ=255)
                                           aGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGaataa                                                                                                       <  2:3146003/149‑1 (MQ=255)
                                               cGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCa                                                                                                   <  2:1586219/149‑1 (MQ=255)
                                                gATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAg                                                                                                  >  1:1644180/1‑149 (MQ=255)
                                                gATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAg                                                                                                  >  2:3811371/1‑149 (MQ=255)
                                                 aTGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGa                                                                                                 >  1:1587839/1‑149 (MQ=255)
                                                  tGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGaa                                                                                                <  2:2869661/149‑1 (MQ=255)
                                                  tGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGaa                                                                                                <  2:2963582/149‑1 (MQ=255)
                                                  tGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGaa                                                                                                <  2:809595/149‑1 (MQ=255)
                                                  tGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGaa                                                                                                <  2:779244/149‑1 (MQ=255)
                                                   gTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAg                                                                                               >  2:478023/1‑149 (MQ=255)
                                                   gTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAg                                                                                               >  1:319179/1‑149 (MQ=255)
                                                   gTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAg                                                                                               <  2:2396080/149‑1 (MQ=255)
                                                    tttAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAgg                                                                                              >  2:3251325/1‑149 (MQ=255)
                                                    tttAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAgg                                                                                              >  2:1546262/1‑149 (MQ=255)
                                                       aaCATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGc                                                                                           <  1:2978623/149‑1 (MQ=255)
                                                       aaCATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGc                                                                                           <  2:1536760/149‑1 (MQ=255)
                                                          aTCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGc                                                                                        <  1:232637/149‑1 (MQ=255)
                                                              aCGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTg                                                                                    >  2:1661593/1‑149 (MQ=255)
                                                                gCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGAt                                                                                  <  2:2071889/149‑1 (MQ=255)
                                                                   aTAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTAc                                                                               >  1:1324313/1‑149 (MQ=255)
                                                                    tAGGCGGTAGTTTCCTTCTAGATCACCTGTATTTGTTTAACGCCAGCATCTACCAGACGTTTCGCTTCTCTTAACACTTCGCCTTTCGGTCGGCTCTCAAGGTCTCCTCGCATAGTCGGAATTTTGCAGAAGGTGCTTCGGTGAtttct                                                                              <  2:1912930/149‑4 (MQ=255)
                                                                    tAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCCTAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACa                                                                              >  1:698713/1‑149 (MQ=255)
                                                                    tAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCACATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACa                                                                              >  1:651275/1‑149 (MQ=255)
                                                                    tAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCACATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACa                                                                              >  1:651026/1‑149 (MQ=255)
                                                                     cggCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTCACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAg                                                                             <  2:1041948/148‑1 (MQ=255)
                                                                      ggCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATAACAGc                                                                            >  1:1101164/1‑149 (MQ=255)
                                                                            aGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAg                                                                      >  1:2207673/1‑149 (MQ=255)
                                                                              tATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGaa                                                                    >  1:3030812/1‑149 (MQ=255)
                                                                              tATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTAAAGCCTTCAGaa                                                                    >  1:3923116/1‑149 (MQ=255)
                                                                               aTCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGaaa                                                                   <  2:2834401/149‑1 (MQ=255)
                                                                                tCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAAt                                                                  >  1:1259079/1‑149 (MQ=255)
                                                                                tCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAAt                                                                  <  2:338850/149‑1 (MQ=255)
                                                                                 ccTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAAtt                                                                 >  1:934395/1‑149 (MQ=255)
                                                                                   tGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAAtttt                                                               <  2:2434176/149‑1 (MQ=255)
                                                                                     cGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCa                                                             <  2:631435/149‑1 (MQ=255)
                                                                                     cGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCa                                                             <  2:3507481/149‑1 (MQ=255)
                                                                                         aTCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTAAGATa                                                         <  2:4023569/149‑1 (MQ=255)
                                                                                            aCCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGc                                                      >  1:1598627/1‑149 (MQ=255)
                                                                                            aCCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGc                                                      >  1:3861330/1‑149 (MQ=255)
                                                                                             ccAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCa                                                     >  2:1592607/1‑149 (MQ=255)
                                                                                              cAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCAt                                                    <  1:1952805/149‑1 (MQ=255)
                                                                                               aGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATa                                                   <  2:568604/149‑1 (MQ=255)
                                                                                               aGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATa                                                   >  1:176353/1‑149 (MQ=255)
                                                                                               aGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATa                                                   <  2:1559748/149‑1 (MQ=255)
                                                                                                gAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGGATaa                                                  <  2:2546572/149‑1 (MQ=255)
                                                                                                gAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATaa                                                  >  1:2854691/1‑149 (MQ=255)
                                                                                                gAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATaa                                                  >  1:2876430/1‑149 (MQ=255)
                                                                                                 aaTCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAAt                                                 <  1:2815800/149‑1 (MQ=255)
                                                                                                   tctcTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGa                                               <  2:3696105/149‑1 (MQ=255)
                                                                                                     tctTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGAcg                                             >  1:3439261/1‑149 (MQ=255)
                                                                                                     tctTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGAcg                                             <  1:3626774/149‑1 (MQ=255)
                                                                                                        ttAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACgcgg                                          >  2:1003/1‑149 (MQ=255)
                                                                                                           aCGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGt                                       >  2:843122/1‑149 (MQ=255)
                                                                                                           aCGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGt                                       >  1:1539141/1‑149 (MQ=255)
                                                                                                               ccGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGt                                   <  1:357809/149‑1 (MQ=255)
                                                                                                                    tCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTcac                              >  1:3466375/1‑149 (MQ=255)
                                                                                                                    tCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTcac                              >  1:2371724/1‑149 (MQ=255)
                                                                                                                      tACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCacac                            >  2:281225/1‑149 (MQ=255)
                                                                                                                      tACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCacac                            >  1:3357982/1‑149 (MQ=255)
                                                                                                                      tACCAGACGTTTCGCTTCACTTAACACATCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCacac                            >  1:2103761/1‑149 (MQ=255)
                                                                                                                         cAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCtt                         <  1:3775517/149‑1 (MQ=255)
                                                                                                                            aCGTTTCGCTTCACTTAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGtt                      <  2:3637579/149‑1 (MQ=255)
                                                                                                                              gTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCt                    >  1:2655587/1‑149 (MQ=255)
                                                                                                                                 tCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGc                 <  1:141700/149‑1 (MQ=255)
                                                                                                                                 tCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGc                 <  1:1546204/149‑1 (MQ=255)
                                                                                                                                 tCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGc                 >  1:1888372/1‑149 (MQ=255)
                                                                                                                                  cGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCa                <  1:1013913/149‑1 (MQ=255)
                                                                                                                                     ttCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCa             >  2:1609789/1‑149 (MQ=255)
                                                                                                                                        aCTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGc          >  1:2391429/1‑149 (MQ=255)
                                                                                                                                         cTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCt         >  1:2650715/1‑149 (MQ=255)
                                                                                                                                          ttAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTc        >  1:3274542/1‑149 (MQ=255)
                                                                                                                                          ttAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTc        <  2:1987012/149‑1 (MQ=255)
                                                                                                                                           tAACACTTAGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCa       <  2:565225/149‑1 (MQ=255)
                                                                                                                                            aaCACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAg      <  1:1003/149‑1 (MQ=255)
                                                                                                                                             acacTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAgg     >  1:1109617/1‑149 (MQ=255)
                                                                                                                                              cacTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAGGa    >  1:741513/1‑149 (MQ=255)
                                                                                                                                               acTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAGGaa   >  2:2480510/1‑149 (MQ=255)
                                                                                                                                               acTTAGCCAATCGGACGGCTCACCGGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAGGaa   <  1:306148/149‑1 (MQ=255)
                                                                                                                                                cTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAGGAAt  >  2:2623644/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GGCTGACCATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCCATAGGCGGAAGTATCCTGCGAGATCACCAGAATCTCTTTAACGCCCGCATCTACCAGACGTTTCGCTTCACTTAACACTTCGCCAATCGGACGGCTCACCAGGTCGCCGCGCATAGACGGAATAATGCAGAAGGTGCAGCGGTGATTACAGCCTTCAGAAATTTTCAGATAGGCATAATGACGCGGCGTCAGTTTCACACCTTGTTCTGGCACCAGGCTCAGGAAT  >  CP009273/872828‑873120

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: