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

CGGGCGGATTATATTCATCACCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCATCCATAAAAACGAG‑CAATACGACGC‑GACCTTGTGT‑AACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGTGTGAAGGCGGTGAAGTCACCTTTATCAAAAAGATGATTGAAGAGAGCA  >  NC_000913/842578‑843052
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 
cgggcGGATTTATTTCATCCCCTTGCCGATTTCCTGGCAGAAGCAAGCGGAAAAATTCCAGCGAATGCCATCATTATGGATACTGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGGGTTTTA‑CCGGTAGTGAAACCAGCCGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGc                                                                                                                                                                                                                                                   <  1:804806/238‑1 (MQ=255)
  ggcGGATTATATTCATCACCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCaaa                                                                                                                                                                                                                                                >  2:279598/1‑239 (MQ=255)
  ggcGGATTATATTCATCACCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCCCCATTGGTCT‑GCGTCGGCaaa                                                                                                                                                                                                                                                >  2:166844/1‑239 (MQ=255)
              tcatccCCTTGCCTATTTACTGGCAGAAGCAAGGGGAACAATTCCCGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTgggg                                                                                                                                                                                                                                    <  1:269627/239‑1 (MQ=255)
              tcaacaCCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTgggg                                                                                                                                                                                                                                    <  2:264163/239‑1 (MQ=255)
                 tcaCCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGACTGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATACGCTGATTGGCGTACATGACTATGGCTGGAGTTTTA‑CCGGTAGGGAAGACAGCGGCCCGGGGTTAAGCAGATGGCAGGCGATTTTCGGATATAAACCAGGGCATAACAGGCCCATTCGGCT‑GCGGCGGAACAAAGAGAGGGGGG‑CGa                                                                                                                                                                                                                                >  1:801309/1‑239 (MQ=255)
                 tcaCCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGGGCAGGGGGTTATCAGTTCTAATCCGGGGTTTAACACCGCCATTC‑TCTGGCTGCGGCaaaaaaaaaaggggggga                                                                                                                                                                                                                                 >  1:136163/1‑226 (MQ=255)
                   aCCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGGGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGGCCCATTCGTCT‑GCGTCGGCAAAAAAAGAGTGGGG‑CGAtt                                                                                                                                                                                                                              >  2:678892/1‑239 (MQ=255)
                       tGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGGGTTAAGCAGTGCGCAGGCGCTTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGAAAAAAGAGAGGGGGG‑CGAGTTTTa                                                                                                                                                                                                                          >  1:582767/1‑239 (MQ=255)
                        tgcgttttCCTGGCAGAAGGAAGGGGAACAATTCCAGCGAATGCCAGCATTCTGGAGAGCGGTGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCGGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTaa                                                                                                                                                                                                                         <  1:861553/234‑1 (MQ=255)
                        gCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTaa                                                                                                                                                                                                                         <  2:582767/239‑1 (MQ=255)
                                         agcaagcGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGTTGATTGGCGTACATGAATATGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCATCCATaaaa                                                                                                                                                                                                        <  2:632373/239‑1 (MQ=255)
                                            aagcGGGACAATTCCAGCGAATGCCAGCATTCGGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCATCCATAAAAAc                                                                                                                                                                                                      <  2:898591/238‑1 (MQ=255)
                                            aaatggAACAATTCCAGCGAATTCCCGCATTCTGGATTGCTGCTTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGGCAGGTCTTAAGCAGTGGTCAGGCAATTATCAGTTCTAATACGGGGCTTAACCGCGCCATTCGTAT‑GCGTCCGCAAAAAGAGAGCGGGG‑CGATTTTTAACGGCATCATCCATAAAAACg                                                                                                                                                                                                     <  1:678892/235‑1 (MQ=255)
                                                ggAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCATCCATAAAAACGAg                                                                                                                                                                                                   <  2:199260/237‑1 (MQ=255)
                                                ggAAAAATTCCAGGGAATGCCAGGATTTGGGAAATCGGGGTTGGTGCAAACTGTATGTAACCGCTGATTGGTGTACATGAATATGGGTGGCGTTTTA‑CCGGTAGGGAAACCAGCAGCCAGGCGTTAAGCAGGTCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCATCCATAAAAACGAg                                                                                                                                                                                                   <  1:529862/237‑1 (MQ=255)
                                                  aaCAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACGTTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCATCCATAAAAACGAG‑Caa                                                                                                                                                                                               <  2:440757/238‑1 (MQ=255)
                                                       ttCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCATCCATAAAAACGAG‑CAATacgac                                                                                                                                                                                         <  2:346443/239‑1 (MQ=255)
                                                        ttcaggcaaTGCCAGTATTCTGGATTTCGGGTTGGGGGCGAACTGTATTTATCCGCTGAGTGGCGTACATGAATATGGCTGGCGTTTCA‑CCGGTAGTGAACCCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCATCCATAAAAACGAG‑CAATacgacg                                                                                                                                                                                        <  1:279598/232‑1 (MQ=255)
                                                        tCCAGCGAATGCCAGCATTCTGGATATCGGCGTTTGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCATCCATAAAAACGAG‑CAATacgacg                                                                                                                                                                                        <  2:676534/239‑1 (MQ=255)
                                                            gggAATGCCAGCATGTTGGATATCCGGGTTTGTGCGAACTGTATTTATCCGATGATTGGCCGACATGAATCTGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCAGCCATAAAAACGAG‑CAATACGAcgc‑gtc                                                                                                                                                                                   <  2:66707/237‑3 (MQ=255)
                                                             tagaagacAGTATTGTGGATATTGGTTGTGGTGGCAACTGTATTTATCCGCTGGTTGGCGTACGTGAATATTGCTGGCGTGTGA‑AGTGTAGTAAACACAGCAGGCAGGCGTTAAGCAGTGTGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTGT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCAGCCATAAAAACGAG‑CAATACGACGC‑Gccc                                                                                                                                                                                  <  1:855057/232‑1 (MQ=255)
                                                                 tGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTCTCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGAGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGTGCCATTCGTCT‑GCGTCGGCAAAAAGGGGGTGGGG‑CGATTTTTAACGGCATCTACCCTAAAAACGAC‑CAAAacaaacggcccttgg                                                                                                                                                                               >  1:860176/1‑226 (MQ=255)
                                                                 tGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCATCCCTAAAAACGAG‑CAATACGACCC‑GACCTtgt                                                                                                                                                                              >  2:566416/1‑239 (MQ=255)
                                                                    ctgttttctgtttttCGGCTGTGGTGCGAATTGATTTTACCGGGTGTTTGGCGTCAATGAATGTGTGTGGGGTTTTC‑CCGGTAGTGAAACAAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGTTTACCAGTTCTAACCCGGGGCTAAACCGCGCTATTCGTCG‑GCGTCGGCAAAAAGAGAGTGGGG‑CGTTTTTTAACGGCATCACCCAAAAAAACGAG‑CAATACGCCGC‑GACCTTGTGT‑a                                                                                                                                                                          <  1:325393/225‑1 (MQ=255)
                                                                       cATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCCAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCATCCATAAAAACGAG‑CAATACGACGC‑GACCTTGTGT‑AAc                                                                                                                                                                        >  2:208403/1‑238 (MQ=255)
                                                                                   cGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATACCCCATAAAAACGGG‑CAAAACGACCC‑GACCTTTGGTAAACCCCCCATccccc                                                                                                                                                            >  1:387804/1‑237 (MQ=255)
                                                                                   cGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAAGCGTTAAGCAGTGCGCAGGCGGTTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGC‑GGCTTTTTAACGGCATCCTCCATAAAAACGGG‑CACAACGCCGC‑GCCCTTGGGT‑AACCGCCCATTCCAca                                                                                                                                                           >  1:490404/1‑238 (MQ=255)
                                                                                              gtacttttttttttcggttgtttGGCGTAAATG‑ATATGGGTGGGGTTTTACCCGGTAGTGAAACCCGCCGCCGGGCGTAAAGCGGGGCGAGGGGGTTTATCATTTTTATTCCGGGGCTTAACCGCCCCATTCGTCG‑GCGTCGGCAAAAAGAGGGGGGGG‑CGTTTTTTAACGGCATCACCCCAAAAAACGAG‑CAATACGACGC‑GACCTTGTGT‑AACCCGCCATTCCACGATTCCGCCGCt                                                                                                                                                <  1:713585/218‑1 (MQ=255)
                                                                                                    ttATTCATCCGCTGTGTGGGATCCTGGAATTTGGGTGGTTTTTTA‑CCGGTAGGGAAACACGCAGCCGTGCTGTAAGCGGTGCGCAGGCGATTACAAGTTCTAACCGGGGGCTTAACCGCGCAATTCGTCC‑GCTTCGGCAAAAGGAGGGTGGGG‑CAATTTTTAACGGCATAATCCATAAAAACGAG‑CAATACGAGGC‑GCCCTTGTGT‑AACCCGCCATTCCACGATTCCGCCGCTGCGGCa                                                                                                                                          <  1:771061/238‑1 (MQ=255)
                                                                                                    tgttttACCCGGTTATTGGCTTCAGTGAATTTGGTTGGCTTTTGA‑CCGGTAGGAAAACCAGCAGCCGGGCTTTAAGCAGGGCCAAGGCGATTACCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCATCCATAAAAACGAG‑CAATACGACGC‑GACCTTGTGT‑AACCCGCCATTCCACGATTCCGCCGCTGCGGCa                                                                                                                                          <  2:136163/236‑1 (MQ=255)
                                                                                                     tATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCATCCATAAAAACGAG‑CAATACGACCC‑GACCTTGTGG‑AACCCGCCATTCCCCGATTCCGCCGCTGCGGccc                                                                                                                                         >  2:635164/1‑237 (MQ=255)
                                                                                                               cTGATTGGCGTACATGAATATGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTCATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCCAAAAGAGCGGGGGG‑CGATTTTTAACGGCATACTCCCTAAAACACAG‑CAAAACGAACC‑GACCCTTGGT‑AACCCACCCTTCCCAAATTCCGCCGCTGCggacgcggccgggaa                                                                                                                               >  1:58866/1‑226 (MQ=255)
                                                                                                               cTGATTGGCGTACATGAATATGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCATCCATaaaa                                                                                                                                                                                                        <  2:437991/169‑1 (MQ=255)
                                                                                                               cTGATTGGCGTACATGAATATGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCATCCATaaaa                                                                                                                                                                                                        >  1:437991/1‑169 (MQ=255)
                                                                                                               cTGATTGGCGTACATGAATATGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCATCCATAAAAACGAG‑CAATACGACCC‑GACCTTGTGT‑AACCCGCCATTCCACGATTCCCCCCCGCCGGCACGGGCAgggaa                                                                                                                               >  1:519236/1‑236 (MQ=255)
                                                                                                               cTGATTGGCGTACATGAATATGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGGGGGG‑CGATTTTTAACGGCCTCATCCGTAAAAACGAG‑CAA‑AAGACCCGGCCCTTGTGT‑AACCCGCCATTCCCACATTCccccgcccgggaccgggggggagg                                                                                                                               >  1:788222/1‑221 (MQ=255)
                                                                                                                 gATTGGCGTACAGGAATATGGCGGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCGGGCTTTAAGCAGGGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGAAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCATCCATAAAAACGAG‑CAATACGACGC‑GACCTTGTGT‑ACCCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTg                                                                                                                             <  2:244312/239‑1 (MQ=255)
                                                                                                                                   ggggggggtttttA‑CGGGTAGTAAACCCGGCCGCGGGGCGTTAAACAGTGCGCAGGGGATTACCAGTTTAAACCGGGGGCTAACCCGCGCCATTCGGTG‑GGGGGGGAAAAAAGAGAGGGGGG‑CGATTTTTAACGCCATCACCCATAAAAACGGG‑CAATCCGCCGC‑GACCTTGTGT‑AACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAAcc                                                                                                           <  1:566416/230‑1 (MQ=255)
                                                                                                                                      ctggcGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCATCCATAAAAACGAG‑CAATACGACGC‑GACCTTGTGT‑AACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCCGGTAGTGAGCGGAAAAGCCGTAACCTgg                                                                                                        >  2:722284/1‑239 (MQ=255)
                                                                                                                                        ggcGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTAGGG‑CGATTTTTAACGGCATCATCCATAAAAACGAG‑CAATACGACGC‑GACCTTGTGT‑AACCCGCCATTCCACGATTCCGCCCCTGGGGCACGGGCAGGAAGTGAGCGTAAACGCCGTAACCCgggg                                                                                                      >  1:67216/1‑239 (MQ=255)
                                                                                                                                         gcGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGGGAGTGGGG‑CGATTTTTAACGGCATCATCCATAAAAACGAG‑CAATACGACCC‑GCCCTTGTGT‑AACCCGCCATTCCCCGCTTCCGCCCCTCCGGCCCGGGCAGGTAGGAGGCGAAAAACCCGGAACCCggggg                                                                                                     >  1:171736/1‑238 (MQ=255)
                                                                                                                                         gcGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCATCCATAAAAACGAG‑CAATACGACGC‑GACCTTGTGT‑AACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGAAGTGAGCGTAAACGCCGTAACCTggg                                                                                                       >  2:548728/1‑237 (MQ=255)
                                                                                                                                                 cccgggtggaaaaCCAGGAGCCGGGGTTTAAGCAGGGCGCGGGGGATTAGCAGTTCAAAGCGGGGGCTTACCCGCGCCAGTCGTCG‑GCGTGGGAAAAAAGAGAGGGGGG‑CGATTTTTAACGCCACCATCCAAAAAAACGAG‑CAATCCGACGC‑GACCTTGTGT‑AACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACaa                                                                                              <  1:296114/229‑1 (MQ=255)
                                                                                                                                                           aaaCCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCATCCATAAAAACGAG‑CAATACGACGC‑GACCTTGTGT‑AACCCGCCATTCCACGATTCCGCCCCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGAcattgc                                                                                     >  2:634171/1‑233 (MQ=255)
                                                                                                                                                                                   gTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCATCCATAAAAACGAG‑CAATACGACGC‑GACCTTGTGT‑AACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGAAGATGCACTGAACTTTGGCGGCCAGCAACAg                                                            >  2:443077/1‑239 (MQ=255)
                                                                                                                                                                                    tGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCATCCATAAAAACGAG‑CAATACGACGC‑GACCTTGTGT‑AACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTCAGCGCAAACGCCGGAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCACCAACAgg                                                           >  2:227239/1‑239 (MQ=255)
                                                                                                                                                                                       gcAGGCTATTATCAGTTCAAAGCCGGGGTTTAACCGCGCCATTCGGCT‑GCGTCGGAAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCACCCATAAAACCGAG‑CAATACGACGC‑GCCCTTGTGT‑AACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAg                                                         <  2:171736/238‑1 (MQ=255)
                                                                                                                                                                                          ggCGATTAACAGTTCTAATCCGGGTCTTAACCGGGCCATTCGTCT‑GCGTCGGAAAAAAGAGAGGGGGG‑CGATTTTAAACGGCATCATCCATAAAAACGAG‑CAATACGACGC‑GACCTTGTGT‑AACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTtg                                                      <  2:591993/238‑1 (MQ=255)
                                                                                                                                                                                                 aTCAGTTATAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGACGGAAAAAAGAGAGGGGGG‑CAATTTTTAACGGCATCACCCATAAAAACGAG‑CAATACGACGC‑GACCTTGTGT‑AACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGTGTGa                                              <  2:387804/239‑1 (MQ=255)
                                                                                                                                                                                                     gTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCATCCATAAAAACGAG‑CAATACGACGC‑GACCTTGTGT‑AACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCACCAACAGGGGTTGTGGTGTGAAGGc                                          >  2:420660/1‑239 (MQ=255)
                                                                                                                                                                                                       tCTAAGCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGGGTGGGG‑CGATTTTTAACGGCATCATCAATAAAAACGAG‑CAATCGGCCGG‑GACCTTGTGT‑AACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGTGTGAAGGCg                                         <  2:490404/238‑1 (MQ=255)
                                                                                                                                                                                                           aTCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCATCCATAAAAACGAG‑CAATACGACGC‑GACCTTGTGT‑AACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGGAGTGAGCGTAAACGCCGTAACCTGGGGCTGGACAAAGACGCTGCACTGAACTTTGGCGCCCAGCAACAGGAGGGGTGGGGTGAAGACGGGGag                                    >  1:583099/1‑238 (MQ=255)
                                                                                                                                                                                                                gggCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCATCCATAAAAACGAG‑CAATACGACGC‑GACCTTGTGT‑AACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACCAAAACGATGACCTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGGGTGAAGGCGGTGAAGTc                                 >  2:17965/1‑237 (MQ=255)
                                                                                                                                                                                                                            gcCATTCGTTG‑GCTTCGGCAAAAAGAGAGGGGGG‑CTATTTTTAACGGCATCATCCAAAAAAACGAG‑CAATGCGACGC‑GACCTTGTGT‑AACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGTGTGAAGGCGGTGAAGTCACCTTTATCaaaaa                   <  2:712672/239‑1 (MQ=255)
                                                                                                                                                                                                                            gcCATTCGTCG‑GGGGCGGAAAAAAGGGAGTAGGG‑CGATTTTAAACGGCATCACCCAAAAAAACGAG‑CAATACGGCGC‑GACCTTGTGT‑AACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGTGTGAAGGCGGTGAAGTCACCTTTATCaaaaa                   <  2:67216/239‑1 (MQ=255)
                                                                                                                                                                                                                                      t‑GCGCCGGAAAAAAG‑GAGTGGGGTCTTTTTTTAACGGATCCCTCAAAAAAAGCGAG‑CATAGCGGCGC‑GCCCTTGGGA‑AACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGTGTGAAGGCGGTGAAGTCACCTTTATCAAAAAGATGATTGaa         <  1:208403/239‑1 (MQ=255)
                                                                                                                                                                                                                                        gctgcGGAAAAAGGAGAGTGGGG‑CGATTTTAAACGGCATAATCCA‑AAAAACGGGACAATCCGCCGG‑GCCCTTGTGT‑AACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGTGTGAAGGCGGTGAAGTCACCTTTATCAAAAAGATGATTGAag        <  2:70083/235‑1 (MQ=255)
                                                                                                                                                                                                                                         cGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCATCCATAAAAACGAG‑CAATACGACGC‑GACCTTGTGT‑AACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGTGTGAAGGCGGTGAAGTCACCTTTTTCAAAAAGATGATTGAagg       >  2:619345/1‑238 (MQ=255)
                                                                                                                                                                                                                                              gCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCATCCATAAAAACGAG‑CAATACGACGC‑GACCTTGTGT‑AACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGGAGTGCGCGTAAACACCGTAACCTGGGGCCAAACAAAGACAAGGCCCGGACCTTTTGCGGCCACAAAAAGGGGATGTGTGGGGACGGAGGGGAGGTCACCTGTTTCAAAAACAGAATTGaaagaagaa  >  1:881792/1‑232 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 
CGGGCGGATTATATTCATCACCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTA‑CCGGTAGTGAAACCAGCAGCCAGGCGTTAAGCAGTGCGCAGGCGATTATCAGTTCTAATCCGGGGCTTAACCGCGCCATTCGTCT‑GCGTCGGCAAAAAGAGAGTGGGG‑CGATTTTTAACGGCATCATCCATAAAAACGAG‑CAATACGACGC‑GACCTTGTGT‑AACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAACAAAGACGATGCACTGAACTTTGGCGGCCAGCAACAGGAGTTGTGGTGTGAAGGCGGTGAAGTCACCTTTATCAAAAAGATGATTGAAGAGAGCA  >  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: ‑