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.96.7% 255.3 / 1.2 62coding (215/993 nt)rbsRDNA‑binding transcriptional dual regulator RbsR
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (33/26);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATGGCGAC  >  NC_000913/3938198‑3938685
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
aCGCATTTTTAGAAAGGCAGAGGTGACGCTTGGATACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTAtt                                                                                                                                                                                                                                               <  2:263447/251‑1 (MQ=255)
    aTTTTTAGACAGGCAGCGGTGACGCTTGGCTACAATGAAAGATGTTGTCCGCCTGGCGGGCGGTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCCAAGTTGAAGAGGCGATTAAAGAACTCAATTACGATCACTCAGCTCTGGCGCGGAGCCTCAAACACAATCAAACACAAACCATTTGCATTTTGATCACTGCCAGTACCACTCTTTCTATTc                                                                                                                                                                                                                                              >  2:31927/1‑248 (MQ=255)
         tAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCGGTACCAATCTTTCTATTCAGAACTg                                                                                                                                                                                                                                       <  1:195977/250‑1 (MQ=255)
                             ttGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTAACTCAACAGTTTCTCACGTTATCAATAAAGACCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCACTCAACAACATACCATTGGCATGTTGATCACTGCCAGAACCAATCTTTCTATTCAGAACTGGTGCATGGCGTTGAACGctgc                                                                                                                                                                                                                  >  1:537013/1‑248 (MQ=255)
                              ggtatACAATGAAAGATGTTGCCCGCCTGGCGGGCGGTTCTACCTCAACAGTTTCTCACTTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACAGACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAgct                                                                                                                                                                                                                 <  2:461333/247‑1 (MQ=255)
                                   aCAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCACAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAgctgct                                                                                                                                                                                                              <  2:183335/249‑1 (MQ=255)
                                                 tGCCCGCCTGTCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCCGTGAAGAGATTATAGCCAAAGTTGAAGCGGCGACTAAAGAACTCCATTACGCTCCCTCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGACTGTTGAGTCCTGCCAGTACCCAACTTACTATTCAGATCTGCTGCGAGGCGTTGAACGCCACTACTTCTAACGCCGTtata                                                                                                                                                                                               >  2:216100/1‑250 (MQ=255)
                                                      gCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCAGGCGCGTAGCCTCCAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTTCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCg                                                                                                                                                                                         >  2:349447/1‑251 (MQ=255)
                                                            cgggcgTTACTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCttt                                                                                                                                                                                   >  2:468071/1‑251 (MQ=255)
                                                                 gTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCCCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAAt                                                                                                                                                                              >  1:474967/1‑251 (MQ=255)
                                                                           cAACAGTTTATCACGTTATCAAAAAAGATCTCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAATCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCGAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAggtg                                                                                                                                                                    <  2:413491/251‑3 (MQ=255)
                                                                                 tttCTCACGTTATCAATAAAGATCGCTTAGTCAGTGACACGATTACCGCCAAAGTTGAAGCGTCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAACCGCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGCACCAATCTTTCTATTCAGAACTAGTGCGTGGCGTTGACCACAGCTGCTTCGAAAGCGGATATAGTCTCGTCCTGTGCACTACCGAAGGCAATGaa                                                                                                                                                              >  1:491402/1‑250 (MQ=255)
                                                                                   tcccactataacAATAAAGATCGTATCGTCAGTGAAGTGATGTCCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACTTACCATTGGCATGTTGATCACTGCCAGTCCCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACa                                                                                                                                                            <  2:141169/239‑1 (MQ=255)
                                                                                      cACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCgg                                                                                                                                                        <  1:390420/251‑1 (MQ=255)
                                                                                        cGTTCTCAATAAAGTTCGCTGCGGCGGTGAGGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCCCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGACCGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGAt                                                                                                                                                      <  2:162661/251‑1 (MQ=255)
                                                                                            aTCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACCGCGGATGAAt                                                                                                                                                  >  2:365501/1‑251 (MQ=255)
                                                                                                   aaGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATc                                                                                                                                           >  1:315016/1‑251 (MQ=255)
                                                                                                     gATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGCTGAACAGCGGATGAATCGCAATCTg                                                                                                                                         >  1:142866/1‑251 (MQ=255)
                                                                                                      aTCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACACAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTAGCAATCTTTCTATTCAGAACTGGTGCGTGGCGGTGAACGCAGCTGCTTCGCACGCGGTTATAGTATCGTCATTTGCAATACCGAAGGCTGAGACCAGCGGCTTAATCGGAATCTgg                                                                                                                                        >  2:130035/1‑251 (MQ=255)
                                                                                                         gtttctccaatGAAGCGAGTACCGCCAAAGTTGAAGCGGCGATTAAATACCTCAAGTACGCGACATCAGCTCTGGCGCGGAGCCTCAAACCCAATCAAACACATACCATTGGCATGTTGATCTCTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGAGGAATCGCATTCTGGaaa                                                                                                                                     <  2:111751/241‑1 (MQ=255)
                                                                                                         gCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGAGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGaa                                                                                                                                      >  1:43857/1‑250 (MQ=255)
                                                                                                            tgggaaGTGAAGCGATGACCGAAAAAGGTGAAGCGGCGATGAAAGAACTCAATTACGCGCAATCAGCTCTGGCGCGTAGCCCCAAACTCAATCAAACACATTCCATTGGCATTTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACg                                                                                                                                   <  2:567917/245‑1 (MQ=255)
                                                                                                            tCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAAACAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAAc                                                                                                                                    <  1:14783/249‑1 (MQ=255)
                                                                                                              gTCAGTGAAGCGATTACCGCCAAAGTTCAAGCGGCGATTAAAGACCTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCACTCAAACACATACCATTGGCATGTTGCTCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGCACGCGGTTATAATCTCGTCCCTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGACACgatg                                                                                                                                >  1:506524/1‑248 (MQ=255)
                                                                                                                   ggAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCGTTGGCGCGTATCCTCAAACTCAATCAAACCCATACCATGGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCa                                                                                                                           <  2:128371/250‑1 (MQ=255)
                                                                                                                      aGCGATTACCGACAAAGGTGAAGCGGAGATTAAAGAACTTAATTACGCGCCATCAGCTCTGGCGCGTAGCCGCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGACGGCGATGAACAGCGGATGAATCGAAATCTTGAAAGGCttgagcaaa                                                                                                                         >  2:100164/1‑242 (MQ=255)
                                                                                                                             cccGCCAAAGATGAATGGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCGGGCCAGTAGCCTCAAACTCAATCAATCACACACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCGTTGCAATACCGAAGGCGATGAACCGCGGATGAATCGCAATCTTGAAACGCTGATGCAAAAACGCGtt                                                                                                                 <  2:243475/250‑1 (MQ=255)
                                                                                                                                cccAAATTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATTCCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGGCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGAt                                                                                                              <  1:529299/250‑1 (MQ=255)
                                                                                                                                   aaaGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGc                                                                                                           >  1:94079/1‑251 (MQ=255)
                                                                                                                                   aaaGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTACTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACAGAAGGCGATGAACAGCGGTTGAATCCCAATCTGGAAACGATGATGCAAAAACGCGTTGATggg                                                                                                           >  2:309209/1‑250 (MQ=255)
                                                                                                                                      gTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGAGTAGCCTCAAACTCAATCACTCACCTCCCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGCACGCAGCTGCATCGCACGCGGTTATAGTCTCGTCCTGTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAAAGCTTATACAAGACCGCGTTGTTGGATTg                                                                                                        >  2:5407/1‑251 (MQ=255)
                                                                                                                                              ggcgATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACt                                                                                                <  1:318781/251‑1 (MQ=255)
                                                                                                                                                 gATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACtgtg                                                                                             >  1:310222/1‑251 (MQ=255)
                                                                                                                                                    gaaaGAACCCATTGACGCTCCATCATCTCTGGCGCGTATCCTAAATCTCAATCAAACACATACCATTGGCATGTTGATCACGTCCAGTACCAATCTTTCTATGCAGAACGGGTTCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGGTATAGTCTCGTCCTTGGCAATACCGAAGGCGATGAACAGCGTATGAATCGCATGCGGGAAACGCTGAGGCAAAAACGCGTTGTTGGCTTGCTGTTACTGTGCa                                                                                           <  2:319532/249‑1 (MQ=255)
                                                                                                                                                      aagatctaaaTTCCGCGCCTTCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACTCATACCATTTGCATGTTGATCCCTGCCAGTACCAATCTTTCTATTCAGAACTGTTGCGTGGCGTTGAACTCTGCTGCGTCGACCGCGGTTATAGTCTCGTCCTTTGCAATACCTAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTTCACCg                                                                                        <  2:304973/243‑1 (MQ=255)
                                                                                                                                                      aaGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATAcc                                                                                                                                                                           >  2:488755/1‑169 (MQ=255)
                                                                                                                                                      aaGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATAcc                                                                                                                                                                           <  1:488755/169‑1 (MQ=255)
                                                                                                                                                          ccTCACTTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGCCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCTCAATCTGGAAACGCTGATGCCAAAACGCGTTGATGGCTTACTGTTACTGTGCACCGAAAc                                                                                    >  2:14402/2‑251 (MQ=255)
                                                                                                                                                               aTTACGCGCCATCAGCTCTGGCGCTTAGCCTCAAACTCAATCAAACCCATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATc                                                                               <  2:14884/251‑1 (MQ=255)
                                                                                                                                                                  aCGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATc                                                                                                                                                 >  2:335929/1‑182 (MQ=255)
                                                                                                                                                                  aCGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATc                                                                                                                                                 <  1:335929/182‑1 (MQ=255)
                                                                                                                                                                   cgcgCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGAt                                                                                                              >  1:361512/1‑216 (MQ=255)
                                                                                                                                                                   cgcgCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGCCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGAt                                                                                                              <  2:361512/216‑1 (MQ=255)
                                                                                                                                                                    gcgcCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCACTCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCt                                                                          >  1:522890/1‑251 (MQ=255)
                                                                                                                                                                    gcactattcgCTCTGGCGCGTAGCCCCAACCTCACTCAAATACATACCATTGGCATGTTGATCACCTCTAGTACCAACATATCTATTATGAACTGGTGCGTGGCGGTGAACGCAGCCGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCTATCCCGAAGGCGATGAACAGCGGATGAATCGCAATCTGTAAACGCGGATGCCAAAACGTGTTGATGGTGTGCTGTTACTGTGTACCGAAACGCATCAACCt                                                                          <  2:156534/242‑1 (MQ=255)
                                                                                                                                                                     cgcCATCAGCTCTGTCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAACGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGATACTGTGCACCGAAACGCATCAACCtt                                                                         >  2:395882/1‑251 (MQ=255)
                                                                                                                                                                        cATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTcgc                                                                      <  1:365501/251‑1 (MQ=255)
                                                                                                                                                                         aTCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGCTCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTcgcg                                                                     >  2:1106/1‑251 (MQ=255)
                                                                                                                                                                                     gcgTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGcc                                                         >  2:483103/1‑250 (MQ=255)
                                                                                                                                                                                         aGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGTAATCATGCAACg                                                      <  1:349447/250‑1 (MQ=255)
                                                                                                                                                                                                        caaACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGGGCACCGAAACGCATCAACCTTCGCGTGAAATAATGCAACGTTATCCGACAGTGccc                                      >  2:503178/1‑250 (MQ=255)
                                                                                                                                                                                                           acacAAACCAGTGGCATTTTGATCCCTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACt                                   <  2:503919/251‑1 (MQ=255)
                                                                                                                                                                                                              caTACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACtgtg                                >  1:120221/1‑251 (MQ=255)
                                                                                                                                                                                                                 aCCATTGGCATTTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGtgatg                             <  1:229457/251‑1 (MQ=255)
                                                                                                                                                                                                                    aTTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGAGCCTACTGtgatgctg                          >  2:403195/1‑248 (MQ=255)
                                                                                                                                                                                                                         cATGTTGATCACTGCCTGTACCAACCTTTCTATGCAGAACTGGTGCGGGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCATCGTTATCCGACAGTGCCTACTGTGATGATGGACTg                     <  1:301916/251‑1 (MQ=255)
                                                                                                                                                                                                                           tGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTggg                   >  1:507387/1‑251 (MQ=255)
                                                                                                                                                                                                                                   cACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGc                                        >  1:396856/1‑223 (MQ=255)
                                                                                                                                                                                                                                    aCTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGCCAGTGCCTACTGTGATGATGGACTGGGCTCCGTTc           >  1:559456/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aCTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTc           >  1:549654/1‑251 (MQ=255)
                                                                                                                                                                                                                                       gCCAGTACCAATCGTTCTAGTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGAt        <  2:94079/251‑1 (MQ=255)
                                                                                                                                                                                                                                              ccAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGGGCCTACTGTGATGATGGACTGGGCTCCGTTCGATGGCGAc  >  1:339480/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
ACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATGGCGAC  >  NC_000913/3938198‑3938685

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