Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,492,470 C→A G19V (GGC→GTC)  idnT ← L‑idonate/5‑ketogluconate/gluconate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,492,4700CA92.5% 102.9 / ‑1.7 41G19V (GGC→GTC) idnTL‑idonate/5‑ketogluconate/gluconate transporter
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base A (15/22);  total (15/26)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

GCGCCCATTGCACGCGTTTTTTACCAAAAGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGG‑GG‑CTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGG  >  NC_000913/4492225‑4492714
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
gcgcCCATTGCCCGCGTTTTTTACCAAAAGTAGCAATCAGCGTAGTGGCGATACGTTGGGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGatc                                                                                                                                                                                                                                                   <  2:51382/251‑1 (MQ=255)
gcgcCCATTGCACGCGTTTTTTACCAAAAGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGatc                                                                                                                                                                                                                                                   <  2:76415/251‑1 (MQ=255)
                    ttACAAAAAGTAGAAATCAGCGTAGTGGCGATCCGTTGTGCCTCACCCGTATCATAAATCAGCTTGCCTATCATGGCCCCGAAACCGAGGATCATGGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTGCGGCAAATCCTACGACGGCAGCTACCAGAACGAGTGCAATAAAGCCGTTAACTTTAAAGTCGATCATCAGGATAAGAATCAgcg                                                                                                                                                                                                                                <  2:346492/250‑1 (MQ=255)
                      aCCAAAAGTAGCAATCAGCGTAGTGGGGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATGGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGTCCTGTGCAACCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGcgcg                                                                                                                                                                                                                              <  1:469032/250‑1 (MQ=255)
                           aaGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGACCCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGTATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTTTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCt                                                                                                                                                                                                                        <  1:417161/251‑1 (MQ=255)
                               aGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCg                                                                                                                                                                                                                    <  1:91601/251‑1 (MQ=255)
                                             tGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAAt                                                                                                                                                                                                      >  1:282802/1‑251 (MQ=255)
                                               gCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATgg                                                                                                                                                                                                    >  1:276729/1‑251 (MQ=255)
                                                  atccgtggtgccccccccGTATCAAAAACCCGCTTGCCTAACATGCCCCCAAAACCGAGGATCATGGCCAGCCCGCCGGGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGc                                                                                                                                                                                                   <  2:282802/234‑1 (MQ=255)
                                                                  ccctttCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCTATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCttat                                                                                                                                                                                 <  1:134784/246‑1 (MQ=255)
                                                                  ccGTATAAGAAATCAGCTTGCCTAACAGGGCCCGGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCttat                                                                                                                                                                                 <  1:272480/251‑1 (MQ=255)
                                                                        cAGAAACCAGCTTCCCAAACATGGCCTCGAAACCGGGGACCATTGCCAGCCCGCCGGGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCAATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTg                                                                                                                                                                            <  2:323952/250‑1 (MQ=255)
                                                                           aaaTCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAg                                                                                                                                                                        <  1:296595/251‑1 (MQ=255)
                                                                            aaTCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGcc                                                                                                                                                                          >  1:82719/1‑248 (MQ=255)
                                                                             aTCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGAc                                                                                                                                                                                                                            <  2:474614/197‑1 (MQ=255)
                                                                             aTCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGAc                                                                                                                                                                                                                            >  1:474614/1‑197 (MQ=255)
                                                                                                               acATTGCCAGCCCGCCGGGCGTGCTGCCTATACCATTTTGTATAGAGTGCAGGACGTCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAgggg‑gggct                                                                                                                                   <  1:94633/250‑4 (MQ=255)
                                                                                                                ccttGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGCGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCACTAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAgggg‑gggctg                                                                                                                                  >  2:31014/3‑247 (MQ=255)
                                                                                                                       aGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGG‑GGGCTGTTAAAc                                                                                                                            >  2:398114/1‑250 (MQ=255)
                                                                                                                                     gctgaCGATACAATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTCAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGAGG‑CTGTTAAACAGCCCCTAAAATgcc                                                                                                             >  2:101645/1‑251 (MQ=255)
                                                                                                                                                    ttGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCaacaaaca                                                                                               <  1:379923/250‑1 (MQ=255)
                                                                                                                                                      gTATAGAGTCCAGGACGGCCTGGGCCTCCACCCCTTCGGCAAACCCGACGACGGCACCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGG‑GG‑CTGTTAAACAGCCACTAACATGCCGCCATCAACAAacaaca                                                                                            <  1:226794/250‑1 (MQ=255)
                                                                                                                                                      gTATAGAGGGCAGGTCGCCTTGTGCATCCATCCCTTCGGTAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGTCGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGG‑GG‑CTGTTAAACAGCCACTAACATGCCGCCATCAACAAacaaca                                                                                            <  1:237704/250‑1 (MQ=255)
                                                                                                                                                        ataGAGTGCAGGACGGCCTGTGCATCCCTCCCTTCGGCAAAACCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAAAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAATAACAg                                                                                           <  1:90376/250‑1 (MQ=255)
                                                                                                                                                                                  catccCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCttt                                                                >  1:252025/1‑251 (MQ=255)
                                                                                                                                                                                  catccCTCCGGCAAATCCGCCGACGGCAGCTACCAGAACGAGGGCAATAAAGCCTTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCttt                                                                <  1:247996/251‑1 (MQ=255)
                                                                                                                                                                                          cggaaaaTCCGACGATGACAGCCACCAGAACGAGGAAAAAAAAGCCGTTAACTTAAAAGACGATCATGAGAAAAAGAAGCAGCGCGACGCCAGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGGGCCAGAAAAGCAGCGCTTGCGCGCTGCACGTAGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCAGCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGaaa                                                        <  2:24854/247‑1 (MQ=255)
                                                                                                                                                                                          cGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAag                                                                                                                                                                   >  1:162794/1‑145 (MQ=255)
                                                                                                                                                                                          cGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTgaagaa                                                         >  2:380067/1‑250 (MQ=255)
                                                                                                                                                                                          cGGCAAACCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAag                                                                                                                                                                   <  2:162794/145‑1 (MQ=255)
                                                                                                                                                                                           ggCAAAACCGACTACGGCAGCTACCAAAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGACCATCAGGATAAGAAGCAGCGCGCCGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGG‑GGGCTGTTAAACAGCCACGAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAg                                                       <  1:138147/251‑1 (MQ=255)
                                                                                                                                                                                                        acgGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCgcacc                                           >  2:449458/1‑248 (MQ=255)
                                                                                                                                                                                                          gcaacctccccaaccGAGGGCATAAAATCCGTTACCTTTAAAGACGACCATCAGGATAAGAACCAGCCCGCCCCCTGCCGCAATAATTATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGG‑GGGCTGTTAACCAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCAcc                                        <  2:252025/237‑1 (MQ=255)
                                                                                                                                                                                                             aGCTACCAGAACGAGGGAAATAAAGCCGTTAACTTTAAAGACGATCATCAGGAAAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAAt                                     <  1:174694/251‑1 (MQ=255)
                                                                                                                                                                                                                  ccAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGtt                                >  1:406007/1‑251 (MQ=255)
                                                                                                                                                                                                                    aGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTcc                              >  2:429212/1‑251 (MQ=255)
                                                                                                                                                                                                                          gggacAATAATGCGGTTAACTTTAAAGACTACCACCAGGATAAGAAGCAGCGGGACGCCGGCCGCAATAAGGATTAATGGCATAGTCTTTATCTCTTATTAGTGCTAGAAAAGCAGCCCTTTGTCTCTGCTCGTAGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCCATCAGTTCCTGCGGa                        <  1:31014/247‑1 (MQ=255)
                                                                                                                                                                                                                           gggCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGAt                       >  1:263674/1‑251 (MQ=255)
                                                                                                                                                                                                                                       cGTTAATTTTAAAGACGATCATAAGGATAAGAGGCAGCGCGACGCCTGCCGCAAAAATGATTAATGGCATAGTCTTTATCTCTTATTAGGGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTg           <  2:263674/251‑1 (MQ=255)
                                                                                                                                                                                                                                               ttAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGAGCTCCCCAGCGTGcggcgg     <  1:177533/249‑1 (MQ=255)
                                                                                                                                                                                                                                                tAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGG‑GGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCggggg  >  1:56637/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
GCGCCCATTGCACGCGTTTTTTACCAAAAGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGG‑GG‑CTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGG  >  NC_000913/4492225‑4492714

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