Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,946,671 A→G intergenic (+64/‑29) rrlC → / → rrfC 23S ribosomal RNA/5S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,946,6710AG100.0% 80.3 / NA 34intergenic (+64/‑29)rrlC/rrfC23S ribosomal RNA/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/19);  total (15/19)

TAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGG  >  NC_000913/3946539‑3946806
                                                                                                                                    |                                                                                                                                       
tAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAGCGCAGaa                                                                                                                                  <  2:46021/140‑1 (MQ=17)
             tAAGCACAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAGCGCAGAAGCGGTCTGATGaa                                                                                                                     >  1:141308/1‑140 (MQ=19)
                    aGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAGCGCAGAAGCGGTCTGATGAAACAGAAt                                                                                                              >  1:364363/1‑140 (MQ=21)
                     gCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATGAAATCAGAGCGGAGGAGCGGTCTGATGAAACAGAAtt                                                                                                             >  1:179154/1‑140 (MQ=255)
                      cGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAGCGCAGAAGCGGTCTGATGAAACAGAAttt                                                                                                            >  1:360172/1‑140 (MQ=21)
                             ttGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGAGTAAATCAGAGCGCAGAAGCGGTCTGATGAAACAGAATTTGCCTggc                                                                                                     <  2:664262/140‑1 (MQ=21)
                                       cGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAGCGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAgcg                                                                                           <  2:369359/140‑1 (MQ=21)
                                         gTACTAATGAACCGTGACGCTTCACCTTACACCGCAGAATATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAGCGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGcgcg                                                                                         <  2:201239/140‑1 (MQ=11)
                                         gTACTAATGAACCGTGAAGCTTGACCTTACACCGCTGAACATGTTTTGACGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAGCGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGcgcg                                                                                         <  2:179154/140‑1 (MQ=11)
                                              aaTGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAGCGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCggtggt                                                                                    <  2:360172/140‑1 (MQ=21)
                                                 gAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAGCGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTccc                                                                                 >  1:455138/1‑140 (MQ=25)
                                                 gAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAGCGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTccc                                                                                 >  2:567183/1‑140 (MQ=25)
                                                       tGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAGCGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGa                                                                           >  1:272123/1‑140 (MQ=21)
                                                         aGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAGCGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGAcc                                                                         >  1:657040/1‑140 (MQ=25)
                                                                    tACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGGTTTTCAGCCTGATACAGATTAAATCAGAGCGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAAc                                                              <  2:272123/140‑1 (MQ=21)
                                                                     aCAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAGCGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACt                                                             <  2:422554/140‑1 (MQ=21)
                                                                        aCGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAGCGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAg                                                          <  2:657040/140‑1 (MQ=21)
                                                                           ccGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAGCGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGTGGTGGTCCCACCTGACCGCATGCCGAACTCAGAAg                                                       >  2:183981/1‑140 (MQ=255)
                                                                                gATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAGCGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGaaa                                                  <  1:570700/140‑1 (MQ=18)
                                                                                gATGTTTTGGCGGATGAGAGAAGATATTCCGCCTGATCCAGATTAAATCAGAGCGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGaaa                                                  <  1:266844/140‑1 (MQ=255)
                                                                                         gCGGATGTGAGAAGATTTTCAGCCTGATACGAATTGAATCAGAGCGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAgc                                         <  2:312737/140‑1 (MQ=11)
                                                                                             aTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAGCGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGAACCCATGCAGAATTCGGAAGTGCAACGCCGTAGCGCCg                                     >  1:254143/1‑140 (MQ=11)
                                                                                                agagaAGATTTTCAGCCTGATACAGATTAAATCAGAGCGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATg                                  <  2:433536/140‑1 (MQ=17)
                                                                                                  agaagaTTTTCAGCCTGATACAGATTAAATCAGAGCGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGt                                <  2:364363/140‑1 (MQ=17)
                                                                                                               ccTGATACAGATTAAATCAGAGCGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGtctc                   >  1:591949/1‑140 (MQ=18)
                                                                                                                   aTACAGATTAAATCAGAGCGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCa               <  2:455138/140‑1 (MQ=17)
                                                                                                                   aTACAGATTAAATCACAGCGCACAAACGGTCTGATGATACAGAATTTGTCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCa               <  1:192486/140‑1 (MQ=11)
                                                                                                                        gATTAAATCAGAGCGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCga          <  1:31122/140‑1 (MQ=17)
                                                                                                                         aTTAAATCAGAGCGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCgag         >  1:563968/1‑140 (MQ=18)
                                                                                                                            aaaTCAGAGCGCAGAAGTGGTCTGATGAAACAGAATTTGGCTGGCGGCTGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGt      <  2:254143/140‑1 (MQ=255)
                                                                                                                             aaTCAGAGCGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTa     >  2:331922/1‑140 (MQ=18)
                                                                                                                             aaTCAGAGCGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGAAGCTCTGACGGTAGTGTGGGGGGTCCCCATTCGAGAGTa     >  2:327531/1‑140 (MQ=11)
                                                                                                                              aTCAGAGCGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGAGCCGATGGTAGTGTGGGGTCTACCCATGCGAGAGTAg    >  1:228786/1‑140 (MQ=255)
                                                                                                                                cAGAGCGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAggg  <  2:141308/140‑1 (MQ=17)
                                                                                                                                    |                                                                                                                                       
TAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGG  >  NC_000913/3946539‑3946806

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑