Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,492,583 (G)6→7 intergenic (‑58/+4) idnT ← / ← idnO L‑idonate/5‑ketogluconate/gluconate transporter/5‑keto‑D‑gluconate 5‑reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,492,5831.G89.2% 102.3 / 7.6 37intergenic (‑58/+4)idnT/idnOL‑idonate/5‑ketogluconate/gluconate transporter/5‑keto‑D‑gluconate 5‑reductase
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base G (16/17);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

TCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GG‑CTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAATACCGTTGACCTGAATATTGTGGCGCGCCAG  >  NC_000913/4492336‑4492826
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   
acATTGCCAGCCCGCCGGGCGTGCTGCCTATACCATTTTGTATAGAGTGCAGGACGTCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑Agggg‑gggct                                                                                                                                                                                                                                                   <  1:94633/250‑4 (MQ=255)
 ccttGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGCGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCACTAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑Agggg‑gggctg                                                                                                                                                                                                                                                  >  2:31014/3‑247 (MQ=255)
        aGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAAc                                                                                                                                                                                                                                            >  2:398114/1‑250 (MQ=255)
                      gctgaCGATACAATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTCAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGGAGG‑CTGTTAAACAGCCCCTAAAATgcc                                                                                                                                                                                                                             >  2:101645/1‑251 (MQ=255)
                                     ttGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCaacaaaca                                                                                                                                                                                                               <  1:379923/250‑1 (MQ=255)
                                       gTATAGAGTCCAGGACGGCCTGGGCCTCCACCCCTTCGGCAAACCCGACGACGGCACCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GG‑CTGTTAAACAGCCACTAACATGCCGCCATCAACAAacaaca                                                                                                                                                                                                            <  1:226794/250‑1 (MQ=255)
                                       gTATAGAGGGCAGGTCGCCTTGTGCATCCATCCCTTCGGTAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGTCGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GG‑CTGTTAAACAGCCACTAACATGCCGCCATCAACAAacaaca                                                                                                                                                                                                            <  1:237704/250‑1 (MQ=255)
                                         ataGAGTGCAGGACGGCCTGTGCATCCCTCCCTTCGGCAAAACCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAAAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAATAACAg                                                                                                                                                                                                           <  1:90376/250‑1 (MQ=255)
                                                                   catccCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCttt                                                                                                                                                                                >  1:252025/1‑251 (MQ=255)
                                                                   catccCTCCGGCAAATCCGCCGACGGCAGCTACCAGAACGAGGGCAATAAAGCCTTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCttt                                                                                                                                                                                <  1:247996/251‑1 (MQ=255)
                                                                           cggaaaaTCCGACGATGACAGCCACCAGAACGAGGAAAAAAAAGCCGTTAACTTAAAAGACGATCATGAGAAAAAGAAGCAGCGCGACGCCAGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGGGCCAGAAAAGCAGCGCTTGCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCAGCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGaaa                                                                                                                                                                        <  2:24854/247‑1 (MQ=255)
                                                                           cGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTgaagaa                                                                                                                                                                         >  2:380067/1‑250 (MQ=255)
                                                                            ggCAAAACCGACTACGGCAGCTACCAAAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGACCATCAGGATAAGAAGCAGCGCGCCGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACGAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAg                                                                                                                                                                       <  1:138147/251‑1 (MQ=255)
                                                                                         acgGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCgcacc                                                                                                                                                           >  2:449458/1‑248 (MQ=255)
                                                                                           gcaacctccccaaccGAGGGCATAAAATCCGTTACCTTTAAAGACGACCATCAGGATAAGAACCAGCCCGCCCCCTGCCGCAATAATTATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAACCAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCAcc                                                                                                                                                        <  2:252025/237‑1 (MQ=255)
                                                                                              aGCTACCAGAACGAGGGAAATAAAGCCGTTAACTTTAAAGACGATCATCAGGAAAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAAt                                                                                                                                                     <  1:174694/251‑1 (MQ=255)
                                                                                                   ccAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGtt                                                                                                                                                >  1:406007/1‑251 (MQ=255)
                                                                                                     aGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTcc                                                                                                                                              >  2:429212/1‑251 (MQ=255)
                                                                                                           gggacAATAATGCGGTTAACTTTAAAGACTACCACCAGGATAAGAAGCAGCGGGACGCCGGCCGCAATAAGGATTAATGGCATAGTCTTTATCTCTTATTAGTGCTAGAAAAGCAGCCCTTTGTCTCTGCTCGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCCATCAGTTCCTGCGGa                                                                                                                                        <  1:31014/247‑1 (MQ=255)
                                                                                                            gggCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGAt                                                                                                                                       >  1:263674/1‑251 (MQ=255)
                                                                                                                        cGTTAATTTTAAAGACGATCATAAGGATAAGAGGCAGCGCGACGCCTGCCGCAAAAATGATTAATGGCATAGTCTTTATCTCTTATTAGGGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTg                                                                                                                           <  2:263674/251‑1 (MQ=255)
                                                                                                                                ttAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGAGCTCCCCAGCGTGcggcgg                                                                                                                     <  1:177533/249‑1 (MQ=255)
                                                                                                                                 tAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCggggg                                                                                                                  >  1:56637/1‑251 (MQ=255)
                                                                                                                                     gacGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCg                                                                                                              >  1:474832/3‑251 (MQ=255)
                                                                                                                                      acGATCATCCGGATCAGCAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTCGGGGG‑GG‑CTGTTAAACAGCCACTAACATGCCGCAATCAACAACCAACAGGTGGCCGTTTACGCAATCAGAGGCTTTTGCAGAAAGGAACAACGCAGCACCAATTAGTTCCTGCGGATCTCCACAGCGTGCGGCGGGGGGCCGt                                                                                                             >  2:442721/2‑251 (MQ=255)
                                                                                                                                          ccatcaGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGc                                                                                                         <  1:429212/250‑1 (MQ=255)
                                                                                                                                              caGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGcacac                                                                                                     >  1:54762/1‑250 (MQ=255)
                                                                                                                                                           agcGCGACGCCTGTCGCAATAATGATTAATGGCCTAGTCTTTCTCACTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGCGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCTTGCGGCGGTGGTGCGTTTGCACACCCACGCGGTgt                                                                                          >  2:196562/1‑248 (MQ=255)
                                                                                                                                                                                atGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCCGGCGGTGAAGGCTTCAGCCTCCACCAGAGCt                                                                   >  1:297382/1‑251 (MQ=255)
                                                                                                                                                                                         ggCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCAtt                                                          <  2:474832/251‑1 (MQ=255)
                                                                                                                                                                                                    tacctctTATAAGTGCCATAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACAGGCCGCCATCAACAAACAACAGGTGGCCGTTTACg                                                                                                                                                                                              <  1:71338/105‑1 (MQ=255)
                                                                                                                                                                                                    tATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACg                                                                                                                                                                                              >  2:71338/1‑108 (MQ=255)
                                                                                                                                                                                                      tctcTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGGGAAGGCGTCGTCCTCCACCCATGCTTTTTTCAttttctttttgtaa                                              >  2:96225/1‑239 (MQ=255)
                                                                                                                                                                                                         cgtattaggtccAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGc                                          <  2:82719/241‑1 (MQ=255)
                                                                                                                                                                                                           tattaGTGCCGGTAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCTCAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGc                                          <  1:288669/249‑1 (MQ=255)
                                                                                                                                                                                                            attaGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCg                                       <  1:195413/251‑1 (MQ=255)
                                                                                                                                                                                                                   ccAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAAt                                <  1:398114/251‑1 (MQ=255)
                                                                                                                                                                                                                                gcTTTCGCGCTGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTACGTTTTCAAAACCAAGCGGTGAAGGCTTCGGCCTCCACCCCTGCTTTTTTAATTTTTGTTTTTAAATAGCCCCgcgccaaaccggtttcctg                    >  2:437605/1‑235 (MQ=255)
                                                                                                                                                                                                                                          tGCACGT‑AGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTACGGCGGGGGTCCGTTTGCACAACCAGGCGGTGGAGGGTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGTAATAGCCCGGCCCAATACCGTTGACCTGAATAtttttgg         >  1:147317/1‑246 (MQ=255)
                                                                                                                                                                                                                                                        gg‑CTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAATACCGTTGACCTGAATATTGTGGCGCGCCAg  <  1:398825/245‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   
TCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGG‑GG‑CTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAATACCGTTGACCTGAATATTGTGGCGCGCCAG  >  NC_000913/4492336‑4492826

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