Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,048,596 G→A G32G (GGC→GGT gfcC ← capsule biosynthesis GfcC family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,048,5960GA98.0% 170.1 / NA 50G32G (GGC→GGTgfcCcapsule biosynthesis GfcC family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (22/27);  minor base C (1/0);  total (23/27)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGA  >  NC_000913/1048357‑1048842
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      
atatTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATAGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGAt                                                                                                                                                                                                                                              >  1:25397/1‑250 (MQ=255)
              tAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTc                                                                                                                                                                                                                               >  2:388200/1‑251 (MQ=255)
                aCTGCTGGCGCCAGGATTTAATCGTTGCCGCCACCTCGTCGTCAGATTATGCTTCCCAGGAAGCAAGTTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCACGGCACCAGGCCACACGAGGCGATCGCGCAGTTGCGGTTGTGTCGCCAACTGGACACCACTTTCCACGGGGGCACCGGCAAGTGTCTTTTGTGCCGCAGGCAGGTAAATACTCCCCGTTCCt                                                                                                                                                                                                                             >  2:234886/1‑251 (MQ=255)
                   gctgGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGCGTCTTTTGTTCACCAGGCAGATAAATAGTCAGCGTTCCTTg                                                                                                                                                                                                                           <  1:162110/250‑1 (MQ=255)
                      ggCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTGTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCGGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCa                                                                                                                                                                                                                       <  1:441218/251‑1 (MQ=255)
                           aCGGATTTAATCGTTGACGCCACATCGTCGTCAGCTTCTGCTTCCCAGGCAGCAAGATGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCCACTGGACAACATTTTCCACGGGGCCAACGGACAGGGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGc                                                                                                                                                                                                                  >  2:284279/1‑251 (MQ=255)
                                 ttAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATggg                                                                                                                                                                                                             <  2:262725/250‑1 (MQ=255)
                                    aTCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAAGAGTCACCGTTCCTTGCGCAAAGGCATGGGgtgt                                                                                                                                                                                                         <  2:25397/251‑1 (MQ=255)
                                     gcgtggccgccTCAGCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACAGGTTGGTAGTCTGTAAGCGCTTGAGCTTTCGCCGCAGAGACGGTCAGCAAAGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTGCCACGGGGCCTACGGAAAGGGTCTGTGGTTCACCAGGCAGATAAATAGTCACCGTTCCGCACGCAAAGGCGTGGGGTGTc                                                                                                                                                                                                        <  2:282465/246‑1 (MQ=255)
                                      cGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTc                                                                                                                                                                                                        <  2:350406/250‑1 (MQ=255)
                                          gccgccACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTAc                                                                                                                                                                                                   <  2:144173/251‑1 (MQ=255)
                                                          cAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGc                                                                                                                                                                                    >  1:262536/1‑250 (MQ=255)
                                                               ccTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGGTCACCAGTCAGATAAATAGTCACCGTACCTTGCGCAAAGGCATGGGGTGTCAGTACGTAAAGTACGCTGGCAATGaa                                                                                                                                                                              <  2:447798/250‑1 (MQ=255)
                                                                cTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGAGCTTTAGCTTTAGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCAACCACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCCATGaa                                                                                                                                                                              >  1:447798/1‑250 (MQ=255)
                                                                    ttCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACg                                                                                                                                                                         >  1:290808/1‑251 (MQ=255)
                                                                     tCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTGAAGCGCTTTAGCTTGCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGGTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGGCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTGCGCTGGCAATGAAATAc                                                                                                                                                                          <  1:473331/249‑1 (MQ=255)
                                                                             agcaagcTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTTTAAtt                                                                                                                                                                <  2:150849/251‑1 (MQ=255)
                                                                                aagatgtgCCATGACCTGTTGGTAGTCTTTACGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCCAGGCACCAGGCCACCCGAGGCGATCGCGCAGTTGCGTTTGTGTCACCAACTGCACAACCTTTTCCACGGGGCCAACGGCAAGTGTCTGTTGTTCACCAGGCAGATAAATCGTCAACGTTCCCTCCGCAACGGCATGGGGTGTCATTACGTAAAGGACGCTGGCAATGAAATACgtaggacatttat                                                                                                                                                             >  2:408496/5‑239 (MQ=255)
                                                                                 agcTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGCAATTTAtt                                                                                                                                                            >  2:479358/1‑251 (MQ=255)
                                                                                 agcTGTGCCATGACATGTTGGTAGTATTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTAtt                                                                                                                                                            >  1:143839/1‑251 (MQ=255)
                                                                                  gcTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAACGTACGCTGGCAATGAAATACGACTGTAATTTATTc                                                                                                                                                           >  1:485202/1‑251 (MQ=255)
                                                                                  gcTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACAAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACCACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTAACCGTTACTTGCGCAAAGGCATGGGGTGTCATTCCGTAAAGTACGCTGGCAATGAAAt                                                                                                                                                                            >  2:429601/1‑234 (MQ=255)
                                                                                  gcTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTGTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGGCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGGAAAGTACGCTGGCAATGAAAt                                                                                                                                                                            <  1:429601/234‑1 (MQ=255)
                                                                                          aTGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCg                                                                                                                                                   <  2:208888/251‑1 (MQ=255)
                                                                                                    ggTAGTCTTTAAGAGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGGTCACCAGGCAGATAAATAGTCAACGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGGACGCTGGCAATGAAAGACGACTGTAATTTTTTCATGGTGCGGGTTTAAGa                                                                                                                                          >  2:190915/1‑250 (MQ=255)
                                                                                                        gTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGt                                                                                                                                     >  1:401975/1‑251 (MQ=255)
                                                                                                          cTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACCACATTTTCCCCGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTCACCAAAGGCGTGGGGGGTCATTACGTAACGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTc                                                                                                                                    >  2:87615/1‑250 (MQ=255)
                                                                                                            ttAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTGGCTGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGGCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAATCGTCa                                                                                                                                   <  2:404905/249‑1 (MQ=255)
                                                                                                                     ttAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCTATCGCGCAGTTGCGGGTGTGTCTCCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGGAAAGTACTCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTGTAAGAAACTTCATTTCTACggg                                                                                                                         <  2:215378/250‑1 (MQ=255)
                                                                                                                              gccgcAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCt                                                                                                               <  1:400923/251‑1 (MQ=255)
                                                                                                                              gccgcAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCt                                                                                                               >  2:29240/1‑251 (MQ=255)
                                                                                                                              gccgcAGAGTCGGGCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAAGAGGCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGGACGCTGGCAATGAAATACGACGGTAATTGATTCATGGGGCGGGTTGAAGAAGCTTCATTTCTACGGGAATAACCCCt                                                                                                               <  2:238812/251‑1 (MQ=255)
                                                                                                                                 gcAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGGCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCTGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTgcg                                                                                                            <  1:55759/251‑1 (MQ=255)
                                                                                                                                         ggTCAGCAAGGCACCAGGCCACCAGAGGCGCTCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCCGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACAt                                                                                                    >  1:75859/1‑251 (MQ=255)
                                                                                                                                            cAGCAAGGCACCATGCCCCCAGAGGCGATCGCTCAGTTGCGTTTGTGTCCACAACTGTACAACATTTTCCACGGGTCCAACGGAAATTGTCTGTTGTTCACCAGTCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACTACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTg                                                                                                 <  1:234886/251‑1 (MQ=255)
                                                                                                                                                       cAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAATAGTCAACGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTAAGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACg                                                                                     >  1:406492/1‑251 (MQ=255)
                                                                                                                                                         ggCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTACACGGGGCCAAAGGACAGTGTATGTTGTTCACCAGGCAGATAAATAGACACAGTTAATTGCGCAAAGGCATGGGGCGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACg                                                                                     >  2:484490/1‑250 (MQ=255)
                                                                                                                                                             cccaGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCTATTCTACGGGAATAACCCCTGCTCTTAACATTTGCCGACTTTGACGAACCt                                                                                <  2:313707/250‑1 (MQ=255)
                                                                                                                                                                             cAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATc                                                                <  2:392778/251‑1 (MQ=255)
                                                                                                                                                                                 tGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTGGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACAGTGGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGAc                                                             <  1:330567/250‑1 (MQ=255)
                                                                                                                                                                                 tGCGGTTGTGCCACCACCTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGAAGATAAATAGTCAACGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTACAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGCGCGGGTTTAAGAAAc                                                                                                                                       >  2:47455/1‑176 (MQ=255)
                                                                                                                                                                                 tGCGGTTGTGCCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGTCAGATAAATAGGCACCGGTCCTTGCGCAAAGGCATGGGGTGGCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAAc                                                                                                                                       <  1:47455/176‑1 (MQ=255)
                                                                                                                                                                                            cACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGtt                                                 <  2:290808/251‑1 (MQ=255)
                                                                                                                                                                                              ccAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGATCCTGGCCAGTAGCCGAATCGACCCAGAATGTGtttt                                               <  1:281127/251‑1 (MQ=255)
                                                                                                                                                                                                   gggACAACATTTTCCACGGGGCCGACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAAGAGTCACCGTTCCTTGCGCAAAGGCATGGGGGGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAGTAACCCCTGCGCCTAACATGTGCCGACTGTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCAc                                          <  2:407303/250‑1 (MQ=255)
                                                                                                                                                                                                                 cACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactg                             >  2:398352‑M1/1‑239 (MQ=255)
                                                                                                                                                                                                                    ggggCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagta                         <  2:262536‑M1/251‑16 (MQ=255)
                                                                                                                                                                                                                             ggAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACTTTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagac                  >  1:259937‑M1/1‑227 (MQ=255)
                                                                                                                                                                                                                              gAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagt               >  1:187912‑M1/1‑226 (MQ=255)
                                                                                                                                                                                                                                           gttCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgat  <  2:360995‑M1/251‑39 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      
ATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGA  >  NC_000913/1048357‑1048842

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