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

AATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTGGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCC  >  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)
                                           cGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGGATGGCTGAAGCGCAAAATGATTCCCTGCTGCCg                                                                                                                                                                                                        >  1:397649/1‑251 (MQ=255)
                                                    aatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCg                                                                                                                                                                                               <  1:653138/251‑1 (MQ=255)
                                                    aatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTGGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTc                                                                                                                                                                                                <  1:318261/250‑1 (MQ=255)
                                                     atCAATGTGGCCTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGCGTGTGAAAAGGATGGCTGAAGCGCAAAATGATCACCCTCTGCCGGGATAccggt                                                                                                                                                                                              >  2:321355/1‑247 (MQ=255)
                                                       cAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGATAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGtt                                                                                                                                                                                             <  1:176268/250‑1 (MQ=255)
                                                       cAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGGCGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGtt                                                                                                                                                                                             <  1:468087/250‑1 (MQ=255)
                                                             gttctatgcTGCCGCGATGATAGTCTCTTTTGTTACGCTCTTTTGTCATGGCTTTGGTCCCGCGTTGTTACTGAATGCTTTAAATAATCGGGGTCACCGGTTTGGTGAGCGAGAAGAGCCAGGAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAAGTGGTTTCTTCTCTGAATTGCGGGAGTATGAAAAGTATGTCTGAAGCGCAAAATGATCCCCTGCTGTCGGGATACTCGATTAATGcc                                                                                                                                                                                      <  2:516151/243‑1 (MQ=255)
                                                                                    actcTTTTGTTACGCGGTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGGTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGcc                                                                                                                                                               <  1:434164/248‑1 (MQ=255)
                                                                                           tGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGGTACTCGTTTAATGCCCATCTGGGGGCGGGGTTAAACCCGATTggg                                                                                                                                                        >  1:46907/1‑249 (MQ=255)
                                                                                             ttACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGGGTATGAAAAGAATGGCTGGAGCGCCAAATGATTCCCTGCTGGCGGGGTACTCGTTTTATGCCCCTCTGGTGGCGGGttttgagccgcttgaggg                                                                                                                                                      >  2:98028/1‑236 (MQ=255)
                                                                                                cgcgTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGTTCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGGTTGAGGccc                                                                                                                                                    >  1:199882/1‑249 (MQ=255)
                                                                                                                    ggTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGATTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGAGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGCTTGTGGCCAACGGGTATTTCGAtttttttt                                                                                                                               >  1:624616/1‑250 (MQ=255)
                                                                                                                        acGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCATGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGACGATTGGGGCCAACGGTTTTCTCGATTTTTTTATcg                                                                                                                            >  1:24829/2‑250 (MQ=255)
                                                                                                                          gCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgacc                                                                                                                         <  2:521364/251‑1 (MQ=255)
                                                                                                                          gCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGGCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAAc                                                                                                                                                                  <  1:407069/210‑1 (MQ=255)
                                                                                                                           cTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTTACGCCGCTTGAGGCCCACGGTTTTTTCGATTTTTTTTTcgaccc                                                                                                                        >  1:389190/1‑250 (MQ=255)
                                                                                                                            tttGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGGGTATGGAAAGTATGGGTGAAGAGCCAAATGATCCCCTGCTGCCGGGATAATCGTTTAATGCCCATCTGGTGGCGGGTTTTACGCCCATTTAGGCCCAAGGGTATCTCGATTTTTTTTTcgcacgg                                                                                                                       >  1:697393/1‑246 (MQ=255)
                                                                                                                             ttGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGGTTGAGGCCCACGGGTATCTCGGTTTTTTTATcgaccgac                                                                                                                      >  1:604919/1‑251 (MQ=255)
                                                                                                                                   cAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTgg                                                                                                                <  1:267106/251‑1 (MQ=255)
                                                                                                                                       aTGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAAt                                                                                                            <  1:98028/251‑1 (MQ=255)
                                                                                                                                                ataaGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTTAAGGCGATTGGGGGCCACGGGTTTCTCGGTTTTTTTTTCGACCGACCGCTgggggatcggaaa                                                                                                     >  1:486778/1‑239 (MQ=255)
                                                                                                                                                aaaaGCGTGGTTACCTGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTggg                                                                                                               <  2:486778/237‑1 (MQ=255)
                                                                                                                                                      ggggTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGtt                                                                                                                                                                      <  2:227290/178‑1 (MQ=255)
                                                                                                                                                      ggggTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGtt                                                                                                                                                                      >  1:227290/1‑178 (MQ=255)
                                                                                                                                                         gttttcGGTTTGGTTAGCGGGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCACTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAAAGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCaa                                                                                       <  2:429565/246‑1 (MQ=255)
                                                                                                                                                            actggTTTGGTTAGAGAGAAGAGCCATTAAAAGACGCAGTTACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAAGGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGTCAACGGTTATCTCGATTTTTTATCGACCGACCGCTGTGAATGAAAGGTTATATTCTCAAtctc                                                                                   <  2:81429/248‑1 (MQ=255)
                                                                                                                                                            accgtttttgTTAGCGAGAATAGCCAGTAAAAGCCGCAGGGACGGTAATGTCTGATGCAATATGGACAATTGGTTTCGTCTCTGAATGGCGGGAGGATGAAAAGTATTGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTGTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAAtct                                                                                    <  2:383132/242‑1 (MQ=255)
                                                                                                                                                            aCCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAAtct                                                                                    <  2:199882/251‑1 (MQ=255)
                                                                                                                                                              cGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGGGGCCAACGGTTAt                                                                                                                                            >  1:189118/1‑196 (MQ=255)
                                                                                                                                                                     gTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGGTAATCGTTTAATGCCCATCTGGGGGCGGGGTTTACGCCGGTTTGGGCCAACGGTTCTCTCGATTTTTTTATTGACCGCACGCGGGGGATGAAAGGGTATTTTTTCAATCTCACCATTcg                                                                           >  2:192308/1‑251 (MQ=255)
                                                                                                                                                                                 agCCAGTTAAAGTCGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGggt                                                               <  2:397649/251‑1 (MQ=255)
                                                                                                                                                                                        agaaGACGCAGTGAAGGTAATGGCTGATGTAATAGGGACGATGTGTTTCTTCTCTTAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAATATGATCCCCTGCTGCCTGGATACTCGTTTAATGCCCATCTGGTGGCGGGGTTAACGCCGAGTGAGGCCAACGGTTATCGCGATTTTTTTATCGAGCGACTGCTGGGGAGGAAAGGTTATATTCGCAATCTCCCCATTCGCGGTCAGGGGGTGGTGaaa                                                       <  2:51452/249‑1 (MQ=255)
                                                                                                                                                                                        aaaaGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGTCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGaaa                                                       <  1:321355/251‑1 (MQ=255)
                                                                                                                                                                                                 aGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATg                                                                                                           >  1:373188/1‑194 (MQ=255)
                                                                                                                                                                                                 aGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATg                                                                                                           <  2:373188/194‑1 (MQ=255)
                                                                                                                                                                                                            tGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGGGGATGAAAAGTATGGCTGAAGCGCAAAAAGATCCCCTGCTGCCGGGATAATCGTTTAATGCACATTTGGTGGAGGGGTTAAAGCCGATTTAGGACCACGGGTCTCTCGAATTTTTTTTCGGCCGTCAGGTGGGAATGAAAAGGTTTATTCTCCAATTAACCAATACAGGTCAGGGGGTGGGTGATAATTCAGGACGAGAAttttt                                   >  1:86873/1‑249 (MQ=255)
                                                                                                                                                                                                              tCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGGGAATGGAAAGTATGGCTGCAGCGCAAAATCATCCCCTGCTGCCGGGATACTCGATTAATGCCCATCTGGTGGCGGGGTTAACGCCGATTGAGGCCACAGGTTTTCTCGTTTTTTTTATTGCACGCCCGCTGGGACTTAAACGTTATTTTCTCCATCTTCCCATTCTCGCGCAGGGGGGGGGGAAAAATCAGGGatcaggatttgttt                                 >  1:40661/1‑238 (MQ=255)
                                                                                                                                                                                                                 gATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGttt                                 <  1:192308/248‑3 (MQ=255)
                                                                                                                                                                                                                 gATGCAATATGGACAATTGGTTTCTTCTATGAATGGCGGGAGTATGAAAAGTATGGCTGCAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCCACGGTTATCTCCGTTTTTTTATCGACCGACCCCTGGGGATGAAAGGTTTTATTCTCAATCTCACCCTTCGCGGTCAGGGGGGGGGGAAAAATTAGGGACGAGAATTTGTTTGcc                              >  1:136549/1‑251 (MQ=255)
                                                                                                                                                                                                                             aCAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCggg                                                                                                                                                                        <  2:593818/105‑1 (MQ=255)
                                                                                                                                                                                                                             aCAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCggg                                                                                                                                                                        >  1:593818/1‑105 (MQ=255)
                                                                                                                                                                                                                                     tttCTTCTCTGAATGGCGGGAGTATGAAAAGGATGGCTGAAGCGCCAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAAAGCCGATTGAAGCCAAAGGTTATCTCGGTTTTTTTTTCGAACGACCGCTGGGAATGAAAGATTATATTCTTCAATTTACCAATTGCGGGCATGGGGTGGTGAAAAATCCGGGGGGAGAATTTTTTTTCCGACCGGGGGATATTTTGCTg          >  1:191411/1‑251 (MQ=255)
                                                                                                                                                                                                                                             ctGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTccc     <  2:604919/248‑1 (MQ=255)
                                                                                                                                                                                                                                               gAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGcc  <  1:11023/249‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       
AATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTGGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCC  >  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: