Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 70,371 T→C intergenic (‑323/‑16) araB ← / → araC ribulokinase/DNA‑binding transcriptional dual regulator AraC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091370,3710TC94.2% 151.1 / 0.6 53intergenic (‑323/‑16)araB/araCribulokinase/DNA‑binding transcriptional dual regulator AraC
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base C (23/26);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.58e-01

AATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTGGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAA‑GGTTATATTCT‑C‑AATCTCACCATTCGCGGTCAGGGGGTGG‑TGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCC  >  NC_000913/70126‑70613
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       
aaTCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTTCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACTCTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAg                                                                                                                                                                                                                                                    <  2:152459/251‑1 (MQ=255)
     tATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGt                                                                                                                                                                                                                                                   <  1:354895/247‑1 (MQ=255)
       gggATAAAAATGCTATGGCATAGCAAAGTGTGTCGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGaaa                                                                                                                                                                                                                                             <  2:421602/250‑1 (MQ=255)
                      tGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATgggttaaggg                                                                                                                                                                                                                              >  1:595927/1‑243 (MQ=255)
                      tGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAgcg                                                                                                                                                                                                                              >  1:521364/1‑251 (MQ=255)
                        gCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACTGGTTGGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCa                                                                                                                                                                                                                            <  1:41845/251‑1 (MQ=255)
                        gCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTATGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACTGGTTGGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCATTATGGACAATTGGTTTCTTCTCTGATTGTGGGGAGTATGAAAAGTATGGTTGAAgcccc                                                                                                                                                                                                                            >  2:41845/1‑248 (MQ=255)
                            aGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCaaaa                                                                                                                                                                                                                         <  1:475968/250‑1 (MQ=255)
                                aaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGTGGATGAAAAGTATGGCTGAAGCGCAAAATGATc                                                                                                                                                                                                                    >  2:109031/1‑251 (MQ=255)
                                           cGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGGATGGCTGAAGCGCAAAATGATT‑CCCTGCTGCCg                                                                                                                                                                                                        >  1:397649/1‑251 (MQ=255)
                                                    aatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCg                                                                                                                                                                                               <  1:653138/251‑1 (MQ=255)
                                                    aatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTGGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTc                                                                                                                                                                                                <  1:318261/250‑1 (MQ=255)
                                                     atCAATGTGGCCTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGCGTGTGAAAAGGATGGCTGAAGCGCAAAATGATCACCCT‑CTGCCGGGATAccggt                                                                                                                                                                                              >  2:321355/1‑247 (MQ=255)
                                                       cAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGATAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGtt                                                                                                                                                                                             <  1:176268/250‑1 (MQ=255)
                                                       cAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGGCGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGtt                                                                                                                                                                                             <  1:468087/250‑1 (MQ=255)
                                                             gttctatgcTGCCGCGATGATAGTCTCTTTTGTTACGCTCTTTTGTCATGGCTTTGGTCCCGCGTTGTTACTGAATGCTTTAAATAATCGGGGTCACCGGTTTGGTGAGCGAGAAGAGCCAGGAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAAGTGGTTTCTTCTCTGAATTGCGGGAGTATGAAAAGTATGTCTGAAGCGCAAAATGATC‑CCCTGCTGTCGGGATACTCGATTAATGcc                                                                                                                                                                                      <  2:516151/243‑1 (MQ=255)
                                                                                    actcTTTTGTTACGCGGTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGGTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGcc                                                                                                                                                               <  1:434164/248‑1 (MQ=255)
                                                                                           tGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGGTACTCGTTTAATGCCCATCTGGGGGCGGGGTTAAACCCGATTggg                                                                                                                                                        >  1:46907/1‑249 (MQ=255)
                                                                                             ttACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGGGTATGAAAAGAATGGCTGGAGCGCCAAATGATT‑CCCTGCTGGCGGGGTACTCGTTTTATGCCCCTCTGGTGGCGGGttttgagccgcttgaggg                                                                                                                                                      >  2:98028/1‑236 (MQ=255)
                                                                                                cgcgTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGTTC‑CCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGGTTGAGGccc                                                                                                                                                    >  1:199882/1‑249 (MQ=255)
                                                                                                                    ggTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGATTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGAGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGCTTGTGGCCAACGGGTATTTCGAtttttttt                                                                                                                               >  1:624616/1‑250 (MQ=255)
                                                                                                                        acGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCATGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGACGATTGGGGCCAACGGTTTTCTCGATTTTTTTATcg                                                                                                                            >  1:24829/2‑250 (MQ=255)
                                                                                                                          gCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgacc                                                                                                                         <  2:521364/251‑1 (MQ=255)
                                                                                                                          gCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGGCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAAc                                                                                                                                                                  <  1:407069/210‑1 (MQ=255)
                                                                                                                           cTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTTACGCCGCTTGAGGCCCACGGTTTTTTCGATTTTTTTTTcgaccc                                                                                                                        >  1:389190/1‑250 (MQ=255)
                                                                                                                            tttGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGGGTATGGAAAGTATGGGTGAAGAGCCAAATGATC‑CCCTGCTGCCGGGATAATCGTTTAATGCCCATCTGGTGGCGGGTTTTACGCCCATTTAGGCCCAAGGGTATCTCGATTTTTTTTTcgcacgg                                                                                                                       >  1:697393/1‑246 (MQ=255)
                                                                                                                             ttGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGGTTGAGGCCCACGGGTATCTCGGTTTTTTTATcgaccgac                                                                                                                      >  1:604919/1‑251 (MQ=255)
                                                                                                                                   cAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTgg                                                                                                                <  1:267106/251‑1 (MQ=255)
                                                                                                                                       aTGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAAt                                                                                                            <  1:98028/251‑1 (MQ=255)
                                                                                                                                                ataaGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTTAAGGCGATTGGGGGCCACGGGTTTCTCGGTTTTTTTTTCGACCGACCGCTgggggatcggaaa                                                                                                     >  1:486778/1‑239 (MQ=255)
                                                                                                                                                aaaaGCGTGGTTACCTGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTggg                                                                                                               <  2:486778/237‑1 (MQ=255)
                                                                                                                                                      ggggTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGtt                                                                                                                                                                      <  2:227290/178‑1 (MQ=255)
                                                                                                                                                      ggggTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGtt                                                                                                                                                                      >  1:227290/1‑178 (MQ=255)
                                                                                                                                                         gttttcGGTTTGGTTAGCGGGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCACTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAAAGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAA‑GGTTATATTCT‑C‑aa                                                                                       <  2:429565/246‑1 (MQ=255)
                                                                                                                                                            actggTTTGGTTAGAGAGAAGAGCCATTAAAAGACGCAGTTACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAAGGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGTCAACGGTTATCTCGATTTTTT‑ATCGACCGACCGCTGTGAATGAAA‑GGTTATATTCT‑C‑AAtctc                                                                                   <  2:81429/248‑1 (MQ=255)
                                                                                                                                                            accgtttttgTTAGCGAGAATAGCCAGTAAAAGCCGCAGGGACGGTAATGTCTGATGCAATATGGACAATTGGTTTCGTCTCTGAATGGCGGGAGGATGAAAAGTATTGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTGTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAA‑GGTTATATTCT‑C‑AAtct                                                                                    <  2:383132/242‑1 (MQ=255)
                                                                                                                                                            aCCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAA‑GGTTATATTCT‑C‑AAtct                                                                                    <  2:199882/251‑1 (MQ=255)
                                                                                                                                                              cGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGGGGCCAACGGTTAt                                                                                                                                            >  1:189118/1‑196 (MQ=255)
                                                                                                                                                                     gTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGGTAATCGTTTAATGCCCATCTGGGGGCGGGGTTTACGCCGGTTTGGGCCAACGGTTCTCTCGATTTTTTTATTGACCGCACGCGGGGGATGAAA‑GGGTATTTTTT‑C‑AATCTCACCATTcg                                                                           >  2:192308/1‑251 (MQ=255)
                                                                                                                                                                                 agCCAGTTAAAGTCGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAA‑GGTTATATTCT‑C‑AATCTCACCATTCGCGGTCAGGGggt                                                               <  2:397649/251‑1 (MQ=255)
                                                                                                                                                                                        agaaGACGCAGTGAAGGTAATGGCTGATGTAATAGGGACGATGTGTTTCTTCTCTTAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAATATGATC‑CCCTGCTGCCTGGATACTCGTTTAATGCCCATCTGGTGGCGGGGTTAACGCCGAGTGAGGCCAACGGTTATCGCGATTTTTTTATCGAGCGACTGCTGGGGAGGAAA‑GGTTATATTCG‑C‑AATCTCCCCATTCGCGGTCAGGGGGTGG‑TGaaa                                                       <  2:51452/249‑1 (MQ=255)
                                                                                                                                                                                        aaaaGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGTCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAA‑GGTTATATTCT‑C‑AATCTCACCATTCGCGGTCAGGGGGTGG‑TGaaa                                                       <  1:321355/251‑1 (MQ=255)
                                                                                                                                                                                                 aGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATg                                                                                                           >  1:373188/1‑194 (MQ=255)
                                                                                                                                                                                                 aGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATg                                                                                                           <  2:373188/194‑1 (MQ=255)
                                                                                                                                                                                                            tGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGGGGATGAAAAGTATGGCTGAAGCGCAAAAAGATC‑CCCTGCTGCCGGGATAATCGTTTAATGCACATTTGGTGGAGGGGTTAAAGCCGATTTAGGACCACGGGTCTCTCGAATTTTTTTTCGGCCGTCAGGTGGGAATGAAAAGGTT‑TATTCT‑CCAAT‑TAACCAATACAGGTCAGGGGGTGGGTGATAATTCAGG‑ACGAGAAttttt                                   >  1:86873/1‑249 (MQ=255)
                                                                                                                                                                                                              tCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGGGAATGGAAAGTATGGCTGCAGCGCAAAATCATC‑CCCTGCTGCCGGGATACTCGATTAATGCCCATCTGGTGGCGGGGTTAACGCCGATTGAGGCCACAGGTTTTCTCGTTTTTTTTATTGCACGCCCGCTGGGACTTAAA‑CGTTATTTTCT‑C‑CATCTTCCCATTCTCGCGCAGGGGGGGG‑GGAAAAATCAGGGatcaggatttgttt                                 >  1:40661/1‑238 (MQ=255)
                                                                                                                                                                                                                 gATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAA‑GGTTATATTCT‑C‑AATCTCACCATTCGCGGTCAGGGGGTGG‑TGAAAAATCAGGGACGAGAATTTGttt                                 <  1:192308/248‑3 (MQ=255)
                                                                                                                                                                                                                 gATGCAATATGGACAATTGGTTTCTTCTATGAATGGCGGGAGTATGAAAAGTATGGCTGCAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCCACGGTTATCTCCGTTTTTTTATCGACCGACCCCTGGGGATGAAA‑GGTTTTATTCT‑C‑AATCTCACCCTTCGCGGTCAGGGGGGGG‑GGAAAAATTAGGGACGAGAATTTGTTTGcc                              >  1:136549/1‑251 (MQ=255)
                                                                                                                                                                                                                             aCAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCggg                                                                                                                                                                        <  2:593818/105‑1 (MQ=255)
                                                                                                                                                                                                                             aCAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCggg                                                                                                                                                                        >  1:593818/1‑105 (MQ=255)
                                                                                                                                                                                                                                     tttCTTCTCTGAATGGCGGGAGTATGAAAAGGATGGCTGAAGCGCCAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAAAGCCGATTGAAGCCAAAGGTTATCTCGGTTTTTTTTTCGAACGACCGCTGGGAATGAAA‑GATTATATTCTTC‑AAT‑TTACCAATTGCGGGCATGGGGTGG‑TGAAAAATCCGGGGGGAGAATTTTTTTTCCGACCGGGGGATATTTTGCTg          >  1:191411/1‑251 (MQ=255)
                                                                                                                                                                                                                                             ctGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAA‑GGTTATATTCT‑C‑AATCTCACCATTCGCGGTCAGGGGGTGG‑TGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTccc     <  2:604919/248‑1 (MQ=255)
                                                                                                                                                                                                                                               gAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAA‑GGTTATATTCT‑C‑AATCTCACCATTCGCGGTCAGGGGGTGG‑TGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGcc  <  1:11023/249‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       
AATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTGGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATC‑CCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAA‑GGTTATATTCT‑C‑AATCTCACCATTCGCGGTCAGGGGGTGG‑TGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCC  >  NC_000913/70126‑70613

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: ‑