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.98.5% 300.0 / ‑1.7 67coding (215/993 nt)rbsRDNA‑binding transcriptional dual regulator RbsR
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (35/31);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.33e-01

ACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATGGCGACAGCGAT  >  NC_000913/3938198‑3938691
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                          
aCGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTAtt                                                                                                                                                                                                                                                    >  1:586263/1‑251 (MQ=255)
      ttttAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTCACCTCAAACTCAATCAAACACATACCATTGGACTGTTGATCACTGCCAGTACCCATCATTCTATTCAGAAc                                                                                                                                                                                                                                              >  2:449061/1‑251 (MQ=255)
           gACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTg                                                                                                                                                                                                                                         >  1:378636/1‑251 (MQ=255)
             cAGGCAGAGGTGACGCTTGACTACAATGAAAGTTGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGATTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATAACATTGGCATGTTGATCACTGCCAGTAACAATCTTTCTATTCAGAACTGGTGCg                                                                                                                                                                                                                                       >  2:500071/1‑251 (MQ=255)
             cAGGCAGAGGAGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCCGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCCTGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCg                                                                                                                                                                                                                                       >  2:500100/1‑251 (MQ=255)
                 cAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGc                                                                                                                                                                                                                                   >  2:269653/1‑251 (MQ=255)
                  agagGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAACCTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCg                                                                                                                                                                                                                                  >  2:155963/1‑251 (MQ=255)
                     ggTGACGCTTGGTTACAATGAAAGATGTTGCCCGCCTGGCGGGGGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCACAGTTGAAGCGGAGATTAAAGAACTCAATTACGCGCCATCAGCTTTGACGCGTAGCCTCAAACTCAATAAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCCGAACTGGTGCGTGGCGTTg                                                                                                                                                                                                                               >  2:229053/1‑251 (MQ=255)
                              tGGCTACAAGGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATGCAGAACTGGTGCGTGGTGTTGAACGCAgc                                                                                                                                                                                                                       <  2:577083/250‑1 (MQ=255)
                                   aCAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTTTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTCCCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGAGCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAgctgct                                                                                                                                                                                                                   <  2:506068/249‑1 (MQ=255)
                                   aCAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTc                                                                                                                                                                                                                 <  1:331396/251‑1 (MQ=255)
                                        gAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAAc                                                                                                                                                                                                             <  1:259766/250‑1 (MQ=255)
                                                          ggcgggcgTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGACCTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTcc                                                                                                                                                                                           <  1:398089/250‑1 (MQ=255)
                                                          ggcgggcgTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTGCTTCAGTGAAGCGATTTCCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTcc                                                                                                                                                                                           <  2:421401/250‑1 (MQ=255)
                                                              ggcgTTTCTACCTCAACAGTTTCTAACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGCTTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCCGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATCCCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGACCGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGc                                                                                                                                                                                      >  2:218144/1‑251 (MQ=255)
                                                               gcgTTTCTACCTCAACATGTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCa                                                                                                                                                                                     <  2:475964/251‑1 (MQ=255)
                                                                        ccTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGa                                                                                                                                                                              >  1:111911/1‑249 (MQ=255)
                                                                            aaCAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACAGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCg                                                                                                                                                                         >  1:56737/1‑250 (MQ=255)
                                                                                   tctcACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAAc                                                                                                                                                                   >  1:581903/1‑249 (MQ=255)
                                                                                   tctcACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACa                                                                                                                                                                  <  1:411049/250‑1 (MQ=255)
                                                                                            aTCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGCACAGCGGATGAAt                                                                                                                                                        >  2:317833/1‑251 (MQ=255)
                                                                                            aTCAATAAAGATCGCTTCGTCAGAGAAGCGCTTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAAt                                                                                                                                                        >  1:161881/1‑251 (MQ=255)
                                                                                             acAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCATCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATc                                                                                                                                                       <  2:468625/250‑1 (MQ=255)
                                                                                              cAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAATTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATGCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATc                                                                                                                                                       <  2:572222/250‑1 (MQ=255)
                                                                                                  aaaGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGCATCGCAAt                                                                                                                                                  >  1:499819/1‑251 (MQ=255)
                                                                                                     gATCGCTTCGTCAGTGAAGCGATTACCGCCACAGTTGAAGCGGCGATTACCGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCCGAACTGGTGCGTGGCGTTGCACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTg                                                                                                                                               >  2:198535/1‑251 (MQ=255)
                                                                                                     gATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGCACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTg                                                                                                                                               >  1:375925/1‑251 (MQ=255)
                                                                                                        cGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCGTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGa                                                                                                                                             <  1:243375/250‑1 (MQ=255)
                                                                                                         gCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGaa                                                                                                                                            >  1:244133/1‑250 (MQ=255)
                                                                                                         gCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCGGGaaa                                                                                                                                           <  1:494768/251‑1 (MQ=255)
                                                                                                            tCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAAGCTGTCTAGTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACg                                                                                                                                         <  2:440264/250‑1 (MQ=255)
                                                                                                                     aaGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGGACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCaa                                                                                                                                >  1:444431/1‑250 (MQ=255)
                                                                                                                       gCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCGCAAACACAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTAGTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCTTCCGGTGCAAGAGCGAAGTCGAAGAACAGCGGATGAATCGCAAGCTGGAAACGcgtatgtaaaaa                                                                                                                             <  2:552304/251‑12 (MQ=255)
                                                                                                                        cGATTACAGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGCAGCCTCCAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAAACGCGGATGAATCGCCATCTGGAAACGCTGATGCAAAAAc                                                                                                                            >  2:205695/1‑251 (MQ=255)
                                                                                                                         gATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAAcg                                                                                                                           >  2:206118/1‑251 (MQ=255)
                                                                                                                                gCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGAt                                                                                                                    >  1:385329/1‑251 (MQ=255)
                                                                                                                                 ccAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATg                                                                                                                   <  1:183058/251‑1 (MQ=255)
                                                                                                                                  cAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCACTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGCACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATgg                                                                                                                  >  1:38713/1‑251 (MQ=255)
                                                                                                                                  cAAAGTAGAAGCGGCGCTTAAAGACCTGAATTACGCTCCTTCTGCTCTAGAGCGTAACAACAAACTCCATCACACACATACCACTGGCATGATGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGAGAGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCACGGCCTTTGCAATACCGAAGGCGATGAAAACCGGGTGAATCGCAATCATGAAACGCTGCAGCAAAAACGCGATGATgg                                                                                                                  >  2:58779/1‑251 (MQ=255)
                                                                                                                                   aaaGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAACCACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGc                                                                                                                 >  1:301650/1‑251 (MQ=255)
                                                                                                                                    aaGTGGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATGCCAGTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATGCAGAACTGGTGCGTGGCGTGGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGGCCTGTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCt                                                                                                                <  2:575142/251‑1 (MQ=255)
                                                                                                                                          aaGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGaa                                                                                                                                                         >  1:259507/1‑204 (MQ=255)
                                                                                                                                          aaGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGaa                                                                                                                                                         <  2:259507/204‑1 (MQ=255)
                                                                                                                                          aaGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGt                                                                                                          >  2:263338/1‑251 (MQ=255)
                                                                                                                                            gcggcgATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGCTCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGCACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGACTCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGtt                                                                                                         >  2:468740/1‑250 (MQ=255)
                                                                                                                                                  aTTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTAtt                                                                                                                                                                                                                                                    >  2:65825/1‑105 (MQ=255)
                                                                                                                                                  aTTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTAtt                                                                                                                                                                                                                                                    <  1:65825/105‑1 (MQ=255)
                                                                                                                                                    tAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGACTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCa                                                                                                 <  1:75371/250‑1 (MQ=255)
                                                                                                                                                     aaaGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATTCCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATTAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCAcc                                                                                               <  1:535376/251‑1 (MQ=255)
                                                                                                                                                     aaaGAACTCAATTACGCGCCAACAGCTCTGGCGCGTCGCCTCAAACTCACTCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGCCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCAcc                                                                                               >  2:289846/1‑251 (MQ=255)
                                                                                                                                                                  aCGCGCCATCAGCTCGGGCGCGTAGCCTCAAACTCAATCAAACACATGCCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGGGGCGTTGAACGCCGCTGCTTCGAACGCGGTTATAGTCTCGTCCGTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAAc                                                                                  <  2:42564/251‑1 (MQ=255)
                                                                                                                                                                     cgcCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCtt                                                                               <  1:271693/251‑1 (MQ=255)
                                                                                                                                                                     cgcCAGCAGCCCTGTCGCGTAGCCTCAAACTCAATCAACCACATCCCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGGTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAAcag                                                                                <  2:397432/250‑3 (MQ=255)
                                                                                                                                                                      gcCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTc                                                                              >  1:91220/1‑251 (MQ=255)
                                                                                                                                                                         aTCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTcgcg                                                                           <  2:56737/251‑1 (MQ=255)
                                                                                                                                                                         aTCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTcgcg                                                                           >  2:511070/1‑251 (MQ=255)
                                                                                                                                                                         aTCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTcgcg                                                                           <  1:488448/251‑1 (MQ=255)
                                                                                                                                                                                ctGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATc                                                                    >  2:503532/1‑251 (MQ=255)
                                                                                                                                                                                  ggCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGAGCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCGTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCTTGAAATCAt                                                                  <  2:296521/251‑1 (MQ=255)
                                                                                                                                                                                       gTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCAAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAGGGCGATGAACAGCGGATGAATCGCAATCTGGGAACGCTGATGCAAAAACGCGTTGATGGCTTGCGGTTACTGTGCACCGAAACGCATCAACCTTCGCGAGAAATCATGCacc                                                             >  1:47436/1‑249 (MQ=255)
                                                                                                                                                                                         aGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACg                                                            <  1:459378/250‑1 (MQ=255)
                                                                                                                                                                                              cAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATcc                                                      <  1:138247/251‑1 (MQ=255)
                                                                                                                                                                                                       ccaaACACATACCATGGGCATGGTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGGTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGcc                                             <  2:135313/250‑1 (MQ=255)
                                                                                                                                                                                                         aaaCACATACCAGTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATGCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCGCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTa                                           <  2:246840/251‑1 (MQ=255)
                                                                                                                                                                                                                   cATTGGCATGTTGATCACTGCCAGTACCACTCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTgatgat                                 >  1:489705/1‑251 (MQ=255)
                                                                                                                                                                                                                       ggCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGa                              <  1:209671/250‑1 (MQ=255)
                                                                                                                                                                                                                                   ccTGCCAGTACCAATCTTTCTATGCAGAACGGGTGCGTGGCGTTGAACGCAGCTGCTACGAACGCGGTTATAGTCTCGTCCTTTGCAATCCCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTGCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTc                 <  2:201846/250‑1 (MQ=255)
                                                                                                                                                                                                                                                aatcTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATGGCGACAGc     >  1:138306/4‑251 (MQ=255)
                                                                                                                                                                                                                                                  tcTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCATCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGACACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATGGCGACAGCg    >  2:312503/2‑250 (MQ=255)
                                                                                                                                                                                                                                                  tcTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGCAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATGGCGACAGCg    >  1:364696/2‑250 (MQ=255)
                                                                                                                                                                                                                                                   cTTTCTCTTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGACCGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGCCAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATGGCGACCGCGAt  >  2:521453/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                          
ACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATGGCGACAGCGAT  >  NC_000913/3938198‑3938691

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: