Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 842,815 C→A 19.1% Q162K (CAA→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_000913842,8150CA19.1% 133.8 / 2.8 53Q162K (CAA→AAA) rlmF23S rRNA m(6)A1618 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (22/20);  new base A (2/9);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 8.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGGCGGATTATATTCATCACCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGTGTGAAGGCGGTGAAGTCACCTTTATCAAAAAGATGATTGAAGAGAGCA  >  NC_000913/842578‑843052
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                
cgggcGGATTTATTTCATCCCCTTGCCGATTTCCTGGCAGAAGCAAGCGGAAAAATTCCAGCGAATGCCATCATTATGGATACTGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGGGTTTTACCGGTAGTGAAACCAGCCGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGc                                                                                                                                                                                                                                                  <  2:804806/238‑1 (MQ=255)
  ggcGGATTATATTCATCACCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCaaa                                                                                                                                                                                                                                               >  1:279598/1‑239 (MQ=255)
  ggcGGATTATATTCATCACCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCCCCATTGGTCTGCGTCGGCaaa                                                                                                                                                                                                                                               >  1:166844/1‑239 (MQ=255)
              tcatccCCTTGCCTATTTACTGGCAGAAGCAAGGGGAACAATTCCCGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTgggg                                                                                                                                                                                                                                   <  2:269627/239‑1 (MQ=255)
              tcaacaCCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTgggg                                                                                                                                                                                                                                   <  1:264163/239‑1 (MQ=255)
                 tcaCCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGACTGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATACGCTGATTGGCGTACATGACTATGGCTGGAGTTTTACCGGTAGGGAAGACAGCGGCCCGGGGTTAAGCAGATGGCAGGCGATTTTCGGATATAAACCAGGGCATAACAGGCCCATTCGGCTGCGGCGGAACAAAGAGAGGGGGGCGa                                                                                                                                                                                                                               >  2:801309/1‑239 (MQ=255)
                   aCCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGGGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGGCCCATTCGTCTGCGTCGGCAAAAAAAGAGTGGGGCGAtt                                                                                                                                                                                                                             >  1:678892/1‑239 (MQ=255)
                       tGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGGGTTAAGCAGTGCGCAGGCGCTTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGAAAAAAGAGAGGGGGGCGAGTTTTa                                                                                                                                                                                                                         >  2:582767/1‑239 (MQ=255)
                        tgcgttttCCTGGCAGAAGGAAGGGGAACAATTCCAGCGAATGCCAGCATTCTGGAGAGCGGTGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCGGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTaa                                                                                                                                                                                                                        <  2:861553/234‑1 (MQ=255)
                        gCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTaa                                                                                                                                                                                                                        <  1:582767/239‑1 (MQ=255)
                                         agcaagcGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGTTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATaaaa                                                                                                                                                                                                       <  1:632373/239‑1 (MQ=255)
                                            aagcGGGACAATTCCAGCGAATGCCAGCATTCGGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAAc                                                                                                                                                                                                     <  1:898591/238‑1 (MQ=255)
                                            aaatggAACAATTCCAGCGAATTCCCGCATTCTGGATTGCTGCTTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGGCAGGTCTTAAGCAGTGGTCAGGCAATTATCAGTTCTAATACGGGGCTTAACCGCGCCATTCGTATGCGTCCGCAAAAAGAGAGCGGGGCGATTTTTAACGGCATCATCCATAAAAACg                                                                                                                                                                                                    <  2:678892/235‑1 (MQ=255)
                                                ggAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAg                                                                                                                                                                                                  <  1:199260/237‑1 (MQ=255)
                                                ggAAAAATTCCAGGGAATGCCAGGATTTGGGAAATCGGGGTTGGTGCAAACTGTATGTAACCGCTGATTGGTGTACATGAATATGGGTGGCGTTTTACCGGTAGGGAAACCAGCAGCCAGGCGTTAAGCAGGTCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAg                                                                                                                                                                                                  <  2:529862/237‑1 (MQ=255)
                                                  aaCAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACGTTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCaa                                                                                                                                                                                              <  1:440757/238‑1 (MQ=255)
                                                       ttCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATacgac                                                                                                                                                                                        <  1:346443/239‑1 (MQ=255)
                                                        ttcaggcaaTGCCAGTATTCTGGATTTCGGGTTGGGGGCGAACTGTATTTATCCGCTGAGTGGCGTACATGAATATGGCTGGCGTTTCACCGGTAGTGAACCCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATacgacg                                                                                                                                                                                       <  2:279598/232‑1 (MQ=255)
                                                        tCCAGCGAATGCCAGCATTCTGGATATCGGCGTTTGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATacgacg                                                                                                                                                                                       <  1:676534/239‑1 (MQ=255)
                                                            gggAATGCCAGCATGTTGGATATCCGGGTTTGTGCGAACTGTATTTATCCGATGATTGGCCGACATGAATCTGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCAGCCATAAAAACGAGCAATACGAcgcgtc                                                                                                                                                                                   <  1:66707/237‑3 (MQ=255)
                                                             tagaagacAGTATTGTGGATATTGGTTGTGGTGGCAACTGTATTTATCCGCTGGTTGGCGTACGTGAATATTGCTGGCGTGTGAAGTGTAGTAAACACAGCAGGCAGGCGTTAAGCAGTGTGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTGTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCAGCCATAAAAACGAGCAATACGACGCGccc                                                                                                                                                                                  <  2:855057/232‑1 (MQ=255)
                                                                 tGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCCTAAAAACGAGCAATACGACCCGACCTtgt                                                                                                                                                                              >  1:566416/1‑239 (MQ=255)
                                                                       cATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCCAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAAc                                                                                                                                                                        >  1:208403/1‑238 (MQ=255)
                                                                                   cGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATACCCCATAAAAACGGGCAAAACGACCCGACCTTTGGTAAACCCCCCATccccc                                                                                                                                                            >  2:387804/1‑237 (MQ=255)
                                                                                   cGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAAGCGTTAAGCAGTGCGCAGGCGGTTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGCGGCTTTTTAACGGCATCCTCCATAAAAACGGGCACAACGCCGCGCCCTTGGGTAACCGCCCATTCCAca                                                                                                                                                           >  2:490404/1‑238 (MQ=255)
                                                                                                    ttATTCATCCGCTGTGTGGGATCCTGGAATTTGGGTGGTTTTTTACCGGTAGGGAAACACGCAGCCGTGCTGTAAGCGGTGCGCAGGCGATTACAAGTTCTAACCGGGGGCTTAACCGCGCAATTCGTCCGCTTCGGCAAAAGGAGGGTGGGGCAATTTTTAACGGCATAATCCATAAAAACGAGCAATACGAGGCGCCCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCa                                                                                                                                          <  2:771061/238‑1 (MQ=255)
                                                                                                    tgttttACCCGGTTATTGGCTTCAGTGAATTTGGTTGGCTTTTGACCGGTAGGAAAACCAGCAGCCGGGCTTTAAGCAGGGCCAAGGCGATTACCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCa                                                                                                                                          <  1:136163/236‑1 (MQ=255)
                                                                                                     tATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACCCGACCTTGTGGAACCCGCCATTCCCCGATTCCGCCGCTGCGGccc                                                                                                                                         >  1:635164/1‑237 (MQ=255)
                                                                                                               cTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATaaaa                                                                                                                                                                                                       >  2:437991/1‑169 (MQ=255)
                                                                                                               cTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATaaaa                                                                                                                                                                                                       <  1:437991/169‑1 (MQ=255)
                                                                                                               cTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACCCGACCTTGTGTAACCCGCCATTCCACGATTCCCCCCCGCCGGCACGGGCAgggaa                                                                                                                               >  2:519236/1‑236 (MQ=255)
                                                                                                                 gATTGGCGTACAGGAATATGGCGGGCGTTTTACCGGTAGTGAAACCAGCAGCCGGGCTTTAAGCAGGGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGAAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTACCCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTg                                                                                                                             <  1:244312/239‑1 (MQ=255)
                                                                                                                                   ggggggggtttttACGGGTAGTAAACCCGGCCGCGGGGCGTTAAACAGTGCGCAGGGGATTACCAGTTTAAACCGGGGGCTAACCCGCGCCATTCGGTGGGGGGGGAAAAAAGAGAGGGGGGCGATTTTTAACGCCATCACCCATAAAAACGGGCAATCCGCCGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAAcc                                                                                                           <  2:566416/230‑1 (MQ=255)
                                                                                                                                      ctggcGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCCGGTAGTGAGCGGAAAAGCCGTAACCTgg                                                                                                        >  1:722284/1‑239 (MQ=255)
                                                                                                                                        ggcGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTAGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCCCTGGGGCACGGGCAGGAAGTGAGCGTAAACGCCGTAACCCgggg                                                                                                      >  2:67216/1‑239 (MQ=255)
                                                                                                                                         gcGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGGGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACCCGCCCTTGTGTAACCCGCCATTCCCCGCTTCCGCCCCTCCGGCCCGGGCAGGTAGGAGGCGAAAAACCCGGAACCCggggg                                                                                                     >  2:171736/1‑238 (MQ=255)
                                                                                                                                         gcGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGAAGTGAGCGTAAACGCCGTAACCTggg                                                                                                       >  1:548728/1‑237 (MQ=255)
                                                                                                                                                cccgggtggaaaaCCAGGAGCCGGGGTTTAAGCAGGGCGCGGGGGATTAGCAGTTCAAAGCGGGGGCTTACCCGCGCCAGTCGTCGGCGTGGGAAAAAAGAGAGGGGGGCGATTTTTAACGCCACCATCCAAAAAAACGAGCAATCCGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACaa                                                                                              <  2:296114/229‑1 (MQ=255)
                                                                                                                                                          aaaCCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCCCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGAcattgc                                                                                     >  1:634171/1‑233 (MQ=255)
                                                                                                                                                                                  gTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGAAGATGCACTGAACTTTGGCGGCCAGCAACAg                                                            >  1:443077/1‑239 (MQ=255)
                                                                                                                                                                                   tGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTCAGCGCAAACGCCGGAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCACCAACAgg                                                           >  1:227239/1‑239 (MQ=255)
                                                                                                                                                                                      gcAGGCTATTATCAGTTCAAAGCCGGGGTTTAACCGCGCCATTCGGCTGCGTCGGAAAAAAGAGAGTGGGGCGATTTTTAACGGCATCACCCATAAAACCGAGCAATACGACGCGCCCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAg                                                         <  1:171736/238‑1 (MQ=255)
                                                                                                                                                                                         ggCGATTAACAGTTCTAATCCGGGTCTTAACCGGGCCATTCGTCTGCGTCGGAAAAAAGAGAGGGGGGCGATTTTAAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTtg                                                      <  1:591993/238‑1 (MQ=255)
                                                                                                                                                                                                aTCAGTTATAATCCGGGGCTTAACCGCGCCATTCGTCTGCGACGGAAAAAAGAGAGGGGGGCAATTTTTAACGGCATCACCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGTGTGa                                              <  1:387804/239‑1 (MQ=255)
                                                                                                                                                                                                    gTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCACCAACAGGGGTTGTGGTGTGAAGGc                                          >  1:420660/1‑239 (MQ=255)
                                                                                                                                                                                                      tCTAAGCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGGGTGGGGCGATTTTTAACGGCATCATCAATAAAAACGAGCAATCGGCCGGGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGTGTGAAGGCg                                         <  1:490404/238‑1 (MQ=255)
                                                                                                                                                                                                          aTCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGGAGTGAGCGTAAACGCCGTAACCTGGGGCTGGACAAAGACGCTGCACTGAACTTTGGCGCCCAGCAACAGGAGGGGTGGGGTGAAGACGGGGag                                    >  2:583099/1‑238 (MQ=255)
                                                                                                                                                                                                               gggCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACCAAAACGATGACCTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGGGTGAAGGCGGTGAAGTc                                 >  1:17965/1‑237 (MQ=255)
                                                                                                                                                                                                                           gcCATTCGTTGGCTTCGGCAAAAAGAGAGGGGGGCTATTTTTAACGGCATCATCCAAAAAAACGAGCAATGCGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGTGTGAAGGCGGTGAAGTCACCTTTATCaaaaa                   <  1:712672/239‑1 (MQ=255)
                                                                                                                                                                                                                           gcCATTCGTCGGGGGCGGAAAAAAGGGAGTAGGGCGATTTTAAACGGCATCACCCAAAAAAACGAGCAATACGGCGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGTGTGAAGGCGGTGAAGTCACCTTTATCaaaaa                   <  1:67216/239‑1 (MQ=255)
                                                                                                                                                                                                                                     tGCGCCGGAAAAAAGGAGTGGGGTCTTTTTTTAACGGATCCCTCAAAAAAAGCGAGCATAGCGGCGCGCCCTTGGGAAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGTGTGAAGGCGGTGAAGTCACCTTTATCAAAAAGATGATTGaa         <  2:208403/239‑1 (MQ=255)
                                                                                                                                                                                                                                      gctgcGGAAAAAGGAGAGTGGGGCGATTTTAAACGGCATAATCCAAAAAACGGGACAATCCGCCGGGCCCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGTGTGAAGGCGGTGAAGTCACCTTTATCAAAAAGATGATTGAag        <  1:70083/235‑1 (MQ=255)
                                                                                                                                                                                                                                       cGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGTGTGAAGGCGGTGAAGTCACCTTTTTCAAAAAGATGATTGAagg       >  1:619345/1‑238 (MQ=255)
                                                                                                                                                                                                                                            gCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGGAGTGCGCGTAAACACCGTAACCTGGGGCCAAACAAAGACAAGGCCCGGACCTTTTGCGGCCACAAAAAGGGGATGTGTGGGGACGGAGGGGAGGTCACCTGTTTCAAAAACAGAATTGaaagaagaa  >  2:881792/1‑232 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                
CGGGCGGATTATATTCATCACCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGTGTGAAGGCGGTGAAGTCACCTTTATCAAAAAGATGATTGAAGAGAGCA  >  NC_000913/842578‑843052

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: