Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,628,807 T→G intergenic (+260/+48) nadR → / ← ettA DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR/energy‑dependent translational throttle protein EttA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,628,8070TG92.3% 122.9 / 0.1 39intergenic (+260/+48)nadR/ettADNA‑binding transcriptional repressor/NMN adenylyltransferase NadR/energy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base G (23/13);  total (25/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

AAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCGGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATCCAGAATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATC  >  NC_000913/4628565‑4629041
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          
aaaGGCGTGGTAAAAAATAAGTTTTTTGATGAGTTGAAACCCTCTCGGGAAGAGTCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACAAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCTTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGc                                                                                                                                                                                                                                   <  1:128432/251‑1 (MQ=255)
   ggCGAGGTATAAAATAAGTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGTTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCgg                                                                                                                                                                                                                                <  1:412455/251‑1 (MQ=255)
          tataAAATCAGTTATTTGATGAGTTGAAAACCTCTCTGGCAGAGGCTGTCAAGCTTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTACGGCTGATGTAGGCCGGAAAAGCCGTTCACGCCGTATCCGGCATTAGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGCGTCTCACAACTTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGc                                                                                                                                                                                                                         >  2:679700/1‑251 (MQ=255)
            aaaaaagaagttttttGAGGGGTTGAAAACCGCTCTGGAAGAGGTGGTCGAGATTAACAAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATGGCGGATATTAAGGCTTTTGTGGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGc                                                                                                                                                                                                                         <  2:410576/240‑1 (MQ=255)
                           ggatgatTTGAAACCCTCTCTGGAAGAGGTTGACGAGATAAAAAAAGGTTTAAAAAACCCTGCACGGGTGACCCGCTACGAAATGGCGGATACAAAGGTTGTTGTAGGCCGAAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                                                                                                                         <  1:619549/250‑1 (MQ=255)
                                    gAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAAGGTATTGAAATTTCTTGGTTTTTAGGCCGGATAAAGCGTTCACGCCGCATCCGGCATTTTAc                                                                                                                                                                                                >  2:461615/1‑251 (MQ=255)
                                         ccTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCAtt                                                                                                                                                                                           >  1:274487/1‑251 (MQ=255)
                                                            cgataatAAACAATGTTTGAAAAAACCGTCACGGGTGCCCCTTTACGAATTTGCGGATACCAAGGCTGTTGTAGCCCGGAAAGGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACACCGTATTGAAATTTCTTGGTTTGTAGGCCGTATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGt                                                                                                                                                                         <  1:379647/244‑1 (MQ=255)
                                                              agaTTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACt                                                                                                                                                                      >  2:33897/1‑251 (MQ=255)
                                                                     aCAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGAAATACGCTTGTACTTGATAcc                                                                                                                                                               >  1:600911/1‑250 (MQ=255)
                                                                       aagttttGAAAAACCCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTc                                                                                                                                                                                                                     <  1:250383/191‑1 (MQ=255)
                                                                                aaaaaCCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTcc                                                                                                                                                     >  1:672002/1‑250 (MQ=255)
                                                                                 aaaaCCAGCAAGGGTGACCCGCTACGAAATTGCGGATACTACGGCTGATGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTAAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCAAGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACACTTCGGCTcc                                                                                                                                                     >  2:449584/1‑249 (MQ=255)
                                                                                  aaaCCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAgt                                                                                                                                                  >  1:533766/1‑250 (MQ=255)
                                                                                     ccTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGtgcg                                                                                                                                               >  1:412526/1‑251 (MQ=255)
                                                                                           cGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCCTTTTACGCATTACTTCCCAATACGCTTGTACTTGATAAGCTTCGGCTCCAGTTCGTCTgcc                                                                                                                                         >  2:498114/1‑250 (MQ=255)
                                                                                                gACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTTCGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCCGTGCGTCTGTGCCCac                                                                                                                                    >  2:212432/1‑250 (MQ=255)
                                                                                                gACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGccc                                                                                                                                      >  1:177115/1‑249 (MQ=255)
                                                                                                 acccgataCGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCCCGACGTATCCGGCATTCGGTGCTCAATGCCTGATGCGCCGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTCTTGAAATTTCTTGGTTTGTATGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTCCTTCGCAATACGCTTGTACTTTATACGCTTCGGCTCCCGTGCGTCTGCGCCCAgc                                                                                                                                   >  2:675468/8‑251 (MQ=255)
                                                                                                        aCGAAATTGCGGATACTAAGGCTTTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAgcgtgcg                                                                                                                              <  1:33897/249‑1 (MQ=255)
                                                                                                         cgcaaTTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCCGCTCCAGTGCGTCTGCGCCCAGCGTGGGttt                                                                                                                           >  2:589034/4‑251 (MQ=255)
                                                                                                                        tAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGt                                                                                                             >  1:568799/1‑250 (MQ=255)
                                                                                                                         aaGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGtcc                                                                                                           >  2:113171/1‑249 (MQ=255)
                                                                                                                                       ccgaaaaaggcgtTCACGCCGTTTCCGGCATCCGGTGCTAAATGCCTGATGCGACGTTTACCGCGGCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTAcc                                                                                             <  2:572039/240‑1 (MQ=255)
                                                                                                                                         ggAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCtt                                                                                           >  1:115462/1‑251 (MQ=255)
                                                                                                                                                              gccGGCATTCTGTGCTCAATGCCTGATGCGACGCTTACCTCGTCTTATCAGGCCTACAGTTTTACAACTTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCGGTAAAGTTACCTTCGAAGAACTCAACTTTAcc                                                                        <  2:600911/248‑1 (MQ=255)
                                                                                                                                                                    aTTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGAGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCt                                                                >  2:380122/1‑251 (MQ=255)
                                                                                                                                                                      gcGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTg                                                               <  1:618194/249‑1 (MQ=255)
                                                                                                                                                                                aTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAgga                                                    <  1:559680/251‑3 (MQ=255)
                                                                                                                                                                                     tGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGc                                               >  1:136786/1‑251 (MQ=255)
                                                                                                                                                                                      gATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCg                                              >  1:182525/1‑251 (MQ=255)
                                                                                                                                                                                                          ttATCAGGCCTGGAGATTCATAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCCCGCCGCATCCGGCATTTTTCGCATTACCTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTAACTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGaa                           <  2:177115/250‑1 (MQ=255)
                                                                                                                                                                                                              caggcctggagattcaCAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAAc                      >  2:265543/16‑251 (MQ=255)
                                                                                                                                                                                                                ggcctggagattcaCAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACgg                    >  1:474189/14‑251 (MQ=255)
                                                                                                                                                                                                                     ggagattcaCAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTg               >  1:315263/9‑251 (MQ=255)
                                                                                                                                                                                                                              cAACGTATTGAAATTTCTTTGTTTGTAGGCCGGATAAGGCGTGCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTTTACTTGATACGCTTCGGCTCCAg                                                                                                                                                   <  2:101603/110‑1 (MQ=255)
                                                                                                                                                                                                                              cAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAg                                                                                                                                                   >  1:101603/1‑110 (MQ=37)
                                                                                                                                                                                                                               aaCGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGAATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAAcc     <  1:652233/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATc  >  1:581496/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          
AAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCGGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATCCAGAATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATC  >  NC_000913/4628565‑4629041

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