Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,867,311 G→T R81S (CGT→AGT)  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,3110GT94.9% 151.0 / ‑3.6 63R81S (CGT→AGT) rpoSRNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base T (21/35);  total (24/39)
Fisher's exact test for biased strand distribution p-value = 5.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

GGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACGCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTG  >  NC_000913/2867066‑2867557
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       
ggTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑CCC‑CAGCGCGCG‑ACtcgcc                                                                                                                                                                                                                                                     >  1:641209/1‑249 (MQ=255)
     ggTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATaa                                                                                                                                                                                                                                               <  1:2549/251‑1 (MQ=255)
     ggTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATa                                                                                                                                                                                                                                                <  1:226887/250‑1 (MQ=255)
         cATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAActt                                                                                                                                                                                                                                           <  1:780957/251‑1 (MQ=255)
             aTCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAActtcttc                                                                                                                                                                                                                                       <  1:393566/251‑1 (MQ=255)
                 cccGTTCAATCGTCTGGCGAATCCACCAGGTTGCTTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑g                                                                                                                                                                                                                                  <  2:605619/251‑1 (MQ=255)
                     ttCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTAATCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGGGCGCG‑ACTCGCAAAAAAAACTTCTTCTTC‑GGCCg                                                                                                                                                                                                                              >  1:11731/1‑251 (MQ=255)
                         aTCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTTACACCAGACGCAAGTTACTCTCGATCCTCCGGCGGCGGGAGGGGACATCTCC‑AAG‑CCGTTCGCG‑GCTCGCAAAAAAAACATCTTCTTC‑GGACGTTa                                                                                                                                                                                                                           >  1:340150/1‑250 (MQ=255)
                                gtcGAATCCACCAGGTTGCGGATGTTGAGAAGCGGAAACCACTTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACTCCATACCACGATTTCCATAACGGCGGGCAATTTTTACCTCCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACCgtggt                                                                                                                                                                                                                   <  1:477941/249‑1 (MQ=255)
                                 gCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACGCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAgtggtg                                                                                                                                                                                                                  <  1:107141/251‑1 (MQ=255)
                                     aTCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATa                                                                                                                                                                                                              <  2:430009/251‑1 (MQ=255)
                                      aacaccaAGTTGCGTATGTAGAGAAGCGGAAACCACGTTCCAGGACAAAAATCTCAACCGCGCGGAGCAGCCCCAGGTTGCCCTCTTCGAAAAGGTCCAGCAACGCCAGACCACGATAGCCATAACGGCGGGCAATTTTTACCACCAGACGAAAGTTACTCTCGATCAACCGGCGGCGAGAGGCGAAAAATCC‑ACG‑AAGTGCGCG‑AAACACAAAAAAAACTTCTTCAGC‑GGCAATTAACAGTGGTGAATaa                                                                                                                                                                                                             <  2:400216/249‑1 (MQ=255)
                                           caggtgtcGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCACTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACTGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAAt                                                                                                                                                                                                        <  2:198233/244‑1 (MQ=255)
                                           caGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCGCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAAt                                                                                                                                                                                                        <  1:375936/251‑1 (MQ=255)
                                            aGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCCATTTTTTACACCCGACGCAAGTTAAACTCCGTCATCCCGCGGCGAGAGGCGACCTCTTC‑CCG‑CCATGCGGGCAC‑CGCAAAATAAAATTCTTTCTC‑GGCCGCTAACAGAGGAGCATAACCaaac                                                                                                                                                                                                       >  1:145040/1‑249 (MQ=255)
                                                tGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGGCATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCCGCGGCGGGAGGCGACATCCCC‑CCG‑CAGGGCGCG‑ACTCGCCAAATAAAATTCTTCTTC‑TGCCGGTTACAGTTGTGAATAACCAAACTCAc                                                                                                                                                                                                   >  1:9691/1‑251 (MQ=255)
                                                                cGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCCAGTTACTCTTGATTATCCGCCGGCGGGAGGCGGCCTCTCC‑CCC‑CCGGGCGAG‑ACTCGCAAAATAAACATTTTCTTT‑GGCCGTTTACAGGGGTGAATAACCAATCTCACCACGGGAAAGCTGAgg                                                                                                                                                                                   >  1:334002/1‑250 (MQ=255)
                                                                 ggAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCCTCCGGCGGCGAGAGGCGACATCTCC‑CCC‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑TGCCGTTAACAGTGGTGAATAACCACTCTCACCAAGGTAAAACTGAGTc                                                                                                                                                                                  >  1:468990/1‑251 (MQ=255)
                                                                 ggAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGAACTTTTTACCACCAGACGCAAGGTACTCTCGCTCATCCGGCGGCGAGGGGCGACAACTCCCAC‑‑CAGTGAGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGAGAACAGTGGGGAATAAACAATCTCAcccaagaaaaactgaggc                                                                                                                                                                                  >  1:506131/1‑235 (MQ=255)
                                                                  gAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGAAGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑CGCCGTTAACAGTGGGGAATAACCAATCTCACCCAGGGAAAGCTGAGTcg                                                                                                                                                                                 >  1:721723/1‑251 (MQ=255)
                                                                   aaaCCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTTCCACCAACCGCCAGTTACTCTCTATCATCCCGCGGCGGGAGGCGCCCACTCC‑CCA‑CACTGCGCG‑CCACGCAAAAAAAACTTCTTCTTT‑GCGCGTTAAACGTGGGGAATAACCAAtactcccccaggaaaacggtgtgcc                                                                                                                                                                                >  2:807675/1‑227 (MQ=255)
                                                                       cACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTTCCCTCTTCGATAAGGTCCAGCAACGCCAGTCCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTcc                                                                                                                                                                            <  2:515034/251‑1 (MQ=255)
                                                                         cGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCGATTTTTTCCCCCCGACACCAGTTACTTTCGGTCATCCGGCGGCGGGAGGCGGCATCTCC‑CCC‑CAGTGCTCG‑ACTCCAAAAATAAACTTCTTCTTC‑TGCCGTTAACACTGGTGAGTAACCAATCTCACCaacagaaaactgagtagtgaccc                                                                                                                                                                           >  2:402053/1‑229 (MQ=255)
                                                                             ccGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGGCCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAacaca                                                                                                                                                                      <  2:281113/251‑1 (MQ=255)
                                                                                   cAAACTACTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGAGTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCCATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCtgt                                                                                                                                                                <  2:1854/251‑1 (MQ=255)
                                                                                               ccGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCCATTTTTATCACCCAACGCAAGTTACACTCGATCATCCGGCGGCGAGAGGCGCCATCTCC‑ACA‑CCCGGCCCG‑CCCCCCCAAATAAACATCTTTTTC‑GGCCGTTTACAGTGGGGAATGAACAATCTCACCCAGGCTAAGCTGTGTTGCGTCCacaccaccccatgtggctcccctc                                                                                                                                                    >  2:668101/1‑228 (MQ=255)
                                                                                                    cgGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAATTTTAACACCCGAAACCAATTACTCTCGGACATCCGGCGGCGGGAGGGGACAGCTCC‑CCGGCCAT‑CGCG‑ACTCGCAAAAAAAACATCCTCTTC‑GGACGCTAAAAGGGGGGGAAAACAAATCTCAACAAGGGAAAGCTGAGGCGCGTCCAaaacacaccggggggccacccgcgata                                                                                                                                                >  2:715506/1‑229 (MQ=255)
                                                                                                        tCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTAACCCCCGACGCAAGATACTCTCGATCATCCGGCGGCGAGAGGCGAACTCTCC‑ACC‑CCATGCGCG‑ACTCGCAAAAAAAACTTCTTCTTCTGCCCGTTA‑CAGTGGGGGATAAACAATCTCACCAAGGTAAAGCTTAGTCGCGTCCCACCCCCCCTGTGTTGCTCCCTCCGATAAccgt                                                                                                                                           >  1:728041/1‑248 (MQ=255)
                                                                                                                       ccTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCaaa                                                                                                                            <  1:48213/251‑1 (MQ=255)
                                                                                                                          ctTCTATACGGTCCAGCCACTCCAGACCACGATTGCCACAACGGCCGTCAATTTTTACCAACAGACGCAAGTTACTCTCCATCATCCGGCGGCGAGAGGCGTCATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAATCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACTCGCTGTGTGGCACCCTGCGATAACAGTTCCTCTTCGTCCAAAt                                                                                                                           <  1:680874/249‑1 (MQ=255)
                                                                                                                          ctTCGATAAGGTCCAGCAACGCCAGGCCACGATTGCCATAACGGCTGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCACAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGATTAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAAt                                                                                                                           <  1:12284/249‑1 (MQ=255)
                                                                                                                            tCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGt                                                                                                                        <  1:703654/250‑1 (MQ=255)
                                                                                                                                        atcAACGCCAGACCACGATTGCCATAACGGCGGTCAATTTTTACCACCAGACGCAAGGTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑AAG‑CAGTGCGCG‑ACTCGCAAAATAACCTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTTTGGCGCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTc                                                                                                           <  2:88510/249‑1 (MQ=255)
                                                                                                                                             cGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCCCCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGtt                                                                                                      <  1:223841/251‑1 (MQ=255)
                                                                                                                                              gCCAGACCACGATGGCCACAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAGCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTc                                                                                                     <  2:721723/251‑1 (MQ=255)
                                                                                                                                                      ccGATTGCCATAACGGCGGGCAATTTTTACCCCCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACT‑CAGTGCGCG‑ACGCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAgg                                                                                             <  2:83862/250‑1 (MQ=255)
                                                                                                                                                           tGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACGCGCAAAATAAACTCCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCtt                                                                                         <  1:426667/250‑1 (MQ=255)
                                                                                                                                                            gcaagaaCGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTAC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTGATCACTGGGTTCCTGTTCTACTAAGGCCtttt                                                                                       <  2:778820/246‑1 (MQ=255)
                                                                                                                                                                         gCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTc                                                                          <  1:428022/251‑1 (MQ=255)
                                                                                                                                                                                ttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCg                                                                   >  1:317622/1‑251 (MQ=255)
                                                                                                                                                                                  accaccAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGtt                                                                 >  1:314060/1‑251 (MQ=255)
                                                                                                                                                                                  accaccAGACGCAAGTTACCCGCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTTCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCCCCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGtt                                                                 <  2:340150/251‑1 (MQ=255)
                                                                                                                                                                                          aCGCAAGTTACGCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCaa                                                         <  1:310066/251‑1 (MQ=255)
                                                                                                                                                                                            gCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAt                                                       >  1:290862/1‑251 (MQ=255)
                                                                                                                                                                                             caatttACTCTCGACCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAtt                                                      <  2:668840/247‑1 (MQ=255)
                                                                                                                                                                                               aGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTcc                                                    <  1:789590/251‑1 (MQ=255)
                                                                                                                                                                                                gtgaCTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCg                                                   <  1:777972/248‑1 (MQ=255)
                                                                                                                                                                                                   aCTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTCTCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCAt                                                >  2:715206/1‑251 (MQ=255)
                                                                                                                                                                                                   aCTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGGTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCAt                                                <  1:487086/251‑1 (MQ=255)
                                                                                                                                                                                                    cgtccGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATc                                               <  2:468990/247‑1 (MQ=255)
                                                                                                                                                                                                     tctcGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCt                                              >  1:462619/1‑251 (MQ=255)
                                                                                                                                                                                                        cGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGGCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCtt                                             <  1:402053/249‑1 (MQ=255)
                                                                                                                                                                                                         gATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGAt                                                                                                                                                 >  1:63354/1‑148 (MQ=255)
                                                                                                                                                                                                         gATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTCC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGAt                                                                                                                                                 <  2:63354/148‑1 (MQ=255)
                                                                                                                                                                                                                cggcggCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAACCCTCAACTCCGTTCTCATCAAATTCCGCCTCTTCCTTTAAATc                                   >  2:312252/1‑251 (MQ=255)
                                                                                                                                                                                                                   cggcgAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATg                                >  1:548595/1‑251 (MQ=255)
                                                                                                                                                                                                                          agGCGACATCTCC‑ACG‑CCGTGCGCG‑ACTTGCAAAATAAACTTCTGCTAC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAGACCGCATCTTCATTTAAATCAGGAACTTTc                         <  2:145040/251‑1 (MQ=255)
                                                                                                                                                                                                                             cGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGt                                                                                                             >  1:638502/1‑164 (MQ=255)
                                                                                                                                                                                                                             cGACATCTCC‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGt                                                                                                             <  2:638502/164‑1 (MQ=255)
                                                                                                                                                                                                                                     cc‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTAt                  <  1:174595/247‑1 (MQ=255)
                                                                                                                                                                                                                                     cc‑ACG‑CAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCAGCTTCATTTAAATCATGAACTTTCAGCGTAt                  <  1:631694/247‑1 (MQ=255)
                                                                                                                                                                                                                                            cAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTg                                                                                                                                                                                      >  2:545248/1‑78 (MQ=255)
                                                                                                                                                                                                                                            cAGTGCGCG‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTg                                                                                                                                                                                      <  1:545248/78‑1 (MQ=255)
                                                                                                                                                                                                                                                   cg‑ACTCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTg  >  1:789454/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       
GGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCC‑ACG‑CAGTGCGCG‑ACGCGCAAAATAAACTTCTTCTTC‑GGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTG  >  NC_000913/2867066‑2867557

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