Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,628,289 G→T R325R (CGG→CGT)  nadR → DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,628,2890GT90.4% 129.2 / 2.4 52R325R (CGG→CGT) nadRDNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
Reads supporting (aligned to +/- strand):  ref base G (2/3);  new base T (25/22);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.38e-01

CCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGAT  >  NC_000913/4628051‑4628531
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  
cccTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTc                                                                                                                                                                                                                                        <  1:269281/251‑1 (MQ=255)
cccTGGTAAACAAACACGCCAATATCTTCAATACCACAAGGGCGTGGGCATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAATGGGGCGTTGATTGACACTGATTTTTTCACTACCCAGGCGTTCTGCAAAAAGTCCGAAGGGCGTGAGCACCCGTTc                                                                                                                                                                                                                                        <  2:326613/251‑1 (MQ=255)
 ccTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGGTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCg                                                                                                                                                                                                                                       <  2:354983/251‑1 (MQ=255)
   tGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACAATGATTTTGTCACTACCCAGGCGTTATGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGt                                                                                                                                                                                                                                      >  1:187155/1‑250 (MQ=255)
    gttAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGAGTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGt                                                                                                                                                                                                                                      <  1:41416/247‑1 (MQ=255)
      tAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTTCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGATTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGtaa                                                                                                                                                                                                                                    >  1:14485/1‑247 (MQ=255)
      tAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCCATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGAACGAAGGGCGTGGGCATCCGTTCGtac                                                                                                                                                                                                                                    >  1:160043/1‑247 (MQ=255)
              cTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGGATTCCGATTACGATAAAATCGCGCTTGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTTACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTtgat                                                                                                                                                                                                                          >  1:11482/1‑251 (MQ=255)
                 gCCAATATCTGCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGttgattga                                                                                                                                                                                                                       <  1:506858/251‑1 (MQ=255)
                   cAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCCTTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGGGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATtgatg                                                                                                                                                                                                                     >  1:269064/1‑251 (MQ=255)
                            cAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGtt                                                                                                                                                                                                            <  1:203025/251‑1 (MQ=255)
                             aataccaccAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCCGGAGTTCTGCAAAAAGAACGAAGGGCGGGAGCATCCGTTCGTACCGGCGTTGGTTGATGAATACCGttt                                                                                                                                                                                                           >  2:196808/4‑251 (MQ=255)
                                  caccaGGGCGTTGGAATATGTACGTGGTTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCTGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCAGAATACATTGATTGTGCAGTGAAATATGCCGATAAAGTGGCGTTTATTGACACTTATTTTTTCACTACCCAGGCGTGCTGCAAAAAGTACGAAGGGCGTTAGCATCTGTTCGTACAGGTGTTGATTGATGAAtttcgctgcgatc                                                                                                                                                                                                      <  2:307202/251‑13 (MQ=255)
                                        ggCGTGGGAATATGGTCGCGATTATTTCTTTTCTCACCCTGGCGGGGATGAGATCGCATTGCAGTATTCCGATTACTATAAAAACGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGa                                                                                                                                                                                                <  2:328977/250‑1 (MQ=255)
                                                   atGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTgg                                                                                                                                                                                       >  1:202623/1‑249 (MQ=255)
                                                     gggcgcgATTATGTCTTGTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGAGAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGAGCTGGTGATCCTGCTGGAGaa                                                                                                                                                                                   <  2:473612/248‑1 (MQ=255)
                                                             ttATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGTTTACGATAAAATCGCGCTGGGCCACGCCAAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACAATGATTTTGTCACTACCCAGGCGTTCTGCAAACAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGTTTGATGAATACCGGTTCGACATGGTGATCCTGCTGGAGAACAACACGc                                                                                                                                                                           >  2:109192/1‑251 (MQ=255)
                                                                  gcTTTTCCCACCGCGGCGGTTATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTTAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATGTATCGATACCGATTTTGGCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAAGACCGTTTCGATCTGGTGATCCTGCTGTATAACAACACGCCgtgg                                                                                                                                                                      <  2:181804/250‑1 (MQ=255)
                                                                     tttCACACCGCGGCGGTGATGAGATCGCATTGCAGTAATCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGAGCTGGTGATCCTGCTGGAGAACAACACGCCgtgggt                                                                                                                                                                    <  1:361968/250‑1 (MQ=255)
                                                                                       atgaGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGc                                                                                                                                                 >  1:483017/1‑251 (MQ=255)
                                                                                             tCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGAGGGTGGCGGATGGTTTACGCAGCCTCGGc                                                                                                                                           >  1:206325/1‑251 (MQ=255)
                                                                                              cGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGcg                                                                                                                                          >  1:190454/1‑250 (MQ=255)
                                                                                                      aGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTg                                                                                                                                  >  1:295725/1‑251 (MQ=255)
                                                                                                           tccgattACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATc                                                                                                                              <  1:196808/244‑1 (MQ=255)
                                                                                                              gattACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGAGCGTGAGCATCCGTTCGTACAGGCGTTGATTGAAGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGACGCCTCGGCAGTTCGGGGGAACGCa                                                                                                                           >  2:376689/4‑250 (MQ=255)
                                                                                                                                  gggCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATTCCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCt                                                                                                        <  1:335392/249‑1 (MQ=255)
                                                                                                                                  gggCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGGCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCtgg                                                                                                      <  1:507006/251‑1 (MQ=255)
                                                                                                                                   ggCCACGCTCCATACCTTGATTTTGCTGTGAAATATGCCAGTTAAGTGGCGTTTATTGACACTGATTTTGGCCCTACCCAGCCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACATGCGTTGATTGATTAATTCCGTTTCGATCTGGTGATCCTGCTGGAGAACCACACGCCGTGGGTGGCGTATGGTTTACGCAGCCTCGGCAGTTCGGTGTATCGCAAAGAGTTCCATAACTTGCtgc                                                                                                      <  2:281415/250‑2 (MQ=255)
                                                                                                                                     actagcacAATACATTTATATTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTGCGGTGGATCGCAAAGAGTTCCAGAACTTGCtggtgg                                                                                                   <  2:12135/244‑1 (MQ=255)
                                                                                                                                           acAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTCCCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGc                                                                                             <  1:68995/250‑1 (MQ=255)
                                                                                                                                           acAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGACCTTGCTGGTGGAGATGc                                                                                             >  1:84350/2‑251 (MQ=255)
                                                                                                                                            cagtACATTGATTTTGCAGTGAAATATGTCAATATCGTGGCGTTTATTGACACTGATTTTGTCAATACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCAGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGCTCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATTGTTTTAGCGGCCTCGGACGTTCGGTGGATCGCAACGAGTTCCAGAACTTGCTGATGGATATGcg                                                                                            >  2:42454/4‑250 (MQ=255)
                                                                                                                                              ctACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTAACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAAAACGCCTTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGGGGCTCGCAAAGAGTTCCAGAACTTGCtggtggggctgctgg                                                                                          >  2:283516/2‑242 (MQ=255)
                                                                                                                                                    tGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAag                                                                                       >  1:2255/1‑248 (MQ=255)
                                                                                                                                                                    tatGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGGGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAGGAGAACAATATCGAATTtgt                                                                    >  2:19146/1‑251 (MQ=255)
                                                                                                                                                                     atGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAAttcgt                                                                    <  1:451299/250‑4 (MQ=255)
                                                                                                                                                                                      tttattgacactgATTTTGTTACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGacgg                                                  >  1:108810/13‑248 (MQ=255)
                                                                                                                                                                                               actgATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGAc                                           >  1:28320/4‑249 (MQ=255)
                                                                                                                                                                                                ctgATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGATCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGCGAACAACACGCCGTGGGTGGCGGATGGTTTACGCTGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATTATGGAGGTGAACATTATCGAATTCGTGCGGGTTGCGGAGGACGATTATGACAgc                                        >  2:127061/3‑250 (MQ=255)
                                                                                                                                                                                                     tttGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGCTCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACCATATCGAATTCGTGAGGGTTGAAGGGGCCGATTATGACAgccggtt                                   >  2:493422/1‑245 (MQ=255)
                                                                                                                                                                                                       tGTCACTACCCAGGCGTTCAGCAAACAGAACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGCTGAATACCGTTTCGATCTGGTGCTCCTGCTGGAGAACAACACGCCGTGGGTGGAGGATGGTTTCCGCAGCCTCGGCAGTTCGGTGGCTCGCAAATAGTTCCAGAACTTGCTGGTGGCGGTGGTGGACGAGAACAATCCCGAATTCGTGCGGGTTGAAGCGGAAGCTTGTGACAGCCGTTTcc                                 >  2:242479/1‑251 (MQ=255)
                                                                                                                                                                                                          ccctaccctggCGTTCTACAAAAATTACGAAGGGCGTTAGCATCCGTTCGTACATGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACTCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTCTGACTGCCGTTTCctc                               <  2:23735/241‑2 (MQ=255)
                                                                                                                                                                                                          cactacccAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCctgc                              >  1:120610/8‑251 (MQ=255)
                                                                                                                                                                                                           actacccAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCctgcg                             <  2:11482/245‑1 (MQ=255)
                                                                                                                                                                                                                ccAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCctgcgctgcg                        <  2:2255/250‑1 (MQ=255)
                                                                                                                                                                                                                  aGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACCAGCCGTGGATGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGGTCGCAAAGAGGTCCAGAAATTGCTGGGGGGCATGCTGGAAGAGAACAATATCGAATTAGTGCGGGCTGAAGGGGACGATTATGGCAGCCGTTTCATGAGCTACGTg                      >  2:136761/1‑251 (MQ=255)
                                                                                                                                                                                                                  aGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTg                      <  1:418928/251‑1 (MQ=255)
                                                                                                                                                                                                                         cggCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATTGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTTCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCGGCGCTGCGTGGAGCTg                <  1:417382/248‑1 (MQ=255)
                                                                                                                                                                                                                           gCAAAAAGTACGAAGGGCGGGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGt              <  2:160043/250‑1 (MQ=255)
                                                                                                                                                                                                                             aaaaaGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCg           >  1:257814/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aCGAAGGGCGTGAGCATCCGTGCGTACAGGGGTTGATGGATGAATACCGTTTCGATCGGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCACGGCAGTTCGGTGGAGCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAgat     <  2:84350/250‑1 (MQ=255)
                                                                                                                                                                                                                                      gAAGGGCGTGAGCATCCGTTAGTACAGGCGTTGCTTGATGAATACCGTTTCGATCTGGTGATCATGCTGGAGATCACCACGCCGTGGGTGGCAGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTACAGAACTTGCTGGTGGAGATGCAGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGTTTATGACAGCCGTTTCCTGAGATGCGTGTAGCTGGTGCGGGAgatgat  >  2:259889/1‑251 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  
CCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGAT  >  NC_000913/4628051‑4628531

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