Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,867,454 T→A Q33L (CAG→CTG)  rpoS ← RNA polymerase sigma factor RpoS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,867,4540TA97.4% 131.4 / NA 38Q33L (CAG→CTG) rpoSRNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (16/21);  minor base C (1/0);  total (17/21)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTG  >  NC_000913/2867222‑2867684
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      
gCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCttt                                                                                                                                                                                                                       <  1:151657/250‑1 (MQ=255)
     aaCGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTc                                                                                                                                                                                                                  <  1:121686/250‑1 (MQ=255)
         gcggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCaaaa                                                                                                                                                                                                              >  2:88405/1‑250 (MQ=255)
                      acctccaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTtc                                                                                                                                                                                               <  2:79287/247‑1 (MQ=255)
                                          ctcGATCATCCGGCGGCGCCAGGCGCCCTCTCCACGCAGTGCGCGACGCGCCCAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGACTAACCAATCTCACCAATGTCAAGCTTAGTCGCGTCCACCACACTCTGTGTGTCTCCCTGCGATACCAGTTCCTCTTCGGCCCCATCGTTATCACTGTGTTCCCGTTCTACTCCGGCCTTTTCGTCAACAACCTCACCTCCGTTCTCTTCAACTTCCGCATCtt                                                                                                                                                                            >  2:51398/1‑251 (MQ=255)
                                           tcGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACACTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTc                                                                                                                                                                           >  1:292485/1‑251 (MQ=255)
                                            cGATCATCCGGCAACGGGGGGCGACATCGCCCCGCATTGCGCGACGCGCAAAATCAACTTCTTCTTAGGCAGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAACGCTGAGTCGCGTCCAACAAACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTGCGGCCAAATCGTTCTCCCTGAGTTCCAGTTCTGCGACGACCTTGTCGTCGAAAACCTCCACTCCGTTCTCATCAAATTCCGCATCTTc                                                                                                                                                                           >  2:146748/1‑250 (MQ=255)
                                             gATCATCCGGCGGCGAGAGGCGACATCTCCACGTAGTGCGCGACGCGCAAAATAAACTGCGTCTGCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCGGTGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTTGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAATCCTCAACTCCGTTCTCATCAAATTCCGCATCTTCAt                                                                                                                                                                         <  2:305797/251‑1 (MQ=255)
                                              atcatcCGGCGGCGAGAGGCGCCATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTCCAtt                                                                                                                                                                        >  2:10762/1‑251 (MQ=255)
                                              atcatcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCAtt                                                                                                                                                                        >  1:194031/1‑251 (MQ=255)
                                                       cggcgAGAGGCGCCATCTCCACGCAGTGCGCGACGCGCAAAATAACCTTCTGCGACGGCCGTTAACAGTGGTGAATAACCAATCTCACCAATGTAAAGCTGATTCGCTTCCAACACACGCTTTGTGGCTCCCTGCGATAACAGTTCCTCTTCTGCCAAATCGTTATCACGGGGTTCCACTTCTACTTAGGCCTTGTCTTCAAAAACCTCAACTCCTTTCTCAACAAATTCCGCATTTTCATTGAAATCATg                                                                                                                                                               <  2:199554/251‑1 (MQ=255)
                                                       cggcgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATg                                                                                                                                                               >  1:351051/1‑251 (MQ=255)
                                                            agagGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACtt                                                                                                                                                          >  1:313592/1‑251 (MQ=255)
                                                                   aCATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGt                                                                                                                                                   >  1:127844/1‑251 (MQ=255)
                                                                         ccACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACACTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTc                                                                                                                                               <  2:292485/249‑1 (MQ=255)
                                                                                     cgacgaGCAAAATAAACTTCTTCTTCGGCCGTCAACAGTGGCGAATAACCAATCTCACCCAGGTAAAACTGAGCCGCGTACAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCCTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACCTTCAGCGTATCCTGCCTCATAAtgtg                                                                                                                                 >  2:389887/1‑251 (MQ=255)
                                                                                            cAAAATAAGCTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCAGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAAGCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTa                                                                                                                          <  2:72293/251‑1 (MQ=255)
                                                                                                                      aaCAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCa                                                                                                <  1:327236/251‑1 (MQ=255)
                                                                                                                            ggtgAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATc                                                                                                                                                                  >  1:354873/1‑179 (MQ=255)
                                                                                                                            ggtgAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCaaa                                                                                          <  1:9953/251‑1 (MQ=255)
                                                                                                                            ggtgAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCaa                                                                                           <  2:277419/250‑1 (MQ=255)
                                                                                                                            ggtgAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCAGCCTCATTTAAATc                                                                                                                                                                  <  2:354873/179‑1 (MQ=255)
                                                                                                                               gAATAACCAATCTCACCAAGGTCAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCtt                                                                                                                                                                                                                        >  1:126951/1‑122 (MQ=255)
                                                                                                                               gAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCtt                                                                                                                                                                                                                        <  2:126951/122‑1 (MQ=255)
                                                                                                                               gAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGc                                                                                       <  1:242320/251‑1 (MQ=255)
                                                                                                                                      cAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTc                                                                                <  1:300444/251‑1 (MQ=255)
                                                                                                                                            cACCAAGGTTAAGCTGAGTCGCGTCCAACACCCTCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAAGCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAATCAAAAGCCTGTGTccgcc                                                                            <  2:360520/249‑1 (MQ=255)
                                                                                                                                                       aGCTGAGTCGCGTCCAACACACGCGGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCACCTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGATGGCTCCTACCCGTGATCCCTTGACGGCACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCtcc                                                               >  2:353561/1‑249 (MQ=255)
                                                                                                                                                             gTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCACGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCaa                                                          >  1:23826/1‑250 (MQ=255)
                                                                                                                                                                     cAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGc                                                   >  1:25528/1‑249 (MQ=255)
                                                                                                                                                                               ggtggggCTCCCTGCGATAACAGTTCCGCTTCGGCCAAAGCGGTATCACTGGGTTCAAGTTCTACTAAGGATGTTGCGTCAAAAACCGCAACTCCGTTCTCATCAAATTCCGCATCGTCATTTAAATCATGAACTTTCAGCGTAGTCGGACTCATAAGGGGGCGCCTAACCGTGATCCCTTGACGGAACATGCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCACATAACGCAGCGGGTTTACg                                       <  1:51398/246‑1 (MQ=255)
                                                                                                                                                                                    ggCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCa                                                                                                                                                                                >  1:227821/1‑109 (MQ=255)
                                                                                                                                                                                    ggCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAACAGCCTGGTTCCGCCGCTTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGAttt                                  >  1:105127/1‑251 (MQ=255)
                                                                                                                                                                                    ggCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTAGCGCTGCGGCAAATAACGCAGCGGGTTTACGGAttt                                  <  2:395706/251‑1 (MQ=255)
                                                                                                                                                                                    ggCTCCCTGCGAGAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCTTCAAAAACCTCAACTCCGCTCTCATCAAATTCCGCa                                                                                                                                                                                <  2:227821/109‑1 (MQ=255)
                                                                                                                                                                                            gCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCGGGTTCCGCCGATTTATCTCTGCGGCAAATAACGCAGCGGGTTGACGGATTTCCCCTTGt                          <  2:18742/251‑1 (MQ=255)
                                                                                                                                                                                                                 ccAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGc     <  2:203067/251‑1 (MQ=255)
                                                                                                                                                                                                                    aaTCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCGCAACTCCGTTCGCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTTGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGGAACGAATTTCAAAATGCAAGCgtg  <  2:46979/251‑1 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      
GCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTG  >  NC_000913/2867222‑2867684

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