Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,827,026 C→T A606V (GCG→GTG)  recG → ATP‑dependent DNA helicase RecG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,827,0260CT97.7% 136.0 / NA 44A606V (GCG→GTG) recGATP‑dependent DNA helicase RecG
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (17/26);  total (18/26)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTAAACAAGGTGAGCTACACCTGCTGGTTGCCACAACCGTTATTGAAGTCGGCGTTGATGTGCCTAACGCCAGTCTGATGATTATCGAAAACCCGGAGCGTCTGGGTCTGGCGCAGTTACACCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGCGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCTGGCGAATTGTTAGGCACGCGTCAGACGGGTAATGCTGAATTTAAAGTGGCGGATTTACTGCGCGATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTGATAGAACGCTGGATGCCGGAGACGG  >  NC_000913/3826780‑3827271
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
acGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTAAACAAGGTGAGCTACACCTGCTGGTTGCCACAACCGTTATTGAAGTCGGCGTTGATGTGCCTAACGCCAGTCTGATGATTATCGAAAACCCGGAGCGTCTGGGTCTGGCGCAGTTACACCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCa                                                                                                                                                                                                                                                    >  1:263024/1‑250 (MQ=255)
         tGAAACCTGCCGAGAAACAGGCGGAGATGGCGTCGTTTAAACAAGGGGAGCTACACCTGCTGGTTGCCACAACCGTTATTGAAGTCGGCGTTGATGTGCCTAACGCCAGTCTGATGATTATCGAAAACCCGGAGCGTCTGGGTCTGGCGCAGTTACACCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTg                                                                                                                                                                                                                                          >  1:63504/1‑251 (MQ=255)
         ggAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTAAACAAGGTGAGCTACACCTGCTGGTTGCCACAACCGTTATTGAAGTCGGCGTTGATGTGCCTAACGCCAGTCTGATGATTATCGAAAACCCGGAGCGTCTGGGTCTGGCGCAGTTACACCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCACTGCGTGCTGCGCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTg                                                                                                                                                                                                                                          <  2:24303/250‑1 (MQ=255)
                 gCCGAGAAACAGGCGGTGATGGCGTCGTTTAAACAAGGTGAGCTACACCTGCTGGTTGCCACAACCGTTATTGAAGTCGGCGTTGATGTGCCTAACGCCAGTCTGATGATTATCGAAAACCCGGAGCGTCTGGGTCTGGCGCAGTTACACCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGtgc                                                                                                                                                                                                                                   <  1:256575/250‑1 (MQ=255)
                  ccGAGAAACAGGCGGTGATGGCGTCGTTTAAACAAGGTGAGCTACACCTGCTGGTTACCACAACCGTTATTGAAGTCGGCGTTGATGTGCCTAACGCCAGTCTGATGATTATCGAAAACCCGGAGCGTCTGGGTCTGGCGCAGTTACACCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGtgctg                                                                                                                                                                                                                                 <  2:233942/251‑1 (MQ=255)
                   cGAGAAACAGGCGGTGATGGCGTCGTTTAAACAAGGTGAGCTACACCTGCTGGTTGCCACAACCGTTATTGAAGTCGGCGTTGATGTGCCTAACGCCAGTCTGATGATTATCGAAAACCCGGAGCGTCTGGGTCTGGCGCAGTTACACCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGtgctgc                                                                                                                                                                                                                                >  2:65071/1‑251 (MQ=255)
                                gTGAAGGCGGCGTTTAAACAAGGTGAGCTACACCTGCGTGTTGCCACAACCGTTATTGAAGTCGGCGTTGATGTGCCTAACGCCAGTCTGATGATTAACGAAAACCCGGGGCGTCTGGGTCTGGCGCAGTTACACCAGCTGCGGGGGCGCGTAGGTCGTTGCGCGGTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAacga                                                                                                                                                                                                                   <  2:301100/251‑1 (MQ=255)
                                 tGATGGCGTCGTTTAAACAAGGTGAGCTCCACCTGCTGGTTGCCACAACCGTTATTGAAGTCGGCGTTGATGTGCCTAACGCCAGTCTGATGATTATCGAAAACCCGGAGCGTCTGGGTCTGGCGCAGTTACACCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAacgac                                                                                                                                                                                                                  >  2:15748/1‑251 (MQ=255)
                                  gATGGCGTCGTTTAAACAAGGTGAGCTACACCTGCTGGTTGCCACAACCGTTATTGAAGTCGGCGTTGATGTGCCTAACGCCAGTCTGATGATTATCGAAAACCCGGAGCGTCTGGGTCTGGCGCAGTTACACCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAacgac                                                                                                                                                                                                                  <  1:140442/250‑1 (MQ=255)
                                  gATGGCGACGTTTAAACAAGGTGAGCTACACCTGCTGGTTGCCACAACCGTTATTGAAGGCGGCGTTGATGTGCCTAACGCCAGTCTGATGATTATCGAAAACCCGGAGCGTCTGGGTCTGGCGCAGTTACACCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGTTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAacgacg                                                                                                                                                                                                                 <  1:271013/251‑1 (MQ=255)
                                            tttAAACAAGGTGAGCTACACCTGCTGGTTGCCACAACCGTTTTTGAAGTCGGCGTTGATGTGCCTAACGCCAGTCTGATGATTATCGAAAACCCGGAGCGTCTGGGTCTGGCGCAGTTACACCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCCATGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGTTGCAAATTCGCCTGCAAGTGCGGCGCGACAGTAACGACGGTTTTGTGAt                                                                                                                                                                                                       >  2:228729/1‑251 (MQ=255)
                                            tttAAACAAGGTGAGCTACACCTGCTGGTTGCCACAACCGTTATTGAAGTCGGCGTTGATGTGCCTAACGCCAGTCTGATGATTATCGAAAACCCGGAGCGTCTGGGTCTGGCGCAGTTACACCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGAt                                                                                                                                                                                                       >  1:228228/1‑251 (MQ=255)
                                            tttAAACAAGGTGAGCTACACCTGATGGTTGCTACAACCGTTATTGAAGTCGGCGTTGATGTGCCTAACGCCAGTCTGATGATTATCGAAAACCCGGAGCGTCTGGGTCTGGCGCAATTACACCAGCTGCGAGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCACTGCGTGCTGCTATACAAAACGCCGCTTTCTAAAACGGCGCACATTCCCCTGCAACGGCTGCGCGACAGTAACGACGGTTTTGTGAt                                                                                                                                                                                                       >  2:228710/1‑251 (MQ=255)
                                              tAAACAAGGTGAGCTACACCTGCTGGTTGCCACAACCGTTATTGAAGTCGGCGTTGATGTGCCTAACGCCAGTCTGATGATTATCGAAAACCCGGAGCGTCTGGGTCTGGCGCAGTTACACCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTg                                                                                                                                                                                                     >  1:273071/1‑251 (MQ=255)
                                                       ggAGCTACACCTGCTGGTTGCCACAACCGTTATTGAAGTCGGCGTTGATGTGCCTAACGCCAGTCTGATGATTATCGAAAACCCGGAGCGTCTGGGTCTGGCGCAGTTACACCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTGCGCAAAAAg                                                                                                                                                                                            <  2:63731/250‑1 (MQ=255)
                                                                  tgctgGTTGCCACAACCGTTATTGAAGTCGGCGTTGATGTGCCTAACGCCAGTCTGATGATTATCGAAAACCCGGAGCGTCTGGGTCTGGCGCAGTTACACCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGGGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAACCGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTGAGCAAAAAGATCTGGAGAtt                                                                                                                                                                                 >  2:101455/1‑251 (MQ=255)
                                                                           ccACAACCGTTATTGAAGTCGGCGTTGATGTGCCTAACGCCAGTCTGATGATTATCGAAAACCCGGAGCGTCTGGGTCTGGCGCAGTTACACCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCt                                                                                                                                                                        >  2:48530/1‑251 (MQ=255)
                                                                                     tATTGAAGTCGGCGTTGATGTGCCTAACGCCAGTCTGATGATTATCGAAAACCCGGAGCGTCTGGGTCTGGCGCAGTTACACCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCTGGCGAATTg                                                                                                                                                               <  2:95489/250‑1 (MQ=255)
                                                                                              cGGCGTTGATGTGCCTAACGCCAGTCTGATGATTATCGAAAACCCGGAGCGTCTGGGTCTGGCGCAGTTACACCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCTGGCGAATTGTTAGGCAcgc                                                                                                                                                     >  2:281156/1‑251 (MQ=255)
                                                                                              cGGCGTTGATGTGCCTAACGCCAGTCTGATGATTATCGAAAACCCGGAGCGTCTGGGTCTGGCGCAGTTACACCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCACTGCGTGCTGCACTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCGGGCGAATTGTTAGGCAcg                                                                                                                                                      <  2:215908/250‑1 (MQ=255)
                                                                                               ggCGTTGATGTGCCTAACGCCAGTCTGATGATTATCGAAAACCCGGAGCGTCTGGGTCTGGCGCAGTTACACCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCTGGCGAATTGTTAGGCAcgcg                                                                                                                                                    >  1:104306/1‑251 (MQ=255)
                                                                                                    tGATGTGCCTAACGCCAGTCTGATGATTATCGAAAACCCGGAGCGTCTGGGTCTGGCGCAGTTACACCATCTGCGCGGGCGCGTAGGTCTTGGCGCGGTGGCTTCTCGCTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGGCAGTAACGACGGTTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCTGGCGAATTGTTAGGCACGCGTCAg                                                                                                                                                <  2:119324/250‑1 (MQ=255)
                                                                                                                           tgatTATCGAAAACCCGGAGCGTCTGGGTCTGGCGCAGTTACACCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCTGGCGAATTGTTAGGCACGCGTCAGACGGGTAATGCTGAATTTAAAGTg                                                                                                                        >  1:157386/1‑251 (MQ=255)
                                                                                                                               tATCGAAAACCCGGAGCGTCTGGGTCTGGCGCAGTTACACCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCTGGCGAATTGTTAGGCACGCGTCAGACGGGTAATGCTGAATTTAAAGTGGCg                                                                                                                     <  1:15748/250‑1 (MQ=255)
                                                                                                                                   gcaaaCCCGGAGCGTCTGGGTCTGCCGCAGATACACCAGATGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCCCACCGCGTGCTGCTCTACAAAACGCCGCTTTCTAAACCGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGCTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCTGGCGAATTGTTAGGCACGCGTCAGACGGGTAATGCTGAATTTAACGTGGCggcgct                                                                                                                >  2:191085/3‑247 (MQ=255)
                                                                                                                                   gAAAACCCGGAGCGTCTGGGTCTGGCGCAGTTACACCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCTGGCGAATTGTTAGGCACGCGTCAGACGGGTAATGCTGAATTTAAAGTGGCGGAttt                                                                                                                <  2:263024/251‑1 (MQ=255)
                                                                                                                                       aCCCGGAGCGTCTGGGTCTGGCGCAGTTACACCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCTGGCGAATTGTTAGGCACGCGTCAGACGGGTAATGCTGAATTTAAAGTGGCGGATTTACTg                                                                                                            >  1:233024/1‑251 (MQ=255)
                                                                                                                                               cGTCTGGGTCTGGCGCAGTTACACCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCTGGCGAATTGTTAGGCACGCGTCAGACGGGTAATGCTGAATTTAAAGTGGCGGATTTACTGCGCGATCa                                                                                                    <  1:228075/251‑1 (MQ=255)
                                                                                                                                                gtctgagtctggCGCGGTTACAACAGCTGCGCGGGAGCGAACGTCGTGGCGCGGTGGCTTCTCACTGCGTGATGCTCTAAAAAACGCCGCTTTCTAACACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCTGGCGAATTGTTAGGCACGCGTCAGACGGGTAATGCTGAATTTAAAGTGGCGGGTTTACTGCGCGAta                                                                                                     >  2:201738/1‑248 (MQ=255)
                                                                                                                                                    gggtctggCGCAGTTACACCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCTGGCGAATTGTTAGGCACGCGTCAGACGGGTAATGCTGAATTTAAAGTGGCGGATTTACTGCGCGATCAgg                                                                                                  <  2:228228/248‑1 (MQ=255)
                                                                                                                                                           gcgcAGTTACCCCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCCTTCTAAAACGGTGCAAAGTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCTGGCGAATTGTTAGGTACGCGTCAGACGGGTAATGCTGAATTTAAAGTGGCGGATTTACTGCGCGATCATGCGATGATccc                                                                                        <  2:275237/251‑1 (MQ=255)
                                                                                                                                                            cgcAGTTACACCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCTGGCGAATTGTTAGGCACGCGTCAGACGGGTAATGCTGAATTTAAAGGGGCGGATTTACTGCGCGATCAGGCGATGATCCCg                                                                                       <  2:198347/251‑1 (MQ=255)
                                                                                                                                                                   acacCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCTGGCGAATTGTTAGGCACGCGTCAGACGGGTAATGCTGAATTTAAAGTGGCGGATTTACTGCGCGATCAGGCGATGATCCCGGAAGttt                                                                                >  1:7820/1‑250 (MQ=255)
                                                                                                                                                                                         gttCGTGGCGCGGTGGATTCTCACTGCGTGCTGCTCTACAAAACGCCTCTTTCTAAAACGGTGCAAATTCGCTTGCAAGTGCTGCGCGACAGTAACGACGTTTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCTGGCGAATTGTTAGGCACGCGTTAGACGGGTAATGCTGAATTGAAAGTGGCGTATTGACTGCGCGATCAGGCGTGGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCa                                                          <  1:286140/249‑1 (MQ=255)
                                                                                                                                                                                                 cgcggtttcTTCTCACTGCGTGCTGCTCTACAAAACTCCGCTTTCTAAAAGGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGGGATTGCGCAAAAAGATCGGGAGACTCGCGGCCCTGGCGAATTGTTAGGCACGCGTCAGACGGGTAATGCTGTATTTAAAGTGGCGGATTTACTGCGCTATCAGGGGATGATGCCGGAGGTTGATCGCGCGGGACGCTATATTCTCGAACGtt                                                  <  2:193955/243‑1 (MQ=255)
                                                                                                                                                                                                  gcgGTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCTGGCGAATTGTTAGGCACGCGTCAGACGGGTAATGCTg                                                                                                                                   <  1:54209/169‑1 (MQ=255)
                                                                                                                                                                                                  gcgGTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCTGGCGAATTGTTAGGCACGCGTCAGACGGGTAATGCTg                                                                                                                                   >  2:54209/1‑169 (MQ=255)
                                                                                                                                                                                                        gCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCTGGCGAATTGTTAGGCACGCGTCAGACGGGTAATGCTGAATTTAAAGTGGCGGATTTACTGCGCGATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCcaca                                           <  1:57391/251‑1 (MQ=255)
                                                                                                                                                                                                        gCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCTGGCGAATTGTTAGGCACGCGTCAGACGGGTAATGCTGAATTTAAAGTGGCGGATTTACTGCGCGATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCcaca                                           <  2:66836/251‑1 (MQ=255)
                                                                                                                                                                                                                         gctCTACAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCTGGCGAATTGTTAGGCACGCGTCAGACGGGTAATGCTGAATTTAAAGTGGCGGATTTACTGCGCGATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCt                            <  1:6109/249‑1 (MQ=255)
                                                                                                                                                                                                                              aCAAAACGCCGCTTTCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCTGGCGAATTGTTAGGCACGCGTCAGACGGGTAATGCTGAATTTAAAGTGGCGGATTTACTGCGCGATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTGATAg                       <  2:302799/249‑1 (MQ=255)
                                                                                                                                                                                                                                           ttCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCTGGCGAATTGTTAGGCACGCGTCAGACGGGTAATGCTGAATTTAAAGTGGCGGATTTACTGCGCGATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTGATAGAACGCTGTATGcc          <  2:203321/249‑1 (MQ=255)
                                                                                                                                                                                                                                           ttCTAAAACGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCTGGCGAATTGTTAGGCACGCGTCAGACGGGTAATGCTGAATTTAAAGTGGCGGATTTACTGCGCGATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTGATAGAACGCTGGATGcc          <  2:85051/249‑1 (MQ=255)
                                                                                                                                                                                                                                                 aaCGGTGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCTGGCGAATTGTTAGGCACGCGTCAGACGGGTAATGCTGAATTTAAAGTGGCGGATTTACTGCGCGATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTGATAGAACGCTGGATGCCGGAGACgg  <  1:218536/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
ACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTAAACAAGGTGAGCTACACCTGCTGGTTGCCACAACCGTTATTGAAGTCGGCGTTGATGTGCCTAACGCCAGTCTGATGATTATCGAAAACCCGGAGCGTCTGGGTCTGGCGCAGTTACACCAGCTGCGCGGGCGCGTAGGTCGTGGCGCGGTGGCTTCTCACTGCGTGCTGCTCTACAAAACGCCGCTTTCTAAAACGGCGCAAATTCGCCTGCAAGTGCTGCGCGACAGTAACGACGGTTTTGTGATTGCGCAAAAAGATCTGGAGATTCGCGGCCCTGGCGAATTGTTAGGCACGCGTCAGACGGGTAATGCTGAATTTAAAGTGGCGGATTTACTGCGCGATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTGATAGAACGCTGGATGCCGGAGACGG  >  NC_000913/3826780‑3827271

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