Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,572,889 C→A 7.6% A337A (GCC→GCA)  pcaB → 3‑carboxy‑cis,cis‑muconate cycloisomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,572,8890CA 7.6% 209.3 / 6.1 66A337A (GCC→GCA) pcaB3‑carboxy‑cis,cis‑muconate cycloisomerase
Reads supporting (aligned to +/- strand):  ref base C (30/31);  new base A (3/2);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

CCGTGTTGATCGGTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGTCGCGACCGTGC  >  NC_002947/1572744‑1573029
                                                                                                                                                 |                                                                                                                                            
ccgTGTTGATCGGTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTg                                                                                                                                           >  3:28304/1‑149 (MQ=255)
 cgTGTTGATCGGTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTgc                                                                                                                                          >  3:40772/1‑149 (MQ=255)
     ttGATCGGTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTgcgc                                                                                                                                        >  8:321090/1‑147 (MQ=255)
     ttGATCGGTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCa                                                                                                                                      <  2:260547/149‑1 (MQ=255)
       gATCGGTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAg                                                                                                                                     <  6:27360/148‑1 (MQ=255)
        aTCGGTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAgg                                                                                                                                    <  3:248953/148‑1 (MQ=255)
        aTCGGTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCAGGGCGCCAga                                                                                                                                    <  5:165609/148‑2 (MQ=255)
          cGGTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGGACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCAGGGCGCCAgaga                                                                                                                                  <  5:240480/148‑4 (MQ=255)
           ggTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTc                                                                                                                                 >  4:309378/1‑148 (MQ=255)
           ggTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTc                                                                                                                                 >  8:238466/1‑148 (MQ=255)
           ggTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGGTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCGGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATCTCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCt                                                                                                                                  <  5:253357/147‑1 (MQ=255)
            gTGCCGCGACCCGCGTGTCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCaa                                                                                                                               >  3:32862/1‑149 (MQ=255)
             tGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAg                                                                                                                              <  2:97449/149‑1 (MQ=255)
               ccGCGACCCGCGTGCCGGGCCTGCTGTCGACACTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCAGGTCTCTGGCGCCCTGCGCCAGGCTCAAGt                                                                                                                             >  1:43501/1‑148 (MQ=255)
                       cgcgTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGcc                                                                                                                      >  1:11984/1‑147 (MQ=255)
                        gcgTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAg                                                                                                                   <  4:55983/149‑1 (MQ=255)
                              cGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATg                                                                                                             <  7:137585/149‑1 (MQ=255)
                              cGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATg                                                                                                             >  4:155804/1‑149 (MQ=255)
                               gggCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGcc                                                                                                                                              <  4:139263/115‑1 (MQ=255)
                               gggCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGcc                                                                                                                                              >  3:139263/1‑115 (MQ=255)
                                 gCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAg                                                                                                          <  1:391194/149‑1 (MQ=255)
                                      ctgTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTgg                                                                                                      <  1:381744/148‑1 (MQ=255)
                                         tCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGc                                                                                                  >  5:81081/1‑149 (MQ=255)
                                           gACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCcg                                                                                                >  2:181444/1‑149 (MQ=255)
                                              gCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCcgcgc                                                                                             <  7:321090/149‑1 (MQ=255)
                                              gCTGTTCGCAGCCATGCCTCACGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCcgcgc                                                                                             <  3:309378/149‑1 (MQ=255)
                                                tGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGcc                                                                                                                      <  7:85558/122‑1 (MQ=255)
                                                tGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGcc                                                                                                                      >  8:85558/1‑122 (MQ=255)
                                                       aGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGTTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGcc                                                                                    >  7:81006/1‑149 (MQ=255)
                                                          cATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGt                                                                                  >  3:249797/1‑148 (MQ=255)
                                                           aTGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTaa                                                                                <  4:40772/149‑1 (MQ=255)
                                                            tGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAAc                                                                               <  6:6619/149‑1 (MQ=255)
                                                               cTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTc                                                                            <  4:32862/149‑1 (MQ=255)
                                                                tCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTc                                                                            >  8:325854/1‑148 (MQ=255)
                                                                 cAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGa                                                                          >  1:250259/1‑149 (MQ=255)
                                                                  aGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGAc                                                                         >  7:127478/1‑149 (MQ=255)
                                                                  aGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGAc                                                                         <  3:357936/149‑1 (MQ=255)
                                                                      gCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGa                                                                     <  7:238466/149‑1 (MQ=255)
                                                                      gCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGa                                                                     >  2:298348/1‑149 (MQ=255)
                                                                        aCGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGAcc                                                                   >  4:226730/1‑149 (MQ=255)
                                                                           aaCGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGAccc                                                                  >  5:180520/1‑147 (MQ=255)
                                                                                 gCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTg                                                          <  1:298348/149‑1 (MQ=255)
                                                                                 gCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTg                                                          <  2:303005/149‑1 (MQ=255)
                                                                                 gCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTg                                                          <  6:81081/149‑1 (MQ=255)
                                                                                  ccTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTgg                                                         <  5:318912/149‑1 (MQ=255)
                                                                                   cTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGTTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGt                                                        <  8:81006/149‑1 (MQ=255)
                                                                                       gCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTg                                                    >  4:76634/1‑149 (MQ=255)
                                                                                        ccTATGGCATGCCGAGTGGGAAACCCTTCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTgg                                                   <  4:28304/149‑1 (MQ=255)
                                                                                              gCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGaa                                              <  1:181444/148‑1 (MQ=255)
                                                                                                    cGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGt                                         >  8:143535/1‑147 (MQ=255)
                                                                                                    cGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGa                                       <  7:325854/149‑1 (MQ=255)
                                                                                                          ggAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATc                                  <  3:76634/148‑1 (MQ=255)
                                                                                                          ggAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCg                                 >  3:10350/1‑149 (MQ=255)
                                                                                                            aaaCCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTc                               <  2:11984/149‑1 (MQ=255)
                                                                                                                ccTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCg                           >  3:369721/1‑149 (MQ=255)
                                                                                                                ccTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCg                           <  8:127478/149‑1 (MQ=255)
                                                                                                                  tGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGcc                         >  6:17320/1‑149 (MQ=255)
                                                                                                                  tGCCGGATATCTGCTGCCTGGTCTCTGGCGCCATGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATTCCGCCCACATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCATAGcc                         >  7:206991/1‑149 (MQ=255)
                                                                                                                        atatCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGt                   >  8:109140/1‑149 (MQ=255)
                                                                                                                          atCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCt                 >  8:281105/1‑149 (MQ=255)
                                                                                                                            ctgctgCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTgg               >  2:263875/1‑149 (MQ=255)
                                                                                                                               ctgcCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTggg              <  3:155804/147‑1 (MQ=255)
                                                                                                                                   cTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGTCGCGa        >  4:139568/1‑149 (MQ=255)
                                                                                                                                      gTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGTCGCGAc       <  6:180520/147‑1 (MQ=255)
                                                                                                                                         cctGGCGCAGGGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGTCGCGACCGTGc  >  5:204048/2‑149 (MQ=255)
                                                                                                                                         cctGGCGCAGGGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGTCGCGACCGTGc  >  7:165244/2‑149 (MQ=255)
                                                                                                                                         cctGGCGCAGGGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGCGTCGCCACCGTGc  >  7:68501/2‑149 (MQ=255)
                                                                                                                                           tGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGTCGCGACCGTg   <  2:250259/146‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
CCGTGTTGATCGGTGCCGCGACCCGCGTGCCGGGCCTGCTGTCGACGCTGTTCGCAGCCATGCCTCAGGAGCACGAACGCAGCCTGGGCCTATGGCATGCCGAGTGGGAAACCCTGCCGGATATCTGCTGCCTGGTCTCTGGCGCCCTGCGCCAGGCTCAAGTGATTGCCGAGGGCATGGAGGTGGATGCCGCGCGCATGCGCCGTAACCTCGACCTGACCCAAGGCCTGGTGCTGGCCGAAGCGGTGAGCATCGTCCTCGCCCAGCGTCTGGGTCGCGACCGTGC  >  NC_002947/1572744‑1573029

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