Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,866,959 A→C I198S (ATC→AGC)  rpoS ← RNA polymerase sigma factor RpoS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,866,9590AC95.2% 118.4 / ‑5.3 42I198S (ATC→AGC) rpoSRNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base C (18/22);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

GAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTC  >  NC_000913/2866718‑2867200
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   
gAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCg                                                                                                                                                                                                                                            >  1:445313/1‑251 (MQ=255)
            cAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTGCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCACTTGGTTCAt                                                                                                                                                                                                                               <  1:68941/251‑1 (MQ=255)
                ccTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCACTTGGTTCATGGTc                                                                                                                                                                                                                           <  1:28202/250‑1 (MQ=255)
                     aCGCTGTTTTGCGTTCAGCACGAACACCCATTTGACGATGATCTGCTTCATCTCGTCATCTTGCGTGGAATCTTCCGGCCCGTTCTCTTTTTCATCGGCCAGGATGACCATCAACGCTTTTTCGGAATCACCACCCAGCGGGGAGTCTACCGAGGTAATGCGCTCGTTAAGCCGAAGCATATGGATGACGCCATCACCTGGCCTATCCAGTTACCCTGCGCTCTCTTCCGCACTGGGTTCCTGGTCCAGCt                                                                                                                                                                                                                      >  2:277918/1‑251 (MQ=255)
                      cGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCACTTGGTTCATGGTCCAGCtt                                                                                                                                                                                                                     <  1:86880/251‑1 (MQ=255)
                             tGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGAc                                                                                                                                                                                                              >  2:122182/1‑251 (MQ=255)
                                cGTTCAGCTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAAc                                                                                                                                                                                                           >  1:216184/1‑251 (MQ=255)
                                       ctagAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTATTTTTCATCGACCAGGATGTCCAGCAACGCTTTTTCGGAATCAAAACACAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGTCTGACGTCATCACAGGGTTTATCCCGTTCCTCTGAGCTCTCTTCCGCGCTTGATTTATGTTCCAGCTTAAGGCACAACTAACGTg                                                                                                                                                                                                    >  2:363351/4‑251 (MQ=255)
                                       cTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTg                                                                                                                                                                                                    >  1:2537/1‑251 (MQ=255)
                                             ccgCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTACAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGCAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCACTTGGATCATGGTCCTGCTTATGGGTCAACTAACGTGCGATTc                                                                                                                                                                                             >  2:334995/3‑251 (MQ=255)
                                             cAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTc                                                                                                                                                                                             >  1:207240/1‑251 (MQ=255)
                                               gCCATTTAACGATGCTCTGCTTCATATCGTAATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCCGCGGGGTGTCTCCCGAGGTAATGCGCTCGTTAAGACGAAACATACGGCTGACGCCATCAAATGGCTTATCCAGGTGCTCTACGATCTCTTCCGACATTGGTTCATGGTCCACCTTTTGGGACAACTCAAGTGCGGttagc                                                                                                                                                                                           >  2:175651/1‑248 (MQ=255)
                                                           tGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGtt                                                                                                                                                                                >  1:338493/1‑250 (MQ=255)
                                                               ctGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCt                                                                                                                                                                           >  1:383077/1‑251 (MQ=255)
                                                                     cATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCttt                                                                                                                                                                      <  2:5476/250‑1 (MQ=255)
                                                                     cATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTa                                                                                                                                                                     <  1:266158/251‑1 (MQ=255)
                                                                         gcGTCATCTTGCGTGGAATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCGCTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGAt                                                                                                                                                                 <  1:396268/250‑1 (MQ=255)
                                                                          cgccATCTTGTGTGGTATCTTCCGGACCGTTCTCTTTTCCATCGGCCAGGCTGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATTCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGAtg                                                                                                                                                                <  2:9985/248‑1 (MQ=255)
                                                                                 ttGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCg                                                                                                                                                         >  1:313687/1‑251 (MQ=255)
                                                                                    cGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCCGCATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGCCGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTAAGATGTGAATCGGc                                                                                                                                                       >  2:229405/1‑250 (MQ=255)
                                                                                     gTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCGCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCaa                                                                                                                                                     <  1:26547/251‑1 (MQ=255)
                                                                                        gTATCTTCCGGACCGTTCTCTTTTTCATAGGCCAGGATGTCCAGCAACGCTTTTACGGAATCACCCCCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACg                                                                                                                                                  <  2:86047/251‑1 (MQ=255)
                                                                                         tATCTTCCGGACCGTTCTCTTTTTCATAGGCCAGGGTGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGCCATCAACTTGCTCATCCAGTTGATCTGCGCTCTCTTCCGCACTATGATCACTGTCCAGATTCTGGGACAACTCCCTGTCGGTTAGCACGTAAAAGTTCAGTCCCTTTActctaggaatcggcacacg                                                                                                                                                  >  2:354318/1‑232 (MQ=255)
                                                                                                             ttttCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGtt                                                                                                                                <  1:91988/248‑1 (MQ=255)
                                                                                                                  aTCGGACAGGAAGACCAGCAAAGCTTTTTCGAAATCCCCCCCCAGCGGGATGTCTCCCGCGGTAATGCGCTCGTTAAGACGCAGCATACGGCTGACGTCATAACCTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGTCAACTCACGTGCGGTTCGCAGGCAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGAACGGGTTTGGTTCATacacg                                                                                                                        >  2:226453/1‑247 (MQ=255)
                                                                                                                   ccGGCCAGGTTGTCCAGCATCTCTTGTTCGGAATCACACCCCAGCGGGGTGTCTACCGAGGTAAGGCGCTCGTTACGACGTAGCATCCGGCTGACGTCATTAACTGGCTTATCCAGTTGCTCTGCGCTCGCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTTCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGc                                                                                                                       <  2:126232/250‑1 (MQ=255)
                                                                                                                          ggATGTCCAGCAACGCTTTTTCGGAATCACCTCCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCa                                                                                                                <  1:332245/251‑1 (MQ=255)
                                                                                                                             gttCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCGCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGAGGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATc                                                                                                             <  1:38155/249‑1 (MQ=255)
                                                                                                                              gTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATGATCGCCCGTTCAATCg                                                                                                            <  1:297008/251‑1 (MQ=255)
                                                                                                                                ccAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTc                                                                                                          <  1:122182/251‑1 (MQ=255)
                                                                                                                                    cAACGCTGTTTCGGAACCACCACCCAGCAGGGTGTCTTCCGAGGTAATTCGCTCGGTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCACTCTTCCGCACTTGTTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGTGTTTGGTTCATAATCGCCCGTTCAAGCGTCTGGc                                                                                                      <  1:88050/251‑1 (MQ=255)
                                                                                                                                           tcttCGGAATCACCACCCAGTGGGTGGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGTATACGGCTGACGGCATCAACTGGCTTATCTAGTTGCTCTGCGCTCTCTTCCTCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCTGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTcttccgtcttgcgaat                                                                                                  <  1:231226/246‑16 (MQ=255)
                                                                                                                                                     caccacCCAGCGGTGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGc                                                                                      <  1:447962/250‑1 (MQ=255)
                                                                                                                                                         ccccAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGGCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAAGCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATg                                                                                 <  1:346499/250‑1 (MQ=255)
                                                                                                                                                            cAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGAGAATCCACCAGGTTGCGTATGttt                                                                              >  2:225437/1‑250 (MQ=255)
                                                                                                                                                                          ccGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAAcc                                                                 >  1:389994/1‑250 (MQ=255)
                                                                                                                                                                                tAATGCGCTCGTTAAGACGAAGCATACGGCTGCCGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGTACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCCATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAAcccagttcc                                                          >  2:228654/1‑244 (MQ=255)
                                                                                                                                                                                               gACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCGCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGGATGTTTAGAAGCGGAAACCACGTGCCGGGTCAAACTTctct                                           <  1:353727/251‑1 (MQ=255)
                                                                                                                                                                                                   aaGCATACGGCTGAAGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACcg                                       >  1:142651/1‑251 (MQ=255)
                                                                                                                                                                                                          cGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGa                                 >  1:395211/1‑250 (MQ=255)
                                                                                                                                                                                                                  gTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCCCTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCa                        <  1:189593/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gCTCTGCGCTCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTc  <  1:330514/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   
GAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTC  >  NC_000913/2866718‑2867200

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