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

CGGGCGGATTATATTCATCACCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGTGTGAAGGCGGTGAAGTCACCTTTATCAAAAAGATGATTGAAGAGAGCA  >  NC_000913/842578‑843052
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 
cgggcGGATTTATTTCATCCCCTTGCCGATTTCCTGGCAGAAGCAAGCGGAAAAATTCCAGCGAATGCCATCATTATGGATACTGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGGGTTTTACCGGTAGTGAAACCAGCCGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGc                                                                                                                                                                                                                                                   <  1:804806/238‑1 (MQ=255)
  ggcGGATTATATTCATCACCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCaaa                                                                                                                                                                                                                                                >  2:279598/1‑239 (MQ=255)
  ggcGGATTATATTCATCACCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCCCCATTGGTCTGCGTCGGCaaa                                                                                                                                                                                                                                                >  2:166844/1‑239 (MQ=255)
              tcatccCCTTGCCTATTTACTGGCAGAAGCAAGGGGAACAATTCCCGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTgggg                                                                                                                                                                                                                                    <  1:269627/239‑1 (MQ=255)
              tcaacaCCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTgggg                                                                                                                                                                                                                                    <  2:264163/239‑1 (MQ=255)
                 tcaCCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGACTGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATACGCTGATTGGCGTACATGACTATGGCTGGAGTTTTACCGGTAGGGAAGACAGCGGCCCGGGGTTAAGCAGATGGCAGGCGATTTTCGGATATAAACCAGGGCATAACAGGCCCATTCGGCTGCGGCGGAACAAAGAGAGGGGGGCGa                                                                                                                                                                                                                                >  1:801309/1‑239 (MQ=255)
                 tcaCCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGGGCAGGGGGTTATCAGTTCTAATCCGGGGTTTAACACCGCCATTCTCTGGCTGCGGCaaaaaaaaaaggggggga                                                                                                                                                                                                                                 >  1:136163/1‑226 (MQ=255)
                   aCCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGGGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGGCCCATTCGTCTGCGTCGGCAAAAAAAGAGTGGGGCGAtt                                                                                                                                                                                                                              >  2:678892/1‑239 (MQ=255)
                       tGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGGGTTAAGCAGTGCGCAGGCGCTTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGAAAAAAGAGAGGGGGGCGAGTTTTa                                                                                                                                                                                                                          >  1:582767/1‑239 (MQ=255)
                        tgcgttttCCTGGCAGAAGGAAGGGGAACAATTCCAGCGAATGCCAGCATTCTGGAGAGCGGTGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCGGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTaa                                                                                                                                                                                                                         <  1:861553/234‑1 (MQ=255)
                        gCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTaa                                                                                                                                                                                                                         <  2:582767/239‑1 (MQ=255)
                                         agcaagcGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGTTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATaaaa                                                                                                                                                                                                        <  2:632373/239‑1 (MQ=255)
                                            aagcGGGACAATTCCAGCGAATGCCAGCATTCGGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAAc                                                                                                                                                                                                      <  2:898591/238‑1 (MQ=255)
                                            aaatggAACAATTCCAGCGAATTCCCGCATTCTGGATTGCTGCTTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGGCAGGTCTTAAGCAGTGGTCAGGCAATTATCAGTTCTAATACGGGGCTTAACCGCGCCATTCGTATGCGTCCGCAAAAAGAGAGCGGGGCGATTTTTAACGGCATCATCCATAAAAACg                                                                                                                                                                                                     <  1:678892/235‑1 (MQ=255)
                                                ggAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAg                                                                                                                                                                                                   <  2:199260/237‑1 (MQ=255)
                                                ggAAAAATTCCAGGGAATGCCAGGATTTGGGAAATCGGGGTTGGTGCAAACTGTATGTAACCGCTGATTGGTGTACATGAATATGGGTGGCGTTTTACCGGTAGGGAAACCAGCAGCCAGGCGTTAAGCAGGTCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAg                                                                                                                                                                                                   <  1:529862/237‑1 (MQ=255)
                                                  aaCAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACGTTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCaa                                                                                                                                                                                               <  2:440757/238‑1 (MQ=255)
                                                       ttCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATacgac                                                                                                                                                                                         <  2:346443/239‑1 (MQ=255)
                                                        ttcaggcaaTGCCAGTATTCTGGATTTCGGGTTGGGGGCGAACTGTATTTATCCGCTGAGTGGCGTACATGAATATGGCTGGCGTTTCACCGGTAGTGAACCCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATacgacg                                                                                                                                                                                        <  1:279598/232‑1 (MQ=255)
                                                        tCCAGCGAATGCCAGCATTCTGGATATCGGCGTTTGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATacgacg                                                                                                                                                                                        <  2:676534/239‑1 (MQ=255)
                                                            gggAATGCCAGCATGTTGGATATCCGGGTTTGTGCGAACTGTATTTATCCGATGATTGGCCGACATGAATCTGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCAGCCATAAAAACGAGCAATACGAcgcgtc                                                                                                                                                                                   <  2:66707/237‑3 (MQ=255)
                                                             tagaagacAGTATTGTGGATATTGGTTGTGGTGGCAACTGTATTTATCCGCTGGTTGGCGTACGTGAATATTGCTGGCGTGTGAAGTGTAGTAAACACAGCAGGCAGGCGTTAAGCAGTGTGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTGTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCAGCCATAAAAACGAGCAATACGACGCGccc                                                                                                                                                                                  <  1:855057/232‑1 (MQ=255)
                                                                 tGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTCTCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGAGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGTGCCATTCGTCTGCGTCGGCAAAAAGGGGGTGGGGCGATTTTTAACGGCATCTACCCTAAAAACGACCAAAacaaacggcccttgg                                                                                                                                                                               >  1:860176/1‑226 (MQ=255)
                                                                 tGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCCTAAAAACGAGCAATACGACCCGACCTtgt                                                                                                                                                                              >  2:566416/1‑239 (MQ=255)
                                                                    ctgttttctgtttttCGGCTGTGGTGCGAATTGATTTTACCGGGTGTTTGGCGTCAATGAATGTGTGTGGGGTTTTCCCGGTAGTGAAACAAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGTTTACCAGTTCTAACCCGGGGCTAAACCGCGCTATTCGTCGGCGTCGGCAAAAAGAGAGTGGGGCGTTTTTTAACGGCATCACCCAAAAAAACGAGCAATACGCCGCGACCTTGTGTa                                                                                                                                                                          <  1:325393/225‑1 (MQ=255)
                                                                       cATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCCAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAAc                                                                                                                                                                        >  2:208403/1‑238 (MQ=255)
                                                                                   cGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATACCCCATAAAAACGGGCAAAACGACCCGACCTTTGGTAAACCCCCCATccccc                                                                                                                                                            >  1:387804/1‑237 (MQ=255)
                                                                                   cGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAAGCGTTAAGCAGTGCGCAGGCGGTTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGCGGCTTTTTAACGGCATCCTCCATAAAAACGGGCACAACGCCGCGCCCTTGGGTAACCGCCCATTCCAca                                                                                                                                                           >  1:490404/1‑238 (MQ=255)
                                                                                              gtacttttttttttcggttgtttGGCGTAAATGATATGGGTGGGGTTTTACCCGGTAGTGAAACCCGCCGCCGGGCGTAAAGCGGGGCGAGGGGGTTTATCATTTTTATTCCGGGGCTTAACCGCCCCATTCGTCGGCGTCGGCAAAAAGAGGGGGGGGCGTTTTTTAACGGCATCACCCCAAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCt                                                                                                                                                <  1:713585/218‑1 (MQ=255)
                                                                                                    ttATTCATCCGCTGTGTGGGATCCTGGAATTTGGGTGGTTTTTTACCGGTAGGGAAACACGCAGCCGTGCTGTAAGCGGTGCGCAGGCGATTACAAGTTCTAACCGGGGGCTTAACCGCGCAATTCGTCCGCTTCGGCAAAAGGAGGGTGGGGCAATTTTTAACGGCATAATCCATAAAAACGAGCAATACGAGGCGCCCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCa                                                                                                                                          <  1:771061/238‑1 (MQ=255)
                                                                                                    tgttttACCCGGTTATTGGCTTCAGTGAATTTGGTTGGCTTTTGACCGGTAGGAAAACCAGCAGCCGGGCTTTAAGCAGGGCCAAGGCGATTACCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCa                                                                                                                                          <  2:136163/236‑1 (MQ=255)
                                                                                                     tATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACCCGACCTTGTGGAACCCGCCATTCCCCGATTCCGCCGCTGCGGccc                                                                                                                                         >  2:635164/1‑237 (MQ=255)
                                                                                                               cTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTCATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCCAAAAGAGCGGGGGGCGATTTTTAACGGCATACTCCCTAAAACACAGCAAAACGAACCGACCCTTGGTAACCCACCCTTCCCAAATTCCGCCGCTGCggacgcggccgggaa                                                                                                                               >  1:58866/1‑226 (MQ=255)
                                                                                                               cTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATaaaa                                                                                                                                                                                                        <  2:437991/169‑1 (MQ=255)
                                                                                                               cTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATaaaa                                                                                                                                                                                                        >  1:437991/1‑169 (MQ=255)
                                                                                                               cTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACCCGACCTTGTGTAACCCGCCATTCCACGATTCCCCCCCGCCGGCACGGGCAgggaa                                                                                                                               >  1:519236/1‑236 (MQ=255)
                                                                                                               cTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGGGGGGCGATTTTTAACGGCCTCATCCGTAAAAACGAGCAAAAGACCCGGCCCTTGTGTAACCCGCCATTCCCACATTCccccgcccgggaccgggggggagg                                                                                                                               >  1:788222/1‑221 (MQ=255)
                                                                                                                 gATTGGCGTACAGGAATATGGCGGGCGTTTTACCGGTAGTGAAACCAGCAGCCGGGCTTTAAGCAGGGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGAAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTACCCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTg                                                                                                                             <  2:244312/239‑1 (MQ=255)
                                                                                                                                   ggggggggtttttACGGGTAGTAAACCCGGCCGCGGGGCGTTAAACAGTGCGCAGGGGATTACCAGTTTAAACCGGGGGCTAACCCGCGCCATTCGGTGGGGGGGGAAAAAAGAGAGGGGGGCGATTTTTAACGCCATCACCCATAAAAACGGGCAATCCGCCGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAAcc                                                                                                           <  1:566416/230‑1 (MQ=255)
                                                                                                                                      ctggcGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCCGGTAGTGAGCGGAAAAGCCGTAACCTgg                                                                                                        >  2:722284/1‑239 (MQ=255)
                                                                                                                                        ggcGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTAGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCCCTGGGGCACGGGCAGGAAGTGAGCGTAAACGCCGTAACCCgggg                                                                                                      >  1:67216/1‑239 (MQ=255)
                                                                                                                                         gcGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGGGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACCCGCCCTTGTGTAACCCGCCATTCCCCGCTTCCGCCCCTCCGGCCCGGGCAGGTAGGAGGCGAAAAACCCGGAACCCggggg                                                                                                     >  1:171736/1‑238 (MQ=255)
                                                                                                                                         gcGTTTTACCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGAAGTGAGCGTAAACGCCGTAACCTggg                                                                                                       >  2:548728/1‑237 (MQ=255)
                                                                                                                                                 cccgggtggaaaaCCAGGAGCCGGGGTTTAAGCAGGGCGCGGGGGATTAGCAGTTCAAAGCGGGGGCTTACCCGCGCCAGTCGTCGGCGTGGGAAAAAAGAGAGGGGGGCGATTTTTAACGCCACCATCCAAAAAAACGAGCAATCCGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACaa                                                                                              <  1:296114/229‑1 (MQ=255)
                                                                                                                                                           aaaCCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCCCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGAcattgc                                                                                     >  2:634171/1‑233 (MQ=255)
                                                                                                                                                                                   gTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGAAGATGCACTGAACTTTGGCGGCCAGCAACAg                                                            >  2:443077/1‑239 (MQ=255)
                                                                                                                                                                                    tGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTCAGCGCAAACGCCGGAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCACCAACAgg                                                           >  2:227239/1‑239 (MQ=255)
                                                                                                                                                                                       gcAGGCTATTATCAGTTCAAAGCCGGGGTTTAACCGCGCCATTCGGCTGCGTCGGAAAAAAGAGAGTGGGGCGATTTTTAACGGCATCACCCATAAAACCGAGCAATACGACGCGCCCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAg                                                         <  2:171736/238‑1 (MQ=255)
                                                                                                                                                                                          ggCGATTAACAGTTCTAATCCGGGTCTTAACCGGGCCATTCGTCTGCGTCGGAAAAAAGAGAGGGGGGCGATTTTAAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTtg                                                      <  2:591993/238‑1 (MQ=255)
                                                                                                                                                                                                 aTCAGTTATAATCCGGGGCTTAACCGCGCCATTCGTCTGCGACGGAAAAAAGAGAGGGGGGCAATTTTTAACGGCATCACCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGTGTGa                                              <  2:387804/239‑1 (MQ=255)
                                                                                                                                                                                                     gTTCTAATCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCACCAACAGGGGTTGTGGTGTGAAGGc                                          >  2:420660/1‑239 (MQ=255)
                                                                                                                                                                                                       tCTAAGCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGGGTGGGGCGATTTTTAACGGCATCATCAATAAAAACGAGCAATCGGCCGGGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGTGTGAAGGCg                                         <  2:490404/238‑1 (MQ=255)
                                                                                                                                                                                                           aTCCGGGGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGGAGTGAGCGTAAACGCCGTAACCTGGGGCTGGACAAAGACGCTGCACTGAACTTTGGCGCCCAGCAACAGGAGGGGTGGGGTGAAGACGGGGag                                    >  1:583099/1‑238 (MQ=255)
                                                                                                                                                                                                                gggCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACCAAAACGATGACCTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGGGTGAAGGCGGTGAAGTc                                 >  2:17965/1‑237 (MQ=255)
                                                                                                                                                                                                                            gcCATTCGTTGGCTTCGGCAAAAAGAGAGGGGGGCTATTTTTAACGGCATCATCCAAAAAAACGAGCAATGCGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGTGTGAAGGCGGTGAAGTCACCTTTATCaaaaa                   <  2:712672/239‑1 (MQ=255)
                                                                                                                                                                                                                            gcCATTCGTCGGGGGCGGAAAAAAGGGAGTAGGGCGATTTTAAACGGCATCACCCAAAAAAACGAGCAATACGGCGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGTGTGAAGGCGGTGAAGTCACCTTTATCaaaaa                   <  2:67216/239‑1 (MQ=255)
                                                                                                                                                                                                                                      tGCGCCGGAAAAAAGGAGTGGGGTCTTTTTTTAACGGATCCCTCAAAAAAAGCGAGCATAGCGGCGCGCCCTTGGGAAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGTGTGAAGGCGGTGAAGTCACCTTTATCAAAAAGATGATTGaa         <  1:208403/239‑1 (MQ=255)
                                                                                                                                                                                                                                        gctgcGGAAAAAGGAGAGTGGGGCGATTTTAAACGGCATAATCCAAAAAACGGGACAATCCGCCGGGCCCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGTGTGAAGGCGGTGAAGTCACCTTTATCAAAAAGATGATTGAag        <  2:70083/235‑1 (MQ=255)
                                                                                                                                                                                                                                         cGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGTGTGAAGGCGGTGAAGTCACCTTTTTCAAAAAGATGATTGAagg       >  2:619345/1‑238 (MQ=255)
                                                                                                                                                                                                                                              gCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGGAGTGCGCGTAAACACCGTAACCTGGGGCCAAACAAAGACAAGGCCCGGACCTTTTGCGGCCACAAAAAGGGGATGTGTGGGGACGGAGGGGAGGTCACCTGTTTCAAAAACAGAATTGaaagaagaa  >  1: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 < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: