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,5960GA97.6% 131.2 / NA 41G32G (GGC→GGTgfcCcapsule biosynthesis GfcC family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (17/23);  minor base C (1/0);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTG  >  NC_000913/1048353‑1048837
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 
aGTGATATTGAGAATGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCAAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCCACGTAAAGTGTCTGTTGTTCACCAGGc                                                                                                                                                                                                                                             >  2:107865/1‑250 (MQ=255)
          agattgcgTAGCTGCTGGCGCTCGGATGTACTCGTTGCCGCCACATCGTCGTCAGCTTCTGCGTCCCATGAAGCAAGCTGTGCCATGACATGTTGGTAGGCTTTAAGCGCTGTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCATAGGCGATCGCGCAGGTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTTTCAGTTGTTCACCATGCAGATAAATAGt                                                                                                                                                                                                                                  <  2:272853/251‑1 (MQ=255)
              tgagTAACTGCAGGCGCACGGCTTTAATCGTTGCCGCCACATCGTCGTAAGCTTCTGATTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACAAGGCCACCAGAGGCGATCGCGCAGTTGAGGTAGTGTCACCAACTGGACACCATTTTCCACGAGGCCACCGGAAAGTGTCTGTTGTTCACAAAGCAGATAAAAAGTCAcc                                                                                                                                                                                                                              >  2:279225/1‑251 (MQ=255)
               gagTAACTGCTGTCGCACGGATTTAATCATTGCCGCCACTTCGTCGTCAGCTTCTGCTCCCCAGGAAGCAAGCTGTGCCATGGCATGTTGTTAGGCATTAAGCGCTGTAGCGTTCGCCGCAGAGTCGGGCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTGGTGTCTCCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGGTGTTCACCAGGCAGATAAATAGTCAcc                                                                                                                                                                                                                              <  2:316894/250‑1 (MQ=255)
                agTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGGAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGt                                                                                                                                                                                                                            >  1:298259/1‑251 (MQ=255)
                  tAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTc                                                                                                                                                                                                                          >  1:55280/1‑251 (MQ=255)
                       tctgGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTg                                                                                                                                                                                                                      <  1:338237/249‑1 (MQ=255)
                       gctgGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCtt                                                                                                                                                                                                                       <  1:8957/249‑1 (MQ=255)
                               aCGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCTGGAAGACAGCTGTGCAATGACATGTTGGTAGTCTTTAAGCGCTTTCGATTTCGACGCAGAGTCGGTCAGCAAGGCACCAGGCCACCCGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCCCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGc                                                                                                                                                                                                             >  1:433579/1‑251 (MQ=255)
                                   aTTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATgg                                                                                                                                                                                                         <  1:264338/251‑1 (MQ=255)
                                       aaTCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGgtg                                                                                                                                                                                                     >  1:432900/1‑251 (MQ=255)
                                        aTCGTTGCCGCCTCATCGTCGGCAGCTTCTGCTTCCCAGGAAGTGAGCTGTGCCATGACAGGTTGGTAGTCTGTAAGCGCTTCAGCTGTCGCCGCAGAGGCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGTCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATATTCACCGTTCCTTGCGCAAATTCATGGGgtgt                                                                                                                                                                                                    <  2:298259/251‑1 (MQ=255)
                                        aTCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGTAAGATGTGCCATGACACGTTGGTCGTCTCTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCACGGCCCCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGgtgt                                                                                                                                                                                                    >  1:469088/1‑251 (MQ=255)
                                              gccgccACATTGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTGCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTAc                                                                                                                                                                                              <  2:349113/251‑1 (MQ=255)
                                              gccgccACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGCGCCATGACATGTTGGCAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTAc                                                                                                                                                                                              >  1:233744/1‑251 (MQ=255)
                                               ccgccACATCGTCGTCAGCTTCTGCGTCCCAGGAAGCAAGCTGGGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACg                                                                                                                                                                                             <  2:90754/251‑1 (MQ=255)
                                                cgccACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTGACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGt                                                                                                                                                                                            <  1:279225/251‑1 (MQ=255)
                                                  ccACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTaa                                                                                                                                                                                          <  1:174256/251‑1 (MQ=255)
                                                    acaTCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGCAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTaa                                                                                                                                                                                          >  1:37406/1‑249 (MQ=255)
                                                                gCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAAt                                                                                                                                                                            <  1:252158/251‑1 (MQ=255)
                                                                                 agcaagcTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTGCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTGCGTAAAGTACGCTGGCAATGAAATACGACTGTAAtt                                                                                                                                                           <  2:352239/251‑1 (MQ=255)
                                                                                        tgtgCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCAt                                                                                                                                                    >  1:409480/1‑251 (MQ=255)
                                                                                                gACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGGCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTGGCGGAAAGTCATGGGTTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCgg                                                                                                                                             <  2:432900/250‑1 (MQ=255)
                                                                                                           aGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGAACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGAGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGCAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGGTTAAGAAACg                                                                                                                                 >  1:203076/1‑251 (MQ=255)
                                                                                                                                  gccgcAGAGTCGGTCAGCAAGGCACCAGGCCACCAGCGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCt                                                                                                          >  2:456893/1‑251 (MQ=255)
                                                                                                                                  gccgcAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAAGAGTCACCGTTCCTTGCGCAAAGGCATGGTGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCAGGGTTCGGGTTGAAGAAACGTCATTTCTACGGGAATAACCCCt                                                                                                          <  2:176397/251‑1 (MQ=255)
                                                                                                                                             ggTCAGCAAGGAACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGCCAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCCTGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGCCTGTAATTTATTCATGGTGCGGGTTTAAGCAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACAt                                                                                               >  2:317650/1‑251 (MQ=255)
                                                                                                                                                cAGCAAGGCACCAGGCCACCAGGGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACCGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCCAAGGCATGGGGAGTCATTACGTAAAGTACGCTGGCAATGAAATACGAATGTAATTTATTCATGGTGCGGGTTTAAGAACCGTCATTTCTCCGGGAATAACCCCTGCGCCTAACATTTg                                                                                            >  2:292546/1‑251 (MQ=255)
                                                                                                                                                        cACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACttt                                                                                    <  2:409480/251‑1 (MQ=255)
                                                                                                                                                          ccAGGCCTCCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTACTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGa                                                                                  <  2:42327/251‑1 (MQ=255)
                                                                                                                                                                caccaGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAAcc                                                                            >  1:1793/1‑251 (MQ=255)
                                                                                                                                                                      cagCGATCGCGCAGTTGCTGTTGTGTCACCTACTGGACAACATTTTCCCCGGGGCCAACGGAAAGTGTCTGTTGGTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCAGACTTTGACGAACCTGGCCa                                                                      <  2:9725/249‑1 (MQ=255)
                                                                                                                                                                                 cAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCAGTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATc                                                           <  1:423377/251‑1 (MQ=255)
                                                                                                                                                                                           ggtgtCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCATGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGAACCAGAAt                                                 <  2:365110/250‑1 (MQ=255)
                                                                                                                                                                                              gtCACCCACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGGAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCAGGGTGCGGGTTTAAGAAACGTCATTGCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTGTGAGTAACGTGGCCAGTAGGCGAATCGACCCAGAAtgt                                               <  2:441270/250‑1 (MQ=255)
                                                                                                                                                                                                                    ccACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgact                         <  1:292546‑M1/250‑11 (MQ=255)
                                                                                                                                                                                                                           gCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGagtgtggactgagtaacc                 <  2:55280‑M1/251‑19 (MQ=255)
                                                                                                                                                                                                                               aCGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagac             <  1:317650‑M1/251‑23 (MQ=255)
                                                                                                                                                                                                                                      gtgtCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCCAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCCCTGattgtggactgagtaaccagacaggtgat      >  2:664‑M1/1‑222 (MQ=255)
                                                                                                                                                                                                                                          cTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCCAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgc  >  1:308644‑M1/1‑218 (MQ=255)
                                                                                                                                                                                                                                          cTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGACGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgc  >  2:115422‑M1/1‑218 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 
AGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTG  >  NC_000913/1048353‑1048837

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