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

GAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTC  >  NC_000913/2866718‑2867200
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   
gAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCg                                                                                                                                                                                                                                            >  1:445313/1‑251 (MQ=255)
            cAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTGCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCAC‑TTGGTTCAt                                                                                                                                                                                                                               <  1:68941/251‑1 (MQ=255)
                ccTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCAC‑TTGGTTCATGGTc                                                                                                                                                                                                                           <  1:28202/250‑1 (MQ=255)
                     aCGCTGTTTTGCGTTCAGCACGAACACCCATTTGACGATGATCTGCTTCATCTCGTCATCTTGCGTGGAATCTTCCGGCCCGTTCTCTTTTTCATCGGCCAGGATGACCATCAACGCTTTTTCGGAATCACCACCCAGCGGGGAGTCTACCGAGGTAATGCGCTCGTTAAGCCGAAGCATATGGATGACGCCATCACCTGGCCTATCCAGTTACCCTGCGCTCTCTTCCGCAC‑TGGGTTCCTGGTCCAGCt                                                                                                                                                                                                                      >  2:277918/1‑251 (MQ=255)
                      cGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCAC‑TTGGTTCATGGTCCAGCtt                                                                                                                                                                                                                     <  1:86880/251‑1 (MQ=255)
                             tGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGAc                                                                                                                                                                                                              >  2:122182/1‑251 (MQ=255)
                                cGTTCAGCTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGACAAc                                                                                                                                                                                                           >  1:216184/1‑251 (MQ=255)
                                       ctagAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTATTTTTCATCGACCAGGATGTCCAGCAACGCTTTTTCGGAATCAAAACACAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGTCTGACGTCATCACAGGGTTTATCCCGTTCCTCTGAGCTCTCTTCCGCGC‑TTGATTTATGTTCCAGCTTAAGGCACAACTAACGTg                                                                                                                                                                                                    >  2:363351/4‑251 (MQ=255)
                                       cTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTg                                                                                                                                                                                                    >  1:2537/1‑251 (MQ=255)
                                             ccgCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTACAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGCAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCAC‑TTGGATCATGGTCCTGCTTATGGGTCAACTAACGTG‑CGATTc                                                                                                                                                                                             >  2:334995/3‑251 (MQ=255)
                                             cAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTc                                                                                                                                                                                             >  1:207240/1‑251 (MQ=255)
                                               gCCATTTAACGATGCTCTGCTTCATATCGTAATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCCGCGGGGTGTCTCCCGAGGTAATGCGCTCGTTAAGACGAAACATACGGCTGACGCCATCAAATGGCTTATCCAGGTGCTCTACGATCTCTTCCG‑ACATTGGTTCATGGTCCACCTTTTGGGACAACTCAAGTG‑CGGttagc                                                                                                                                                                                           >  2:175651/1‑248 (MQ=255)
                                                           tGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTCGCAGGTAAACGtt                                                                                                                                                                                >  1:338493/1‑250 (MQ=255)
                                                               ctGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTCGCAGGTAAACGTTCAGCt                                                                                                                                                                           >  1:383077/1‑251 (MQ=255)
                                                                     cATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTCGCAGGTAAACGTTCAGCTCCttt                                                                                                                                                                      <  2:5476/250‑1 (MQ=255)
                                                                     cATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTCGCAGGTAAACGTTCAGCTCCTTTa                                                                                                                                                                     <  1:266158/251‑1 (MQ=255)
                                                                         gcGTCATCTTGCGTGGAATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCGC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTCGCAGGTAAACGTTCAGCTCCTTTACGAt                                                                                                                                                                 <  1:396268/250‑1 (MQ=255)
                                                                          cgccATCTTGTGTGGTATCTTCCGGACCGTTCTCTTTTCCATCGGCCAGGCTGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATTCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTCGCAGGTAAACGTTCAGCTCCTTTACGAtg                                                                                                                                                                <  2:9985/248‑1 (MQ=255)
                                                                                 ttGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCg                                                                                                                                                         >  1:313687/1‑251 (MQ=255)
                                                                                    cGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCCGCATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGCCGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTCGCAGGTAAACGTTCAGCTCCTTTAAGATGTGAATCGGc                                                                                                                                                       >  2:229405/1‑250 (MQ=255)
                                                                                     gTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCGCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCaa                                                                                                                                                     <  1:26547/251‑1 (MQ=255)
                                                                                        gTATCTTCCGGACCGTTCTCTTTTTCATAGGCCAGGATGTCCAGCAACGCTTTTACGGAATCACCCCCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACg                                                                                                                                                  <  2:86047/251‑1 (MQ=255)
                                                                                         tATCTTCCGGACCGTTCTCTTTTTCATAGGCCAGGGTGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGCCATCAACTTGCTCATCCAGTTGATCTGCGCTCTCTTCCGCAC‑TATGATCACTGTCCAGATTCTGGGACAACTC‑CCTGTCGGTTAGCACGTAAAAGTTCAGTCCCTTTActctaggaatcggcacacg                                                                                                                                                  >  2:354318/1‑232 (MQ=255)
                                                                                                             ttttCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGtt                                                                                                                                <  1:91988/248‑1 (MQ=255)
                                                                                                                  aTCGGACAGGAAGACCAGCAAAGCTTTTTCGAAATCCCCCCCCAGCGGGATGTCTCCCGCGGTAATGCGCTCGTTAAGACGCAGCATACGGCTGACGTCATAACCTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGTCAACTCACGTG‑CGGTTCGCAGGCAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGAACGGGTTTGGTTCATacacg                                                                                                                        >  2:226453/1‑247 (MQ=255)
                                                                                                                   ccGGCCAGGTTGTCCAGCATCTCTTGTTCGGAATCACACCCCAGCGGGGTGTCTACCGAGGTAAGGCGCTCGTTACGACGTAGCATCCGGCTGACGTCATTAACTGGCTTATCCAGTTGCTCTGCGCTCGCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTCGCAGGTAAACGTTCAGCTTCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGc                                                                                                                       <  2:126232/250‑1 (MQ=255)
                                                                                                                          ggATGTCCAGCAACGCTTTTTCGGAATCACCTCCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCa                                                                                                                <  1:332245/251‑1 (MQ=255)
                                                                                                                             gttCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCGCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTCGCAGGTAAACGTTCAGCTCCTTTACGAGGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATc                                                                                                             <  1:38155/249‑1 (MQ=255)
                                                                                                                              gTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATGATCGCCCGTTCAATCg                                                                                                            <  1:297008/251‑1 (MQ=255)
                                                                                                                                ccAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTc                                                                                                          <  1:122182/251‑1 (MQ=255)
                                                                                                                                    cAACGCTGTTTCGGAACCACCACCCAGCAGGGTGTCTTCCGAGGTAATTCGCTCGGTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCACTCTTCCGCAC‑TTGTTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGTGTTTGGTTCATAATCGCCCGTTCAAGCGTCTGGc                                                                                                      <  1:88050/251‑1 (MQ=255)
                                                                                                                                           tcttCGGAATCACCACCCAGTGGGTGGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGTATACGGCTGACGGCATCAACTGGCTTATCTAGTTGCTCTGCGCTCTCTTCCTCAC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CTGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTcttccgtcttgcgaat                                                                                                  <  1:231226/246‑16 (MQ=255)
                                                                                                                                                     caccacCCAGCGGTGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGc                                                                                      <  1:447962/250‑1 (MQ=255)
                                                                                                                                                         ccccAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGGCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAAGCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATg                                                                                 <  1:346499/250‑1 (MQ=255)
                                                                                                                                                            cAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGAGAATCCACCAGGTTGCGTATGttt                                                                              >  2:225437/1‑250 (MQ=255)
                                                                                                                                                                          ccGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAAcc                                                                 >  1:389994/1‑250 (MQ=255)
                                                                                                                                                                                tAATGCGCTCGTTAAGACGAAGCATACGGCTGCCGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGTACAACTCACGTG‑CGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCCATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAAcccagttcc                                                          >  2:228654/1‑244 (MQ=255)
                                                                                                                                                                                               gACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCGCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGGATGTTTAGAAGCGGAAACCACGTGCCGGGTCAAACTTctct                                           <  1:353727/251‑1 (MQ=255)
                                                                                                                                                                                                   aaGCATACGGCTGAAGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACcg                                       >  1:142651/1‑251 (MQ=255)
                                                                                                                                                                                                          cGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGa                                 >  1:395211/1‑250 (MQ=255)
                                                                                                                                                                                                                  gTCATCAACTGGCTTATCCAGTTGCTCTGCGCTCTCTTCCGCCC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCa                        <  1:189593/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gCTCTGCGCTCTCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTc  <  1:330514/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   
GAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCAC‑TTGGTTCATGGTCCAGCTTATGGGACAACTCACGTG‑CGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTC  >  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: ‑