Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,938,441 Δ1 bp coding (215/993 nt) rbsR → DNA‑binding transcriptional dual regulator RbsR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,938,4410C.97.8% 196.3 / ‑3.6 45coding (215/993 nt)rbsRDNA‑binding transcriptional dual regulator RbsR
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (27/17);  total (28/17)
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.52e-01

CGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTC  >  NC_000913/3938196‑3938676
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                           
cGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCt                                                                                                                                                                                                                                        >  1:54658/1‑250 (MQ=255)
                      agGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATGCAGAACGGGTGCGTGGCGtt                                                                                                                                                                                                                 <  2:150411/251‑1 (MQ=255)
                         tGACGCTTGCCTACAATGAAGAAGGTTGCCCGCCTGGCGGGAGTTTCTACCTCAACAGTTTATCACGTTATCAATCAAGATCGATTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGCTTAAAGAACTCAATTACGCGCCATCCGCTCTGGCGCGTAGCCTCAACCTCAATCCAACAAATACCATTGGCATGGTGATCACTGCCAGTAC‑AAAACGTTATATTCAGAACTGGTGCGTGGCGtcga                                                                                                                                                                                                               >  2:35836/1‑247 (MQ=255)
                           aCGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGTCGTTGAACg                                                                                                                                                                                                            <  1:343870/251‑1 (MQ=255)
                                      cAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACAGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACGCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTCCCAATC‑TTTCTATTCAGAACTGATGCGTGGCGTTGCACGCCGCTGCTTCg                                                                                                                                                                                                 >  1:422413/1‑251 (MQ=255)
                                       aaTGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCCGAACTGGTGCGTGGCGTTGAACGCCGCTGCTTCg                                                                                                                                                                                                 >  2:399029/1‑250 (MQ=255)
                                       aaTGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCg                                                                                                                                                                                                 >  1:14588/1‑250 (MQ=255)
                                           aaaGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTCCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAAcgc                                                                                                                                                                                            <  1:35940/251‑1 (MQ=255)
                                                 gTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCACTCAAACACATAACATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGGGGCGTTGAACGCAGCTGCTTCGAACGCGGTtat                                                                                                                                                                                      >  1:230194/1‑251 (MQ=255)
                                                          cTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCCATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTACTTCGAACGCGGTTATAGTCTCGTc                                                                                                                                                                             >  2:164187/1‑251 (MQ=255)
                                                           gggcgggcgTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGAGCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTcc                                                                                                                                                                            <  2:138501/250‑1 (MQ=255)
                                                              cgggcgTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATATTCTCGTCCttt                                                                                                                                                                         <  1:90144/251‑1 (MQ=255)
                                                               gggcgTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTg                                                                                                                                                                        >  1:392964/1‑251 (MQ=255)
                                                               gggcgTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTg                                                                                                                                                                        <  1:89829/251‑1 (MQ=255)
                                                                 gcgattaccaatcaccAGTTACTCACGTTATCAATAAAGATCGCTTCGTCCGTGAAGCGATTACCGCCACAGTTGAAGAGACGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCACACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTCCCAATC‑TTTCTATTCAGCACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGc                                                                                                                                                                       >  2:417222/16‑250 (MQ=255)
                                                                  cgTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCaa                                                                                                                                                                     <  2:14588/251‑1 (MQ=255)
                                                                             cAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCg                                                                                                                                                          >  2:279128/1‑251 (MQ=255)
                                                                                cAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATg                                                                                                                                                       <  1:39569/251‑1 (MQ=255)
                                                                                       tcACGTTATCAATAAAGATCGCTTCGTCAGTGCAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGCAGGCGATGAACAGCg                                                                                                                                                >  2:48286/1‑251 (MQ=255)
                                                                                                                gacAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTAAGCGCCATCAGATATAGCGCGGCGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTCTTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGAGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCtt                                                                                                                       >  2:240533/3‑250 (MQ=255)
                                                                                                                gTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTg                                                                                                                       >  1:434717/1‑251 (MQ=255)
                                                                                                                                 cGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTCGCGCGTAGCCTCAAACTCAATCAAACACATAACATTGGCAAGTTGATCACTGCAAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCACAATCTGGAAACGCTGCTGCAAAAACGCGTTg                                                                                                       >  1:24987/1‑250 (MQ=255)
                                                                                                                                  gCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGAt                                                                                                     <  2:158263/251‑1 (MQ=255)
                                                                                                                                       aGTTGAAACGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCACTCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGAGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATggatt                                                                                                >  2:265179/1‑248 (MQ=255)
                                                                                                                                         ttGAAGCGGCGATTACAGCACTAAATTACGCGCCCTCAGCTCTGGCGCGTAGCCTCAAACTAAATCAAACAAAAACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTATATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGACGGCGATGAACAGCGGATGACTCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGCTGGCTTGc                                                                                              >  2:392498/1‑251 (MQ=255)
                                                                                                                                                                          cATCAGCTCTGGCACGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGCACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAACCGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGCAACGCATCAACCTTcgc                                                             >  2:387666/1‑251 (MQ=255)
                                                                                                                                                                            tCAGCTCTGGCGCGTAGCCTAAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACTGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTACTGTTACTGTACACCGAAACGCATCAACCTTCGCGt                                                           >  1:299192/1‑251 (MQ=255)
                                                                                                                                                                             cAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTAGTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTTCTGTGCTCCGAAACGCATCAACCTTCGCGTg                                                          <  2:362002/251‑1 (MQ=255)
                                                                                                                                                                              aGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGa                                                                                                                              <  1:64102/182‑1 (MQ=255)
                                                                                                                                                                              aGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGa                                                                                                                              >  2:64102/1‑182 (MQ=255)
                                                                                                                                                                               gCTCTCGCGAGTAGCCTCACACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGCACTGGTACGTGGCGTTGACCACAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACAGAAGTCGATGAACAGCGGATGAATCGCAATCTGGAACCGCTGATGCAAAAACGCGTTGATGGCTTGCTGATACTGTGCAACGAAACGCATCAACCATCGCGTGa                                                         >  2:251232/1‑250 (MQ=255)
                                                                                                                                                                                  ctGGCGCGTAGCCTCAAACTCAATCAAACACATACCCTTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATc                                                     >  2:215888/1‑251 (MQ=255)
                                                                                                                                                                                      cgcgTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAGAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGc                                                 <  1:427908/251‑1 (MQ=255)
                                                                                                                                                                                      cgcgTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGc                                                 <  2:343905/251‑1 (MQ=255)
                                                                                                                                                                                         gTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCCGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAAc                                              >  2:201251/1‑251 (MQ=255)
                                                                                                                                                                                             ccTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGtt                                           >  2:95506/1‑250 (MQ=255)
                                                                                                                                                                                                cAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTAt                                         <  2:383268/249‑1 (MQ=255)
                                                                                                                                                                                                        atcaaACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGc                               >  1:424132/1‑251 (MQ=255)
                                                                                                                                                                                                         tcaaACACATACCATTGGCCTGTTGATCACTGCCAGTACCAATC‑TTCCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGcc                              >  1:165259/1‑251 (MQ=255)
                                                                                                                                                                                                          caaACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCt                             >  2:43071/1‑251 (MQ=255)
                                                                                                                                                                                                            aaCACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTAc                           >  2:256511/1‑251 (MQ=255)
                                                                                                                                                                                                                              gTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGc         >  2:210651/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGtt   <  2:137263/251‑1 (MQ=255)
                                                                                                                                                                                                                                    cACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCGCCGtt   <  2:212748/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aCTGCCAGTACCAATC‑TTTCTAGTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTc  <  1:169346/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                           
CGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTC  >  NC_000913/3938196‑3938676

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