Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 843,246 T→A 8.6% N305K (AAT→AAA)  rlmF → 23S rRNA m(6)A1618 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913843,2460TA 8.6% 2389.0 / 35.3 782N305K (AAT→AAA) rlmF23S rRNA m(6)A1618 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (277/413);  new base A (88/2);  total (365/417)
Fisher's exact test for biased strand distribution p-value = 2.44e-29
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 790 total aligned reads displayed.

CGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTAC  >  NC_000913/843097‑843395
                                                                                                                                                     |                                                                                                                                                     
    gtgAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGcc                                                                                                                                                   <  2:13422251/150‑1 (MQ=255)
       aaaaCTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGCGAAGGTGGTTAAAAAAGAGATGACCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGc                                                                                                                                                >  2:11726198/1‑150 (MQ=255)
         aaCTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGt                                                                                                                                              <  2:7841856/150‑1 (MQ=255)
            ttACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTATTGCCTGGACCTTTATGAACGACGAGCCGCGCCGCCGTTTTGTCAAACGCCAGCGGTaa                                                                                                                                           >  2:11489197/1‑150 (MQ=255)
                  ccgTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGtt                                                                                                                                     <  2:12104212/150‑1 (MQ=255)
                   cgTTGTATCGTGCCCTGTCGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTTTCAATCGCCAGCGTTAAAGCGTTg                                                                                                                                    <  1:3181714/150‑1 (MQ=255)
                   cgTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTg                                                                                                                                    <  1:14730727/150‑1 (MQ=255)
                     ttGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTATTGCCGGGACCATTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCACCGTTAAAGCGTTGGt                                                                                                                                  >  2:15264737/1‑150 (MQ=255)
                       gTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTTGCTGGACCATTATGAAAGACGAGCAGCGCCGCCGGTTTTTCAATAGCCAGCGGTTAAGCGTTgggtg                                                                                                                                >  2:2879583/1‑147 (MQ=255)
                        tATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTTCCTGGACCATTATTTACGACGAGCAGCCCCCCCGATTTGTAAAAAGCCACCGGTAAATCGGTGGttac                                                                                                                               >  1:1289179/1‑148 (MQ=255)
                         aTCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTTCCTGGACCGTTATGAACGACGAGCAGCGCCGCCGTTTTTTCAATCGCCAGCGGTAAAGCGTTGGTTgcc                                                                                                                              >  2:10113456/1‑150 (MQ=255)
                             tGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCcgg                                                                                                                          <  2:14472354/150‑1 (MQ=255)
                               cccTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTTTTTCCTGGACCATTATGTACCACGAGCAGCGCCGCCGTTTTGTGAAACGCCAGCCGTAAAGCGGTGGTTGCCCccgcg                                                                                                                        >  1:13720920/1‑150 (MQ=255)
                               cccTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGGTATTTCCGGGACCTATATGAACGACCAGCCGCGCCCCCGGTTTTTCAAAACCCAGCGGTAAACCGGTGGTTGCCCccgcg                                                                                                                        >  2:4399675/1‑150 (MQ=255)
                                ccTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCcggcgg                                                                                                                       <  1:4040899/150‑1 (MQ=255)
                                ccTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCGGGACCTTTATGAACCACGAGCCGCGCCGCCGTTTTGTCAAAAGCCACCGTTAAAGCGTTTGTTACCCccgcgg                                                                                                                       >  2:7593795/1‑150 (MQ=255)
                                 cTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTTCCGGGACCGTTATGAACGACGAGCAGCGCCGCCGGTTTGTCAATAGCCAGCGTTAAAGCGGTGGTTCCCCcggcggg                                                                                                                      >  2:8254465/1‑149 (MQ=255)
                                    acggacgTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTGTATTGCCGGGACCGTTATTAACCACCAGCAGCCCCGCCGGTTTGTCAATAGCCCGCGTTAAAGCGATGGTTACCCCCGCGGGAAc                                                                                                                   >  2:12299985/1‑150 (MQ=255)
                                    acggacgTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTTTTTGCTGGGCCGTTTTGAAAGACGAGCCGCGCCCCCGGGTTGTAAAAAGCCAACGCTAAATCGTTGGTTGCACccgcggggac                                                                                                                   >  1:5917611/1‑146 (MQ=255)
                                     cggacgTGGGCGCGGTGAATGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACg                                                                                                                  <  2:6305608/150‑1 (MQ=255)
                                     cggacgGGGGCGCGGTGAAGGTGGTTAACAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACg                                                                                                                  <  2:14978877/150‑1 (MQ=255)
                                      ggacgTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACgg                                                                                                                 <  1:3777670/150‑1 (MQ=255)
                                        acgTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCGGGACCATTAAGAAACACGAGCACCGCCCCCCGTTTTTAAAAAGCCACCCGTAAATCGTTGGTTACCCCCGCGGGAAAgggt                                                                                                               >  2:154660/1‑148 (MQ=255)
                                            gggCGCGGTGAAGGGGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGCACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTACAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGc                                                                                                           <  2:5478607/150‑1 (MQ=255)
                                               cgcgGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCgg                                                                                                        <  2:8798989/150‑1 (MQ=255)
                                                gcgGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCggg                                                                                                       <  1:15052129/150‑1 (MQ=255)
                                                gcgGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGGTCCCCCGGCGGTAACGGTTCAGCCggg                                                                                                       <  2:10682093/150‑1 (MQ=255)
                                                  ggTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTATTGCCAGGACCTTTATGAACCACCAGCAGCGCCGCCCTTTTTTCAATCGCCAGCCTTAAAGCGTTGGTTCCCCCCGCGGTAACCGTTCAGCCCGGGt                                                                                                     >  2:5131633/1‑150 (MQ=255)
                                                        ggtggtTAAAAAAGAGATAGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGcc                                                                                               <  2:11011078/150‑1 (MQ=255)
                                                         gtggtTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGccc                                                                                              <  2:14697633/150‑1 (MQ=255)
                                                               aaaaaaGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAgca                                                                                        <  2:5950067/150‑1 (MQ=255)
                                                                 aaaaGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAgcagg                                                                                      <  2:8091052/150‑1 (MQ=255)
                                                                  aaaGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTTCCGGGACCTTTATGAAAGACGAGCAGCGCCGCCGTTTTGTCAAACGCCAGCGTTAAAGCGTGGGTTACCCCGGCGGTAAAGGTTCAGCCGGGGTCGCGCCCGGAgcaggc                                                                                     >  2:3653267/1‑150 (MQ=255)
                                                                  aaaGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTATTGCCGGGAACAATAATAAAAAAGATCAACGCCGCCGTTTTTTCAATAGACAGCGTTAAATAGTTTGTTGCCCCCGCCGTAAAGGTTCATCCGGGGTCGCGCCCGGTgccggc                                                                                     >  2:1989239/1‑150 (MQ=255)
                                                                   aaGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAgcaggca                                                                                    <  2:2620140/150‑1 (MQ=255)
                                                                       gaTGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCAc                                                                                <  2:14058259/150‑1 (MQ=255)
                                                                        aTGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCAcc                                                                               <  2:14371819/150‑1 (MQ=255)
                                                                            cccAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGAt                                                                           <  2:7163611/150‑1 (MQ=255)
                                                                             ccAGGGGCAAAAGCAGAGTCGCTTTATTTCCTGGACCTTTATGAACGACGAGCCGCGCCGCCGTTTTGTCAATCGCCAGCGGTAAAGCGTTTGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCCGAGCAGGCAGCACCAGATa                                                                          >  2:10410750/1‑150 (MQ=255)
                                                                             ccAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATa                                                                          <  1:7477225/150‑1 (MQ=255)
                                                                              cAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATaa                                                                         <  1:212361/150‑1 (MQ=255)
                                                                              cAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATaa                                                                         <  1:1342174/150‑1 (MQ=255)
                                                                              cAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATaa                                                                         <  1:9217708/150‑1 (MQ=255)
                                                                                  ggCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCt                                                                     <  2:6610275/150‑1 (MQ=255)
                                                                                  ggCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCt                                                                     <  1:12499350/150‑1 (MQ=255)
                                                                                    cAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGc                                                                   <  1:4326264/150‑1 (MQ=255)
                                                                                             agTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGgcgcg                                                          <  2:952595/150‑1 (MQ=255)
                                                                                               tCGCTTTATTGCCTGGACCTTTATTAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCCGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGcgcgcg                                                        >  1:10968836/1‑150 (MQ=255)
                                                                                                 gCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGaa                                                      >  1:14190649/1‑150 (MQ=255)
                                                                                                  cTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAAc                                                     <  2:13128495/150‑1 (MQ=255)
                                                                                                   tttATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAAcg                                                    >  1:1615231/1‑150 (MQ=255)
                                                                                                      aTTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAAcgccc                                                 >  1:6293485/1‑150 (MQ=255)
                                                                                                       ttGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAAcgcccg                                                >  2:6148819/1‑150 (MQ=255)
                                                                                                       ttGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAAcgcccg                                                >  2:13943281/1‑150 (MQ=255)
                                                                                                       ttGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAAcgcccg                                                >  1:4753504/1‑150 (MQ=255)
                                                                                                        tGCCTGGACCTTTATGAACGACGAGCCGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGGGTTGGTTCCCCCGGCGGTACCGGTTCAGCGGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAAcgcccgc                                               <  1:9630940/150‑1 (MQ=255)
                                                                                                          ccTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCa                                             >  1:2661262/1‑150 (MQ=255)
                                                                                                          ccTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCa                                             >  1:6672197/1‑150 (MQ=255)
                                                                                                           cTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAg                                            >  1:5376718/1‑150 (MQ=255)
                                                                                                            tGGACCTTTAAGAACTACGAGCATCTCAGCCGTTTATTCAATCTCAAGCGTTAAAGCGTTTGTTCCACCCGCGGTAACGGATCAGCCGGGCTCGCGCCCGGAGCAGGCAGCACCTGAGAAGTCTGCACTGGCGCGCGAACGCCCGCCCgg                                           >  2:865025/1‑150 (MQ=255)
                                                                                                             ggACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAgcaggca                                                                                    <  1:13933015/108‑1 (MQ=255)
                                                                                                                 cTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGaa                                      <  2:9172209/150‑1 (MQ=255)
                                                                                                                       gAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAATGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGTTCGCGCCCGGAGCTGGCAGCACCTGATATGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGttt                                <  1:6982483/150‑1 (MQ=255)
                                                                                                                       gAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGGCGAAATGttt                                >  1:13218557/1‑150 (MQ=255)
                                                                                                                       gAACGACGAGCAGCGCCGCCGTTTTGTCAATAGCCAGCGTTAAATCGTTGGTTACCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGttt                                >  1:774533/1‑150 (MQ=255)
                                                                                                                        cccgacgaGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGTGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGtttt                               <  1:5163493/148‑1 (MQ=255)
                                                                                                                         acgacgAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGttttt                              >  1:1918140/1‑150 (MQ=255)
                                                                                                                          cgacgaGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGCAATGtttttt                             >  1:15221884/1‑150 (MQ=255)
                                                                                                                           gacgaGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTc                            >  2:97035/1‑150 (MQ=255)
                                                                                                                           gacgaGCAGCGCCGCCGGTTTGTCAATCGCCAGCGTTAAAGCGGTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATACGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTc                            >  2:13846440/1‑150 (MQ=255)
                                                                                                                              gaGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCAcc                         <  2:12881895/150‑1 (MQ=255)
                                                                                                                               agcagcGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCt                        <  1:10467064/150‑1 (MQ=255)
                                                                                                                               agcagcGCCGCCGTTTTGTCAATCGCCAGCGGTAAAGCGTTGGTTCCCCCCGCGGTAACGGGTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCGGATAAGTATGCACTGGCGCGCGAACGCCCCCCACGTCGAAATTTTTTTTCACCt                        >  1:6098341/1‑150 (MQ=255)
                                                                                                                                gcagcGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTg                       >  2:5789258/1‑150 (MQ=255)
                                                                                                                                 cagcGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTgg                      <  1:10096304/150‑1 (MQ=255)
                                                                                                                                 cagcGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTgg                      <  1:3359262/150‑1 (MQ=255)
                                                                                                                                 cagcGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTgg                      >  2:11286627/1‑150 (MQ=255)
                                                                                                                                  agcGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGc                     <  1:8932540/150‑1 (MQ=255)
                                                                                                                                    cgccgcCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTg                   >  2:5691965/1‑150 (MQ=255)
                                                                                                                                    cgccgcCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGAGAAGTCTGCACTGGCGCGCGACCGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTg                   >  2:3429438/1‑150 (MQ=255)
                                                                                                                                       cgccgTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCg                >  2:7430494/1‑150 (MQ=255)
                                                                                                                                        gccgTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGa               <  1:14588618/150‑1 (MQ=255)
                                                                                                                                        gccgTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGa               >  1:9968134/1‑150 (MQ=255)
                                                                                                                                        gccgTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGa               >  1:8082560/1‑150 (MQ=255)
                                                                                                                                         ccgTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAg              <  2:12650168/150‑1 (MQ=255)
                                                                                                                                         ccgTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAg              <  1:15425421/150‑1 (MQ=255)
                                                                                                                                          cgTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAgg             <  2:6823988/150‑1 (MQ=255)
                                                                                                                                           gTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAggg            <  1:11964297/150‑1 (MQ=255)
                                                                                                                                                gTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGtt                                 <  2:8422336/124‑1 (MQ=255)
                                                                                                                                                gTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGtt                                 >  1:8422336/1‑124 (MQ=255)
                                                                                                                                                gTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCTCCTGGCTGTCGAGGGCAAAg       >  2:3068494/1‑150 (MQ=255)
                                                                                                                                                 tCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGc      <  2:4633947/150‑1 (MQ=255)
                                                                                                                                                   aaTCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCAc                          >  1:11643982/1‑128 (MQ=255)
                                                                                                                                                   aaTCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTc                 <  1:10202919/137‑1 (MQ=255)
                                                                                                                                                   aaTCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGt    >  2:2274308/1‑150 (MQ=255)
                                                                                                                                                   aaTCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGt    >  2:3937916/1‑150 (MQ=255)
                                                                                                                                                    aTCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTa   <  2:4940429/150‑1 (MQ=255)
                                                                                                                                                    aTCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTa   <  1:11264654/150‑1 (MQ=255)
                                                                                                                                                    aTCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTa   <  2:1834915/150‑1 (MQ=255)
                                                                                                                                                     tCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTAc  >  1:8305431/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTAC  >  NC_000913/843097‑843395

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑