Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,304,762 T→C 7.0% intergenic (+369/+346) eco → / ← mqo ecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,304,7620TC 7.0% 175.3 / 4.3 57intergenic (+369/+346)eco/mqoecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base T (29/24);  new base C (2/2);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.59e-02

TATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGG  >  NC_000913/2304495‑2305033
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                               
tatgcggCATACACCGCTGTGAAGGGCGGAACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATaa                                                                                                                                                                                                                                                                           <  2:280963/269‑1 (MQ=255)
                cTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCTGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGCTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCAGGAGGTTGGATCAGGCGGTTACGCCGCAt                                                                                                                                                                                                                                                           >  2:370295/1‑274 (MQ=255)
                cTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTCGGTTGGATAAGGCGGTTACGCCGCAt                                                                                                                                                                                                                                                           >  2:262957/1‑274 (MQ=38)
                cTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCAt                                                                                                                                                                                                                                                           >  2:345656/1‑274 (MQ=38)
                cTGTGAAGTGCGGCACCCCGTAGGACAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCCTGCGCGTCTTATCAGGCCTACCCCGCTGTTAAGTGCTCCCCCCCGTAGGTTGGATAAGGCGGTTACGCCGCaa                                                                                                                                                                                                                                                           >  2:445700/1‑273 (MQ=255)
                                                ggCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATg                                                                                                                                                                                                                          <  1:102602/275‑1 (MQ=35)
                                                                             cccGAGCGGGTTGCCTGATGCGACGCTGGCGCTTCTTATAAGGCCTACACCGCTTGGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                                                                                                                                                             <  2:350378/275‑1 (MQ=255)
                                                                                       ttGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGa                                                                                                                                                                                    >  2:490768/1‑274 (MQ=39)
                                                                                            cgATGCGCCGGTGGCCGGTCTGACCCGGCCTAACCCGCTGGGAAGTGCGCCACCCCGTAGGTCGGATAAGGCGGTTCGGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCt                                                                                                                                                                              <  4:345655/274‑1 (MQ=14)
                                                                                                                                 cTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGCTGAGGCGGTTACGCCGCAt                                                                                                                                          >  2:474609/1‑274 (MQ=255)
                                                                                                                                           agccACCCCGTGTAGCGGATAAGGCGGTTACGCTGCAGCCGAAATCAAACGCCCAAGGCGGTTGCCTGATGCGAGGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTACACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGGCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTa                                                                                                                                <  2:409822/272‑1 (MQ=14)
                                                                                                                                              cACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACAGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCTCATCCCACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGACGTGCTCCACCCCGTCGGTCGGCTAAGGCGGTTGCGCCGCATCCGACCTCTAAc                                                                                                                              >  1:108045/1‑273 (MQ=255)
                                                                                                                                                ccccGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTc                                                                                                                                                                                                           >  2:240520/1‑194 (MQ=34)
                                                                                                                                                                      ttACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCg                                                                                                                                                                                                                                          >  2:42621/1‑141 (MQ=37)
                                                                                                                                                                                    aCATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACg                                                                                                                                                                                                                     <  7:1048139/148‑1 (MQ=39)
                                                                                                                                                                                                 aaGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACa                                                                                                                                                                                                                                                     >  4:490342/1‑103 (MQ=32)
                                                                                                                                                                                                       gTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGgc                                                                                                                                                                                                                >  2:198778/1‑134 (MQ=37)
                                                                                                                                                                                                           ccTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                                                                                                                                            <  4:198777/134‑1 (MQ=37)
                                                                                                                                                                                                            ctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATa                                                                                                                                                                                                                                                                            <  6:322631‑M1/46‑1 (MQ=255)
                                                                                                                                                                                                            ctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATa                                                                                                                                                                                                                                                                            >  5:322631‑M1/24‑69 (MQ=255)
                                                                                                                                                                                                            ctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAgg                                                                                                                                                                                                                                                                         <  6:1111835‑M1/49‑1 (MQ=255)
                                                                                                                                                                                                            ctgataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAgg                                                                                                                                                                                                                                                                         >  5:1111835‑M1/24‑72 (MQ=255)
                                                                                                                                                                                                               ataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGt                                                                                                                                                                                                                                   <  7:359513‑M1/87‑1 (MQ=255)
                                                                                                                                                                                                               ataagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGt                                                                                                                                                                                                                                   >  8:359513‑M1/21‑107 (MQ=255)
                                                                                                                                                                                                                taagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCaa                                                                                                                                                                                                                                                >  8:587065‑M1/20‑93 (MQ=255)
                                                                                                                                                                                                                taagacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCaa                                                                                                                                                                                                                                                <  7:587065‑M1/74‑1 (MQ=255)
                                                                                                                                                                                                                   gacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTgg                                                                                                                                                                                                                                                                               >  7:804039‑M1/17‑59 (MQ=255)
                                                                                                                                                                                                                   gacgcgccagcgtcgcATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTgg                                                                                                                                                                                                                                                                               <  8:804039‑M1/43‑1 (MQ=255)
                                                                                                                                                                                                                             cgTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTg                                                                                                                                                                                                                             >  6:27506/1‑99 (MQ=255)
                                                                                                                                                                                                                             cgTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTg                                                                                                                                                                                                                             <  5:27506/99‑1 (MQ=255)
                                                                                                                                                                                                                             cgTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                                                                                                                                            >  8:144013/1‑116 (MQ=255)
                                                                                                                                                                                                                             cgTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                                                                                                                                            <  7:144013/116‑1 (MQ=39)
                                                                                                                                                                                                                               tCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                                                                                                                                            <  5:83122/114‑1 (MQ=40)
                                                                                                                                                                                                                               tCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                                                                                                                                            >  6:83122/1‑114 (MQ=40)
                                                                                                                                                                                                                                cTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGcc                                                                                                                                                                                                                                      <  7:1028536/87‑1 (MQ=255)
                                                                                                                                                                                                                                cTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGcc                                                                                                                                                                                                                                      >  8:1028536/1‑87 (MQ=255)
                                                                                                                                                                                                                                     CAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGgcgta                                                                                                                                                                                                                                                                     >  8:193607‑M1/1‑49 (MQ=255)
                                                                                                                                                                                                                                     CAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGgcgta                                                                                                                                                                                                                                                                     <  7:193607‑M1/51‑3 (MQ=255)
                                                                                                                                                                                                                                            acacCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCtgtg                                                                                                                                                                                     >  7:152788/1‑124 (MQ=40)
                                                                                                                                                                                                                                            acacCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGt                                                                                                                                                                                 >  7:88246/1‑128 (MQ=255)
                                                                                                                                                                                                                                            acacCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGt                                                                                                                                                                                 <  8:88246/128‑1 (MQ=255)
                                                                                                                                                                                                                                            acacCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCATACACCGCtgtg                                                                                                                                                                                     <  8:152788/124‑1 (MQ=255)
                                                                                                                                                                                                                                            acacCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCAGCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTGATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTAcg                               <  1:121371/274‑1 (MQ=38)
                                                                                                                                                                                                                                                cGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCa                                                                                                                                                                                                     <  7:37639/104‑1 (MQ=255)
                                                                                                                                                                                                                                                cGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCa                                                                                                                                                                                                     >  8:37639/1‑104 (MQ=255)
                                                                                                                                                                                                                                                        aGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCa                                                                                                                                                                                                                                                 <  7:1072914/52‑1 (MQ=255)
                                                                                                                                                                                                                                                        aGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCa                                                                                                                                                                                                                                                 >  8:1072914/1‑52 (MQ=255)
                                                                                                                                                                                                                                                          tGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCtg                                                                                                                                                                                       <  6:647219/108‑1 (MQ=255)
                                                                                                                                                                                                                                                          tGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCtg                                                                                                                                                                                       >  5:647219/1‑108 (MQ=255)
                                                                                                                                                                                                                                                               cACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTc                                                                                                                                                                                                           <  8:249077/83‑1 (MQ=255)
                                                                                                                                                                                                                                                               cACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTc                                                                                                                                                                                                           >  7:249077/1‑83 (MQ=255)
                                                                                                                                                                                                                                                               cACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACg            >  2:180241/1‑274 (MQ=38)
                                                                                                                                                                                                                                                                 ccccGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCtg                                                                                                                                                                                       >  3:123358/1‑101 (MQ=255)
                                                                                                                                                                                                                                                                   ccGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGt                                                                                                                                                                                                                                   >  6:116724/1‑55 (MQ=255)
                                                                                                                                                                                                                                                                   ccGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGt                                                                                                                                                                                                                                   <  5:116724/55‑1 (MQ=255)
                                                                                                                                                                                                                                                                     gTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGaa                                                                                                                                                                                   <  1:123358/101‑1 (MQ=255)
                                                                                                                                                                                                                                                                         gTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGTGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACCCCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCTGTTACGCCGCCTCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACCCCGCTGTTAAGTGCTCTCCCCCGACGGTTGGAAAAGGCGGTTCAGCCGGCCCCCACATCTAAccacgcagccgc  >  4:340861/1‑263 (MQ=11)
                                                                                                                                                                                                                                                                          ttGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCa                                                                                                                                                                                                     <  7:578746/78‑1 (MQ=255)
                                                                                                                                                                                                                                                                          ttGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCa                                                                                                                                                                                                     >  8:578746/1‑78 (MQ=255)
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                               
TATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGG  >  NC_000913/2304495‑2305033

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: