Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,115,847 G→T 7.6% A17E (GCG→GAG)  wcaM ← colanic acid biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,115,8470GT 7.6% 265.7 / 10.3 79A17E (GCG→GAG) wcaMcolanic acid biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base G (36/37);  new base T (2/4);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.47e-01

ACCTGGCTGCCGTCCTGTAAGACAAATCGCCCTCTGCCGTTGCCGCGCAAGCTTCCAAGAATGTGTAACGTTTTACCGGGAGGGATAAAGATGCCGGTGTTGATATTGTCACAAACCAATCCGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCGTTGCGCCAGTGCGCGAGCATCGTTCTCAGGCACCAGCCAGCCGGATTTGTCAGCCTCCACCAGTTCCGGTATTCCAC  >  NC_000913/2115586‑2116118
                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                               
aCCTGGCTGCCGTCCTGTAAGACAAATCGCCCTCTGCCGTTGCCGCGCAAGCTTCCAAGAATGTGTAACGTTTTACCGGGAGGGATAAAGATGCCGGTGTTGATATTGTCACAAACCAATCCGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGt                                                                                                                                                                                                                                                                    >  2:420286/1‑275 (MQ=255)
                         aTCGCCCTCTGCCGTTGCCGCGCAAGCTTCCAAGAATGTGTAACGTTTTACCGGGAGGGATAAAGATGCCGGTGTTGATATTGTCACAAACCAATCCGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGGTGGGAGGGtttt                                                                                                                                                                                                                                           >  2:107215/1‑275 (MQ=255)
                                  ttccGTTGCCGCGCAAGCTTCCAAGAATGTGTTACGTTTTACCGGGAGGGATAAAGATTCCGGTGTTGATATTGTCACAAACCAATCCGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGc                                                                                                                                                                                                                                  <  2:142562/273‑1 (MQ=255)
                                   gCCGTGGCCGCTGAAGCGTTCAAGAGTGTGTAACGTAGTACCGGGAGGGATAAAGATGCCGGGGTCGAGATATGCACGAACCAATCCGTCAGGGACGTGGACTGGTTGACCTTAGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTGGTGAATGTTAACGCCGTGTCGGGCGGGAAGCGGGCGGGCGAAAGGGGTATGGGGGAAGGCTAGCGCCGGGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGc                                                                                                                                                                                                                                  <  4:339209/274‑1 (MQ=255)
                                    ccGTTGCCGCGCAAGCTTCCAAGAATGTGTAACGTTTTACCGGGAGGGATAAAGATGCCGGTGTTGATATTGTCACAAACCAATCCGGCAGGCACGACGACTTTTTTCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGGGGGTTGTAGTCGTTACTGTTAACGATTGGTCGGGCGGGAAGCCCGGGGGCGAAAGGGGAATGGAGGAAGGAAAGCGCCGAACTTGCCGTCAGGAACGTGCGACGGGCAGGTTTTTTAAATGGCAt                                                                                                                                                                                                                                >  4:235470/1‑275 (MQ=255)
                                                                   aCGTTTTACCGGGAGGGATAAAGATGCCGGTGTTGATATTGTCACAAACCAATCCGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCtat                                                                                                                                                                                                                  <  2:16267/258‑1 (MQ=255)
                                                                              ggagggaTAAAGATGCCGGTGTTGATATTGTCACAAACCAATCCGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTtc                                                                                                                                                                                        >  2:121076/1‑273 (MQ=255)
                                                                               gagtgaTAAAGAGGCCTGTGTTGAGATTGTCACAAACCAATCCGGCAGGCAGGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTtctc                                                                                                                                                                                      <  2:229693/274‑1 (MQ=255)
                                                                                gtggaTAAATGTGCCGGTGTGGATATTGTCAAAAACCAATCCGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGa                                                                                                                                                                                    <  1:69259/273‑1 (MQ=255)
                                                                                agggaTAAAGATGCCGGTGTTGATATTGTCACAAACCAATCCGGCAGGCACGACGACTGTTTGCCCTTCACTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGa                                                                                                                                                                                    >  2:98057/1‑275 (MQ=255)
                                                                                  ggaAAAAGATTCCGGTGTTGATATTGTCACAAACCAATCCGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGAtt                                                                                                                                                                                  <  2:403816/275‑1 (MQ=255)
                                                                                            gCCGGTGGGGAGATGGTCACAAACCAATCCGGGAGGAAAGACGACTGTTTGCCCTTCGCTGAAGGATTTGTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAAGGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTTCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctg                                                                                                                                                                         <  4:420281/274‑1 (MQ=255)
                                                                                                                     aaTCCGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTCGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGa                                                                                                                                                                                                                                                                             <  8:389622/149‑1 (MQ=255)
                                                                                                                               gcaCGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCTCGGCGCTTGCCTTc                                                                                                                                                                                                                                                                   <  6:303051/149‑1 (MQ=255)
                                                                                                                                caCGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGACGGCAAGCTCGGCGCTTGCCTTcc                                                                                                                                                                                                                                                                  <  8:288212/149‑2 (MQ=255)
                                                                                                                                    acgacTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTga                                                                                                                                <  2:185540/275‑1 (MQ=255)
                                                                                                                                         tGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGc                                                                                                                                                                                                                                                         >  6:225584/1‑149 (MQ=255)
                                                                                                                                           tttGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAGAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCANGAACGGGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAACCTGCAGCAAGCTGTCGAGTTCTCGATTTNAACCCTGCTGGTTAAAATCGCGTTCNACTTTTTCGTGCGCTCCttttcagcccgg                                                                                                                         >  3:93066/1‑266 (MQ=255)
                                                                                                                                              gCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCgg                                                                                                                                                                                                                                                    <  7:150987/149‑1 (MQ=255)
                                                                                                                                               cccTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCggg                                                                                                                                                                                                                                                   >  7:267550/1‑149 (MQ=255)
                                                                                                                                               cccTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCggg                                                                                                                                                                                                                                                   <  6:455648/149‑1 (MQ=255)
                                                                                                                                                ccTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGga                                                                                                                                                                                                                                                  <  8:242456/149‑1 (MQ=255)
                                                                                                                                                ccTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGga                                                                                                                                                                                                                                                  <  7:27532/149‑1 (MQ=255)
                                                                                                                                                   tCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTTCAGCAAGCTGGCGAGTtct                                                                                                                                                                                       >  2:81484/1‑205 (MQ=255)
                                                                                                                                                      cTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGttt                                                                                                                                                                                                                                            <  6:257080/149‑1 (MQ=255)
                                                                                                                                                       tGAAGGCTTGTTTAAATGAGGCGATCCAGGCGTGTGGGTTGTAGTCGTTAATGGTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATtcgt                                                                                                              <  1:48952/274‑1 (MQ=255)
                                                                                                                                                        gAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGttttt                                                                                                                                                                                                                                          <  8:420934/149‑1 (MQ=255)
                                                                                                                                                         aaGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGtttttt                                                                                                                                                                                                                                         >  5:130733/1‑149 (MQ=255)
                                                                                                                                                            gCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTaaa                                                                                                                                                                                                                                      >  6:462239/1‑149 (MQ=255)
                                                                                                                                                             cTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAAt                                                                                                                                                                                                                                     <  5:462239/149‑1 (MQ=255)
                                                                                                                                                                          ggCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGCAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTtctc                                                                                                                                                                                                                        <  7:50040/149‑1 (MQ=255)
                                                                                                                                                                           gCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTtc                                                                                                                                                                                        >  1:24840/1‑180 (MQ=255)
                                                                                                                                                                           gCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAAGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAAcg                                                                                         >  2:86145/1‑275 (MQ=255)
                                                                                                                                                                               tCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTCTAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAAcgccgc                                                                                     >  2:29681/1‑275 (MQ=255)
                                                                                                                                                                               tCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCCTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTAGGGAGAGTTTTTTAAATTGCATGCGTTCTCCTCta                                                                                                                                                                                                                   >  8:288751/1‑149 (MQ=255)
                                                                                                                                                                               tCCAGGTGTGTGGGTTGTAGTCGTTAATGTTAACGCTGTGTNGGGCGGNAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGa                                                                                                                                                                                    <  3:24840/180‑1 (MQ=255)
                                                                                                                                                                                       tgtgGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCt                                                                                                                                                                                                           <  8:267550/149‑1 (MQ=255)
                                                                                                                                                                                            gTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTctcctc                                                                                                                                                                                                                     >  3:3004/1‑134 (MQ=255)
                                                                                                                                                                                             ttGTGGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCtata                                                                                                                                                                                                                 <  7:207816/137‑1 (MQ=255)
                                                                                                                                                                                             ttGTGGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCtata                                                                                                                                                                                                                 >  8:207816/1‑137 (MQ=255)
                                                                                                                                                                                              tGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCaa                                                                                                                                                                                                    <  5:271628/149‑1 (MQ=255)
                                                                                                                                                                                                tAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCtata                                                                                                                                                                                                                 <  1:3004/134‑1 (MQ=255)
                                                                                                                                                                                                  gTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTg                                                                                                                                                                                                >  6:162461/1‑149 (MQ=255)
                                                                                                                                                                                                  gTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTg                                                                                                                                                                                                >  6:61880/1‑149 (MQ=255)
                                                                                                                                                                                                     gTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCa                                                                                                                                                                             >  2:328926/1‑165 (MQ=255)
                                                                                                                                                                                                      ttAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGa                                                                                                                                                                                            >  7:405831/1‑149 (MQ=255)
                                                                                                                                                                                                      ttAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGa                                                                                                                                                                                            >  5:45219/1‑149 (MQ=255)
                                                                                                                                                                                                         aTGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctg                                                                                                                                                                         <  4:328922/165‑1 (MQ=255)
                                                                                                                                                                                                                cgcttgggggggggggAAGCGCGGGGTGGGATGGGGGATGGAGGAAGGCAAGCGCGTGGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGAATGCGTTCTGCGCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATTCCCTGTTGGTTAAACTCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCGTTGCGCCAGTGCGCGAGCATCGTtct                                                     <  4:272989/265‑1 (MQ=255)
                                                                                                                                                                                                                cGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGAtt                                                                                                                                                                                  >  5:95078/1‑149 (MQ=255)
                                                                                                                                                                                                                   tttGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGAt                                                                                                                                                                               <  5:61880/149‑1 (MQ=255)
                                                                                                                                                                                                                   tttGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGAt                                                                                                                                                                               <  8:189754/149‑1 (MQ=255)
                                                                                                                                                                                                                     tGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAACCGCCGCCAGTCGTTGCGCCAGTTCGCGAGCATCGTTCTCAGGCa                                               >  2:297862/1‑275 (MQ=255)
                                                                                                                                                                                                                     tGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCGTTGCGCCAGTGCGCGAGCATCGTTCTCAGGCa                                               >  2:319321/1‑275 (MQ=255)
                                                                                                                                                                                                                            cgggAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctgctg                                                                                                                                                                      <  8:405831/149‑1 (MQ=255)
                                                                                                                                                                                                                                 aGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCGTTGCGCCAGTGCGCGAGCATCGTTCTCAGGCACCAGCCAGCCgg                                   >  2:151124/1‑275 (MQ=255)
                                                                                                                                                                                                                                 aGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCGTTGCGCCAGTGCGCGAGCATCGTTCTCAGGCACCAGCAAGCCg                                    >  2:479263/1‑274 (MQ=255)
                                                                                                                                                                                                                                        gggCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCgtg                                                                                                                                                          <  5:225584/149‑1 (MQ=255)
                                                                                                                                                                                                                                         ggCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCgtgt                                                                                                                                                         <  6:95078/149‑1 (MQ=255)
                                                                                                                                                                                                                                         ggCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCGTTCCGCCAGTGCGCGAGCATCGCTTTCAGGCACCAGCCAGCCGGGTTTGtg                            >  2:436794/1‑273 (MQ=255)
                                                                                                                                                                                                                                            gAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTCAAAATCGTGTTCa                                                                                                                                                      >  8:301687/1‑149 (MQ=255)
                                                                                                                                                                                                                                                  ggTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCGTTGCGCCAGTGCGCGAGCATCGTTCTCAGGCACCAGCCAGCCGGATTTGTCAGCCTccac                   <  1:93066/274‑1 (MQ=255)
                                                                                                                                                                                                                                                    tATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTc                                                                                                                                              <  8:4903/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    tATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTc                                                                                                                                              <  5:162461/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    tATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTc                                                                                                                                              >  7:222668/1‑149 (MQ=255)
                                                                                                                                                                                                                                                       ggaggaAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTcgcg                                                                                                                                           >  5:221454/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        gaggaAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGGGTCGGGAGAGTTTTTTAATTGGCATGCATTCTCCCCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCGTTGCGCCAGTGCGCGAGCATCGTTCTCAGGCACCAGCCAGCCGGATTTGTCAGCCTCCACCAGTTc             <  2:335947/274‑1 (MQ=255)
                                                                                                                                                                                                                                                        gaggaAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGGTCACCTGCTGGTTAAAATCGTGTTCAACTTTTTcgcgc                                                                                                                                          >  8:36864/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        gagcaAGGAAAGCGCGGCGCGGGCCGTCAGGAGCGTGCGTTGGTAGAGTTTTTTAGATGGCGTCCGTTCTCCTATATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTGTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCGTTGCGCCAGTGCGCGAGCATCGTTCTCAGGCACCAGCCAGCCGGATTTGTCAGCCTCCACCAGTTc             <  4:189341/274‑1 (MQ=255)
                                                                                                                                                                                                                                                          tgacggCAAGCTCGGCGCTTGCCTTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACtttt                                                                                                                                                >  7:317627/5‑141 (MQ=255)
                                                                                                                                                                                                                                                          tgacggCAAGCTCGGCGCTTGCCTTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACtttt                                                                                                                                                <  8:317627/137‑1 (MQ=255)
                                                                                                                                                                                                                                                          ggagtgCAAGGGCCGTGGTTGGCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTGTAAAGCCGGCAGGAAGCTGGCGAGTTCTCGATTGATCACCGGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCGTTGCGCCAGTGCGCGAGCATCGTTCTCAGGCACCAGCCAGCCGGATTTGTCAGCCTCCACCAGTTCCg           <  4:319317/269‑1 (MQ=255)
                                                                                                                                                                                                                                                            aaGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCGTTGCGCCAGTGCGCGAGCATCGTTCTCAGGCACCAGCCAGCCGGATTTGTCAGCCTccacca                 >  2:189343/1‑266 (MQ=255)
                                                                                                                                                                                                                                                             aGGCGAGCGCCGACCTTGCCGTGCGGAGAGTGGATAGGGAGAGATGGTTAAATGGCTGTCTCTGTCCTCTACAAAACGTTCGGTAAGTCTGCAACGTTTCGTGTGATCACCTTCTGCTTAAAGTTGTGTTTGACTTTTACGCGCTAGAGGTTGACAACAGGTGCCGAGTCGTCGGTGTCCAGTTGTGTAAAGGCCGCCAGTCGTTTAGACAGTGCGCGGGCATCGTTCTCAGGCACCAGCCAGCCGGATTGGTCAGCCTCCACCAGTTCCGGTa        <  4:436789/274‑1 (MQ=255)
                                                                                                                                                                                                                                                              ggCAAGCTCGGCGCTTGCCGTCAGGAACGTGCGTCGGTAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGt                                                                                                                                    >  5:389623/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              ggCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCACTTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCGTTGCGCCAGTGCGCGAGCATCGTTCTCAGGCACCAGCCAGCCGGATTTGTCAGCCTCCACCAGTTCCGGTAtt      >  2:456296/1‑275 (MQ=255)
                                                                                                                                                                                                                                                              ggCAAGCGCCGAGCTTGCAGCCAGGAACGTGCGTCGGGAGAGGTTTTTTAATGGCATGCGTTCTTCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGCTCACCTGCTGGTTAAACCCGGGTTCAACTCTTTCGCGCGCGCGTTTGACCACCGGAGCCACTTCGTCGGTTTCCAGTTCGCTAAGCCCCCCCAGCCCTTGCGACCGTGCGCGAGCATCGGTCTCAGGCCCCATCCCACCGGATTTTTCCCCCCCCACCCGCTCTGGCAtt      >  4:248320/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                cAAGCTCGGCGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGttt                                                                                                                                  <  6:389623/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAACTCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCGTTGCGCCAGTGCGCGAGCATCGGTCTCAGGCACCAGCCAGCCGGATTTGTCAGCCTCCACCAGTTCCGGTATTCCAc  >  2:160028/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                     gcCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTgacga                                                                                                                             >  6:83775/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                               
ACCTGGCTGCCGTCCTGTAAGACAAATCGCCCTCTGCCGTTGCCGCGCAAGCTTCCAAGAATGTGTAACGTTTTACCGGGAGGGATAAAGATGCCGGTGTTGATATTGTCACAAACCAATCCGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCGTTGCGCCAGTGCGCGAGCATCGTTCTCAGGCACCAGCCAGCCGGATTTGTCAGCCTCCACCAGTTCCGGTATTCCAC  >  NC_000913/2115586‑2116118

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑