Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 89,471 T→G intergenic (+439/‑163) 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,4710TG96.1% 180.7 / ‑2.4 51intergenic (+439/‑163)cra/mraZDNA‑binding transcriptional dual regulator Cra/DNA‑binding transcriptional repressor MraZ
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base G (27/22);  total (29/22)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

CGATGAATTGCCGATTTAACAAACACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGTAACAAGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTTCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCAGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTATCGGGAACAGCTGCTTG  >  NC_000913/89241‑89715
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
agATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATa                                                                                                                                                                                                                                  <  1:661965/250‑1 (MQ=255)
           cGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTg                                                                                                                                                                                                                        >  2:608709/1‑250 (MQ=255)
             aTTTAACAAGCACTTTTCTTTGATTCTGCGCAAACCCGCTGGCATCAAGCGCCACACAGTCTCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTAACAGCTATACATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGata                                                                                                                                                                                                                     <  2:639485/251‑1 (MQ=255)
                       cacTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTaa                                                                                                                                                                                                           <  1:233832/251‑1 (MQ=255)
                                  gCTTTTGCGCAAACCCGCTGGCATCAATCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGttt                                                                                                                                                                                                >  1:657176/1‑251 (MQ=255)
                                     tttGCGCAAACACGCTGGCATCAAGAGCCACTCAGACGTCACAAGGACTGTCAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGAGATGACGATGAGGGATTTTTTCTTACCGCTATTCATAACGTTACTTTGCTTCGCACGTTGGCCGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGtttttt                                                                                                                                                                                             >  2:126034/1‑251 (MQ=255)
                                         cgcAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGAAATTAGCCGTAAACATCGGGTTTTTTACCt                                                                                                                                                                                         >  2:553308/1‑251 (MQ=255)
                                              aCCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGt                                                                                                                                                                                     >  1:487235/1‑250 (MQ=255)
                                                ccGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATg                                                                                                                                                                                  >  1:144617/1‑251 (MQ=255)
                                                              cgcCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCt                                                                                                                                                                    <  1:183955/251‑1 (MQ=255)
                                                               gcCACACAGACGCAACAGGGACTGTTAACAGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCtt                                                                                                                                                                   >  1:257400/1‑251 (MQ=255)
                                                                     caGACGCAACAGGGACTGTTATCCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACGCGTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCTCACTTTGGACCTAAAATAAACAACGCGGATATTAGCCGTAAACATCGGGTTTTGTACCTCGGTATGCCTTTTGACTGGCTTGACAAg                                                                                                                                                             <  2:818669/251‑1 (MQ=255)
                                                                     caGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAg                                                                                                                                                             <  1:22558/251‑1 (MQ=255)
                                                                       gACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGCGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCACGTGATGCCTTGTTCCTGGCATGACCAGCt                                                                                                                                                           >  2:275766/1‑251 (MQ=255)
                                                                         cgcaaCAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCttt                                                                                                                                                         <  1:853206/248‑1 (MQ=255)
                                                                          gcaaCAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCtttt                                                                                                                                                        >  1:541498/3‑251 (MQ=255)
                                                                            aaCAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTcc                                                                                                                                                      >  1:635300/1‑251 (MQ=255)
                                                                               agggACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGTCTGGCTTGACAAGCTTTGCCTcc                                                                                                                                                   <  2:763115/249‑2 (MQ=255)
                                                                                  gACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTAGTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCTGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCAGCt                                                                                                                                                <  1:126034/251‑1 (MQ=255)
                                                                                   aCTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGAAAACTGAAACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGAACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCAGCTc                                                                                                                                               >  1:35562/1‑251 (MQ=255)
                                                                                                                 tGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTTGGGAAAATCGGGCGTCGTTTTAAAACGGGGAGGCCGATGAGGGATTTTTTCTTACAGCTATTAATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACATCGCTGAGATTAGCCGGATACATCGGGTTGTTTACCTCGGTATGCCTTGTGACGGGTTTGACAAGGTTTTACTACGCTCCGTAAAGTTCTTTCAGTgggcatatttggg                                                                                                                 >  2:196012/1‑241 (MQ=255)
                                                                                                                             cGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGTAACGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGGGGGAAATTGTGGGGCTAAGTGGGa                                                                                                      >  2:156699/1‑250 (MQ=255)
                                                                                                                                aaaCTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGAATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTAAGTGGGAAATTGGGGGGCAAAGTGGGAATa                                                                                                   >  2:592527/1‑250 (MQ=255)
                                                                                                                                   cTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTTCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCAGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAggg                                                                                               >  1:367809/1‑251 (MQ=255)
                                                                                                                                   cTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCCGTGGGAAATTGTGGGGCAAAGGGGGTATAAggg                                                                                               >  2:850827/1‑251 (MQ=255)
                                                                                                                                      ccccTTTATATTTGCTGTGGAAAATAGTGAGTCAGTTTAAAACTGTGAAGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAAGTTGCTGCGCCCGTTGGACGTAAAATAAACATCGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCACAGTGGGAATAAGGGGt                                                                                             <  2:64366/247‑1 (MQ=255)
                                                                                                                                            tataTTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCt                                                                                       <  2:385195/250‑1 (MQ=255)
                                                                                                                                                       tgGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAACCAACGCTGATATTAGCCGTCAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAACTTTTCCTCCGCCCCGTTAACTCCTTTCCGTGGGGAATTGTGGGGGAAAGTGGGAGTAAGGGGGGAGGATGGCATGTTCCgg                                                                           >  2:327949/1‑251 (MQ=255)
                                                                                                                                                       tgGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCgg                                                                           >  1:80902/1‑251 (MQ=255)
                                                                                                                                                         gAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATAGTAGCCGGAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAAGTGTGGGGCAAAGTGGGAATAAGGGGGGAGGCTggaaggggcagggg                                                                         <  1:829848/251‑13 (MQ=255)
                                                                                                                                                              tAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGCTTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCa                                                                     >  1:277350/1‑250 (MQ=255)
                                                                                                                                                               aGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTTACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATtgt                                                                                                                    >  2:584369/1‑202 (MQ=255)
                                                                                                                                                               aGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATtgt                                                                                                                    <  1:584369/202‑1 (MQ=255)
                                                                                                                                                                  gAGTCATTTTAAAACGGTGATGACGATGAGGGTTTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACCAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGTGAAATTGTGGGGGAAGGTGGGAATAAGGGGGGTGGCGGGGATGTTCCggggggacacgtt                                                                >  2:427254/1‑242 (MQ=255)
                                                                                                                                                                                cGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTCTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACa                                                  <  1:50818/251‑1 (MQ=255)
                                                                                                                                                                                  gTGATGATGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGc                                                >  1:620207/1‑251 (MQ=255)
                                                                                                                                                                                   tgatgaCGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGc                                                >  1:273766/1‑250 (MQ=255)
                                                                                                                                                                                     atgaCGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCaaa                                             <  1:156699/251‑1 (MQ=255)
                                                                                                                                                                                     atgaCGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCaa                                              >  1:288453/1‑250 (MQ=255)
                                                                                                                                                                                            tGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGGGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGACAATCTCGACAGCAAAGGGCGCt                                      >  1:91449/1‑251 (MQ=255)
                                                                                                                                                                                                  aTTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAg                                <  1:742098/251‑1 (MQ=255)
                                                                                                                                                                                                   ttttttCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTgg                                                                                                                            >  2:363487/1‑158 (MQ=255)
                                                                                                                                                                                                   ttttttCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTgg                                                                                                                            <  1:363487/158‑1 (MQ=255)
                                                                                                                                                                                                           tACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAATAAACAACGCTGATATTAGCCGTAACCATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATACGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCACACTCGACAGAAAAGGGCGCTTATCAGTGCCTACCCg                      >  2:689660/1‑251 (MQ=255)
                                                                                                                                                                                                               gCTATTCATAACGTTAATTTGCTGCGCAAGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGTCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCGTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTa                   <  2:169844/251‑1 (MQ=255)
                                                                                                                                                                                                                  aTTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTACAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTacc                 <  1:608709/250‑3 (MQ=255)
                                                                                                                                                                                                                     cATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGa             <  1:553308/251‑1 (MQ=255)
                                                                                                                                                                                                                       tAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAAc           >  1:437964/1‑251 (MQ=255)
                                                                                                                                                                                                                         aCGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTTATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAg         <  1:804025/251‑1 (MQ=255)
                                                                                                                                                                                                                                tttGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTg  >  1:237163/1‑251 (MQ=255)
                                                                                                                                                                                                                                 ttGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTg  <  1:327949/250‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
CGATGAATTGCCGATTTAACAAACACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGTAACAAGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTTCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCAGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTATCGGGAACAGCTGCTTG  >  NC_000913/89241‑89715

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