Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,718,081 A→G 6.7% L171P (CTT→CCT)  paaG ← 2‑(1,2‑epoxy‑1,2‑dihydrophenyl)acetyl‑CoA isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,718,0810AG 6.7% 171.3 / 2.9 60L171P (CTT→CCT) paaG2‑(1,2‑epoxy‑1,2‑dihydrophenyl)acetyl‑CoA isomerase
Reads supporting (aligned to +/- strand):  ref base A (22/34);  new base G (2/2);  total (24/36)
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.03e-01

GCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCA  >  NC_002947/3717937‑3718229
                                                                                                                                                 |                                                                                                                                                     
gCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGgc                                                                                                                                                    <  5:245399/148‑1 (MQ=255)
  aTTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCt                                                                                                                                                 >  5:367513/1‑149 (MQ=255)
  aTTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCt                                                                                                                                                 <  1:317148/149‑1 (MQ=255)
      aGGCGGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑Gcc                                                                                                                                            <  7:250049/149‑1 (MQ=255)
         cTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAg                                                                                                                                         <  1:118614/149‑1 (MQ=255)
            cgcTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCa                                                                                                                                       <  6:134899/148‑1 (MQ=255)
            cgattGATCTGTGTGAGGGCGTAGGTGGGCTGGGCGGCGAGGTGGCGTGCGAGGGTGAGGGCTTTATCACGCATCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGc                                                                                                                                        <  8:367776/144‑1 (MQ=255)
                 gATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCa                                                                                                                                 >  2:198563/1‑149 (MQ=255)
                  aTCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAg                                                                                                                                <  6:350976/149‑1 (MQ=255)
                    cAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAgcg                                                                                                                              >  2:33141/1‑149 (MQ=255)
                      gtgtCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTG‑GCCATGCTGGGCGAGCGcc                                                                                                                            <  4:348088/149‑1 (MQ=255)
                       tgtCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGTGTGAGGTCCTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTG‑GCCATGCTGGGCGAGCGCCt                                                                                                                           <  8:317890/149‑1 (MQ=255)
                       tgtCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTG‑GCCATGCTGGGCGAGCGCCt                                                                                                                           <  5:127026/149‑1 (MQ=255)
                         tCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTg                                                                                                                         >  1:224923/1‑149 (MQ=255)
                         tCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTg                                                                                                                         <  4:398796/149‑1 (MQ=255)
                           aGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGc                                                                                                                       >  1:224/1‑149 (MQ=255)
                             gCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGccc                                                                                                                     >  5:321351/1‑149 (MQ=255)
                              ccGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAg                                                                                                                                <  4:70742/137‑1 (MQ=255)
                              ccGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAg                                                                                                                                >  3:70742/1‑137 (MQ=255)
                               cGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCgg                                                                                                                   >  7:132588/1‑149 (MQ=255)
                                     tgggctgggcGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATg                                                                                                             <  1:68378/149‑1 (MQ=255)
                                        gcggggcGGCGATGTGGCGGGCGCGGGGGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCCGCCTGTTCGGCGCCAACGCGCTCTGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGcc                                                                                                           <  1:270523/146‑1 (MQ=255)
                                            ggcggcGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGAcc                                                                                                       <  8:132588/148‑1 (MQ=255)
                                                   aGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCgg                                                                                               <  2:224923/149‑1 (MQ=255)
                                                        gcgggcgAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGgcagca                                                                                          >  4:372023/1‑149 (MQ=255)
                                                          gggcgAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCagcagc                                                                                        >  5:235327/1‑149 (MQ=255)
                                                                     aGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCTCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTAccgcc                                                                             <  1:198563/149‑1 (MQ=255)
                                                                      gggCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTAccgccg                                                                            <  6:367513/149‑1 (MQ=255)
                                                                             atcaCGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCg                                                                     >  8:146836/1‑149 (MQ=255)
                                                                              tcaCGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCAt                                                                                                              <  2:4695/107‑1 (MQ=255)
                                                                              tcaCGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCAt                                                                                                              >  1:4695/1‑107 (MQ=255)
                                                                              tcaCGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCgg                                                                    <  3:372023/149‑1 (MQ=255)
                                                                              tcaCGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCgg                                                                    >  5:334301/1‑149 (MQ=255)
                                                                                     gcgcTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAgg                                                             <  2:224/149‑1 (MQ=255)
                                                                                           catcatCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATc                                                       <  6:321351/149‑1 (MQ=255)
                                                                                           catcatCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATc                                                       >  2:388359/1‑149 (MQ=255)
                                                                                               atcCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTg                                                    <  1:1723/148‑1 (MQ=255)
                                                                                               atcCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTg                                                    <  7:146836/148‑1 (MQ=255)
                                                                                                   accacGCGGTAGATCAGCCCCCATTGCTTCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAg                                                 >  1:229384/1‑148 (MQ=255)
                                                                                                          ggTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTg                                         >  6:225317/1‑148 (MQ=255)
                                                                                                              gATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGGACCAGGCCGATCTTGCAGAAGGCCTGGATg                                     <  1:33141/148‑1 (MQ=255)
                                                                                                               aTCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGa                                    >  8:381976/1‑148 (MQ=255)
                                                                                                                  aGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCgg                                                                                               >  4:171338/1‑86 (MQ=255)
                                                                                                                  aGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCgg                                                                                               <  3:171338/86‑1 (MQ=255)
                                                                                                                    cccccATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTgg                              <  6:235327/149‑1 (MQ=255)
                                                                                                                     ccccATTGC‑TCGGCCTGTTCGGCGCCAGGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGgcagca                                                                                          >  4:9238/1‑88 (MQ=255)
                                                                                                                     ccccATTGC‑TCGGCCTGTTCGGCGCCAGGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGgcagca                                                                                          <  3:15016/88‑1 (MQ=255)
                                                                                                                     ccccATTGC‑TCGGCCTGTTCGGCGCCAGGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGgcagca                                                                                          >  4:15016/1‑88 (MQ=255)
                                                                                                                     ccccATTGC‑TCGGCCTGTTCGGCGCCAGGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGgcagca                                                                                          <  3:9238/88‑1 (MQ=255)
                                                                                                                     ccccATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGc                             <  6:334301/149‑1 (MQ=255)
                                                                                                                      cccATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGcc                            <  3:347947/149‑1 (MQ=255)
                                                                                                                       ccATAGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCATGGCCCGGGCAATGCCCACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCTCAAGGCCTGGATGAAACTTGCCg                           >  7:360470/1‑149 (MQ=255)
                                                                                                                             c‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAgcgg                      <  1:388359/148‑1 (MQ=255)
                                                                                                                                cGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCg                                                          >  2:77334/1‑110 (MQ=255)
                                                                                                                                cGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCg                                                          <  1:77334/110‑1 (MQ=255)
                                                                                                                                    cTGTTCGGCGCCAACGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGGCa               <  3:384428/149‑1 (MQ=255)
                                                                                                                                         cGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTgg                              >  2:55464/1‑129 (MQ=255)
                                                                                                                                         cGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTgg                              <  1:55464/129‑1 (MQ=255)
                                                                                                                                         cGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCAc           <  7:381976/148‑1 (MQ=255)
                                                                                                                                               cAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCg    >  4:248676/1‑149 (MQ=255)
                                                                                                                                                aaGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCg    >  7:106236/1‑148 (MQ=255)
                                                                                                                                                 aGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCa  <  5:225317/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GCATTGAGGCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGC‑TCGGCCTGTTCGGCGCCAAGGCGCTC‑GCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCA  >  NC_002947/3717937‑3718229

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