Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 89,103 G→C intergenic (+71/‑531) cra → / → mraZ DNA‑binding transcriptional dual regulator Cra/DNA‑binding transcriptional repressor MraZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091389,1030GC94.4% 179.7 / ‑1.4 55intergenic (+71/‑531)cra/mraZDNA‑binding transcriptional dual regulator Cra/DNA‑binding transcriptional repressor MraZ
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base C (19/32);  total (21/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

GAACTGCTCGACTTCTTACAGTGTCCGGTGCTGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGCGATGAATTGCCGATTTAACAAACACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGTAACAAGGACTGTTAACCGGGGAAGATATGTCCT  >  NC_000913/88856‑89349
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
gAACTGCTCGACTTCTTACAGTGTCCGGTGCTGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCtt                                                                                                                                                                                                                                                     <  1:459199/251‑1 (MQ=255)
  acggCTCGACTTCTTACAGTGTCCGGTGCTGGAAGTGGCTAAACGTCACCGCGATGGCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCGATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCGCTGGCGCGTAGAGTACGGGAATGGACAGCAATATGCTTaa                                                                                                                                                                                                                                                   <  2:724662/248‑1 (MQ=255)
    tGCTCGACTTCTTACAGTGTCCGGTGCTGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTaa                                                                                                                                                                                                                                                   >  1:774132/1‑249 (MQ=255)
    tGCTCGACTTCTTACAGTGTCCGGAGCTGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTCCGGGAATGGACCTCAATATCCTTaa                                                                                                                                                                                                                                                   >  2:436915/1‑249 (MQ=255)
         gACTTCTTACAGGGTCCGGTGCTGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAAt                                                                                                                                                                                                                                            <  1:96537/251‑1 (MQ=255)
           cttcttACAGTGTCCGGTGCTGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATa                                                                                                                                                                                                                                           >  1:766889/1‑250 (MQ=255)
                tACAGTGTCCGGTGCTTGCAGTGGTTCAACGTCCCCGCGATTGCGCATGGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCTTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCATCCGTAGCTAAGCTGCGAACAAAAATACGCGCCAG‑TGTAATTTCACTCTGGCTCGTAGAGTACGGGAATTGACATCAATATCCTTAAAGTAAAGAAGACt                                                                                                                                                                                                                                      <  2:524778/250‑1 (MQ=255)
                tACAGTGTCCGGTGCTGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGa                                                                                                                                                                                                                                        <  1:594090/248‑1 (MQ=255)
                    gtgtCCGGTGCTGGCAGTGGCTCACCGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAATCCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACAGCAATATCCTTAAAGTAAATAAGACTATTcc                                                                                                                                                                                                                                 <  2:476617/251‑1 (MQ=255)
                          ggTGCTGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGCCTATTCCTgcctt                                                                                                                                                                                                                           >  1:346025/1‑247 (MQ=255)
                                ggCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTg                                                                                                                                                                                                                     <  1:436915/251‑1 (MQ=255)
                                   aGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGAt                                                                                                                                                                                                                   >  1:397238/1‑250 (MQ=255)
                                     tGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATa                                                                                                                                                                                                                  >  1:359102/1‑249 (MQ=255)
                                      ggCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATaagt                                                                                                                                                                                                               <  1:561376/251‑3 (MQ=255)
                                       gCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATaagtt                                                                                                                                                                                                              >  1:465036/1‑248 (MQ=255)
                                           aaCGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTCGAGTACGGGAATGGCCATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATacgttcttt                                                                                                                                                                                                          >  2:233832/1‑243 (MQ=255)
                                           aaCGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATaagttcttt                                                                                                                                                                                                          >  2:818029/1‑244 (MQ=255)
                                                cACCGCGATGTCGCAGAGCTGGTGCTGGAGATTGTCCTGGTAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAAc                                                                                                                                                                                                     <  2:343393/251‑1 (MQ=255)
                                                   cgcgATGTCGCAGGGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAAcccg                                                                                                                                                                                                  <  2:81791/251‑1 (MQ=255)
                                                      gATGTCGCAGAGTGTGGGCTTGAGATTGTCCTGGCAAGCATGGACGACCCGCGTAAGCCAAAACCTGGTTTAACTCGCATTAAACGTAAGCTCGATCGCGGCGGCGTGCTCAGCCGGAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATGTCCCTATGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTAGTGATAAGTTCTTTTAAAcccgccc                                                                                                                                                                                               <  2:827764/251‑1 (MQ=255)
                                                        tGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAACACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGGAGAGTACGGGAATGGACATCAAGATCCTTAAAGTAAATAAGACTATTACTGACAAATATTGaaaaggtctttaaaaccagacagg                                                                                                                                                                                             >  2:661965/1‑228 (MQ=255)
                                                                     gtgCTGGAGATTGTCCTGGCAATCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCAcc                                                                                                                                                                                <  2:263692/251‑1 (MQ=255)
                                                                            agaTTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAATCGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGaa                                                                                                                                                                         <  1:16643/251‑1 (MQ=255)
                                                                                   ccTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTcac                                                                                                                                                                   <  1:699800/250‑1 (MQ=255)
                                                                                      ggCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACaat                                                                                                                                                                <  1:862926/250‑1 (MQ=255)
                                                                                       gCAAGCCAGGACACACCGCGTAAGCAAAAACCTTGTTTAACGCGCATTAAACGTAATCTCTAACGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCAAACAAAAATACGCGCCATGTG‑AATTTCCCTCTTGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTTATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTTAAATTCACaataa                                                                                                                                                              <  2:375728/251‑1 (MQ=255)
                                                                                        cAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATACGACTATTCCTGACTATTATTGCTAAGTTCTTTTAAACCCGCCCGTTAATTAACTCCCCAGCTGAAATTCACattaat                                                                                                                                                             >  2:183955/1‑251 (MQ=255)
                                                                                         aaGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAAtt                                                                                                                                                            <  1:29695/251‑1 (MQ=255)
                                                                                                 aCGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCTTAGAGTTCGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGata                                                                                                                                                    <  2:848709/251‑1 (MQ=255)
                                                                                                    aaCCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAAGTCACAATAATTAAGTGATATCg                                                                                                                                                 <  2:289382/251‑1 (MQ=255)
                                                                                                      ccGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGAc                                                                                                                                               <  1:315910/251‑1 (MQ=255)
                                                                                                       cgcgTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑ACTTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTTAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTAACCAGCTGAAATTCCCAATAATTAAGTGATATCGaa                                                                                                                                               >  2:22558/1‑249 (MQ=255)
                                                                                                           tAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAgcgc                                                                                                                                          >  1:84650/1‑251 (MQ=255)
                                                                                                                  aaaCCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTg                                                                                                                                   <  1:850668/251‑1 (MQ=255)
                                                                                                                  aaaCCTGGTTTAACACGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTg                                                                                                                                   <  2:766889/251‑1 (MQ=255)
                                                                                                                     ccTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCTCGGCGTGCTCAGCCGTAGCTAAGCCGAGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGATATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATAGCGACATCGCGTTTTTGCat                                                                                                                                <  2:780517/251‑1 (MQ=255)
                                                                                                                        ggTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCatta                                                                                                                              <  1:332828/250‑1 (MQ=255)
                                                                                                                             aaCGGGCATTTAACGTAAACTCCATCGCCGCGGCGTGCTTAGCCGTAGCTAAGCCGCGAAAAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTGTTTGCATTATTTTGt                                                                                                                        <  2:346025/251‑1 (MQ=255)
                                                                                                                                          cGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCCTTATTTTGTTACATgcag                                                                                                               >  2:853206/1‑245 (MQ=255)
                                                                                                                                          cGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATgcagagat                                                                                                           <  1:453537/251‑7 (MQ=255)
                                                                                                                                                cgctATCGCCGCGGCGTGCTCAGCCGTAGCGAAGCCGCGAACAAAAAGACGCGCCAGGTT‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAAtt                                                                                                      <  2:667460/248‑1 (MQ=255)
                                                                                                                                                  ctATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGCGTACGGGAATGGACATCAATATACTTAAAGGAAATAAGCCTATTCCTGACTATTATTGTTACGTTCTTTTACACCCGCCCGTTAATTAACTCACCAGCTGAAATTTACAATACTTATGTGATATCGCCAGCGCGTTTTTGCATTATTTTGTTACATGCCAAGATGAATTGcc                                                                                                   >  2:865544/1‑251 (MQ=255)
                                                                                                                                                          gcggcgTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATc                                                                                                                                                  <  1:246892/196‑1 (MQ=255)
                                                                                                                                                          gcggcgTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATc                                                                                                                                                  >  2:246892/1‑196 (MQ=255)
                                                                                                                                                            ggcgTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTaaca                                                                                          >  1:169844/1‑250 (MQ=255)
                                                                                                                                                            ggcgTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGAG‑AATTTCCATCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGCTACGTTCTTTTAAACCCGCCCGTTAATTAAATCACCAGCTGTAATTCACAATAATTAAGTGATATCGAAAGCTCGTTTTTGCATTATTTTGTTACATACAGAGATGAATTgacgcttaaact                                                                                          >  2:118913/1‑239 (MQ=255)
                                                                                                                                                                    cAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATGTCCCTCTGGCGCGTAGAGTACGTGAATGGTCATCAATATCCTTAAAGTAAATAAGACTATTCCTGTCTATTATGGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGGGATATCGGCAGCGCGTTTTTGCAGTATGTTGTTACGTGCAGAGATGAATTGCCGATTGAACAAGCACtt                                                                                   <  2:846675/249‑1 (MQ=255)
                                                                                                                                                                                    gcgAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCa                                                                  >  1:64366/1‑250 (MQ=255)
                                                                                                                                                                                      gAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAAc                                                               >  2:193337/1‑251 (MQ=255)
                                                                                                                                                                                                aCGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCa                                                      <  1:865540/250‑1 (MQ=255)
                                                                                                                                                                                                        ggTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGcc                                             <  1:786963/251‑1 (MQ=255)
                                                                                                                                                                                                                  cccTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGAc                                     <  1:603214/250‑1 (MQ=255)
                                                                                                                                                                                                                             gTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACt                         <  1:818029/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  ccatatCCTTAAAGTAAATAAGACTAGTCCTGACTATTATTGATATGTTCTTTTAAACCCGCCCGTTAAGTAACTCACCAGATGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCCCTGTTCTTTGCTTTTGCGCTAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTc    <  2:226494/249‑1 (MQ=255)
                                                                                                                                                                                                                                                    atatCCTTCAAGTAAATAAGACTATTCCTGCCTATTATTGATAAGTTCTTTTAAACCCGACCGTTAATTAAATCACCAGCTGCAATTAACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGAGATGAATAACCGATTTAACAAGCACTTTTCTTTCCTTTTGCGCAAACCCGCTAGCATCAAGCGCCACACAGACGCAACATGGACTGTTAACCGGGAAAGatatttgct  >  2:80093/1‑246 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
GAACTGCTCGACTTCTTACAGTGTCCGGTGCTGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTG‑AATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGCGATGAATTGCCGATTTAACAAACACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGTAACAAGGACTGTTAACCGGGGAAGATATGTCCT  >  NC_000913/88856‑89349

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