Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,002,505 T→C 6.8% T478A (ACC→GCC)  amaB ← L‑piperidine‑6‑carboxylate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,002,5050TC 6.8% 165.1 / 2.7 59T478A (ACC→GCC) amaBL‑piperidine‑6‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (31/24);  new base C (2/2);  total (33/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 = 8.86e-01

CAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCCTTCTCGCCACCAAACGCACCACCAATCTCGGCGCCGCTGGTACCGATGTTGACGTTGGCGATGCCGCAGTCGCTGCCCGAGGCGCT  >  NC_002947/6002357‑6002652
                                                                                                                                                    |                                                                                                                                                   
cAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGt                                                                                                                                                     <  5:339022/149‑1 (MQ=255)
   gggTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTg                                                                                                                                                  >  7:3454/1‑149 (MQ=255)
   gggTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGCCTg                                                                                                                                                  >  7:371499/1‑149 (MQ=255)
    ggTCGTGGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTgg                                                                                                                                                 <  8:156669/149‑1 (MQ=255)
    ggTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTgg                                                                                                                                                 >  6:19540/1‑149 (MQ=255)
        ggTGTCGCTTGGACGAACCCTGTAGGAACGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGAc                                                                                                                                             <  8:178052/149‑1 (MQ=255)
              gCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATg                                                                                                                                        >  3:375223/1‑148 (MQ=255)
               cTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTa                                                                                                                                      <  8:24242/149‑1 (MQ=255)
                 gtgACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAg                                                                                                                                     <  6:121089/146‑1 (MQ=255)
                   gACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGcc                                                                                                                                   <  2:364555/148‑1 (MQ=255)
                    aCGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCtt                                                                                                                                 <  5:84963/149‑1 (MQ=255)
                      gAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTc                                                                                                                               >  1:208488/1‑149 (MQ=255)
                        aCCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTc                                                                                                                               >  8:217506/1‑147 (MQ=255)
                        aCCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCa                                                                                                                             <  8:19602/149‑1 (MQ=255)
                        aCCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCa                                                                                                                             <  6:309674/149‑1 (MQ=255)
                         cccTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAg                                                                                                                            <  2:208488/149‑1 (MQ=255)
                               agcaGCGCGTCTACCCGAGAAACGGGGGCCAGCACCCCGGCCAAGCAGCAACAATACATCAGTCGAACACGAAACCATGCCCGAGCACCAGCACCCGCGAGTAGTAAACGGTGTAGCTCTGGAGACCCAAGAAGCCTTTACAGGGATcc                                                                                                                      >  4:325695/4‑149 (MQ=255)
                                  aGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCCAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAg                                                                                                                   <  2:348235/149‑1 (MQ=255)
                                           aCCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCa                                                                                                          >  1:161029/1‑149 (MQ=255)
                                            cccGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTg                                                                                                                                                  <  4:381687/108‑1 (MQ=255)
                                            cccGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTg                                                                                                                                                  >  3:381687/1‑108 (MQ=255)
                                              cgcgAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACACATGTAGCCTTTCCAGGCATCCGAGCCCGACTc                                                                                                           >  2:250321/1‑145 (MQ=255)
                                              cgcgAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACACATGTAGCCTTTCCAGGCATCCGAGCCCGACTc                                                                                                           <  1:250321/145‑1 (MQ=255)
                                               gcgAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGc                                                                                                      >  5:205753/1‑149 (MQ=255)
                                               gcgAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGc                                                                                                      >  2:351051/1‑149 (MQ=255)
                                                    ccGGGGGCCTGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCTAGCCCGACTCACGGCCAccg                                                                                                 >  6:373243/1‑149 (MQ=255)
                                                     cGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGt                                                                                                                                                     >  4:284750/1‑96 (MQ=255)
                                                     cGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGt                                                                                                                                                     <  3:284750/96‑1 (MQ=255)
                                                          gCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGtt                                                                                           <  8:3454/149‑1 (MQ=255)
                                                            cAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTc                                                                                         <  1:28780/149‑1 (MQ=255)
                                                            cAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTc                                                                                         <  6:205753/149‑1 (MQ=255)
                                                            cAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTc                                                                                         <  5:19540/149‑1 (MQ=255)
                                                             aGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTcc                                                                                        >  5:330499/1‑149 (MQ=255)
                                                                 ccccGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCCTtct                                                                                    >  3:84997/1‑149 (MQ=255)
                                                                 ccccGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGCCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCCTtct                                                                                    <  8:371499/149‑1 (MQ=255)
                                                                  cccGGCCATGCAGCATCAATACCTCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCCTtctc                                                                                   <  5:306266/149‑1 (MQ=255)
                                                                  cccGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCCTtctc                                                                                   <  1:252843/149‑1 (MQ=255)
                                                                                aTCAATAAATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCAATCAAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCCTTCTCGCCACCACACGcac                                                                     >  6:266277/1‑149 (MQ=255)
                                                                                     tACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCCTTCTCGCCACCAAACGCaccacc                                                                 >  2:131425/1‑148 (MQ=255)
                                                                                       cATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCCTTCTCGCCACCAAACGCACCACCaa                                                               <  4:84997/148‑1 (MQ=255)
                                                                                           aGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTGGTTCACGGTGTTGGCCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCCTTCTCGcca                                                                               >  1:272881/1‑128 (MQ=255)
                                                                                           aGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGCCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCCTTCTCGcca                                                                               <  2:272881/128‑1 (MQ=255)
                                                                                             tCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCCTTCTCGCCACCAAACGCACCACCAATCTCGgc                                                        <  2:161029/149‑1 (MQ=255)
                                                                                                aaCACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCCTTCTCGCCACCAAACGCACCACCAATCTCGGcgcc                                                     >  7:301641/1‑149 (MQ=255)
                                                                                                   acGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCCTTCTCGCCACCAAACGCACCACCAATCTCGGCGCCGCt                                                  >  3:275320/1‑149 (MQ=255)
                                                                                                    cGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCCTTCTCGCCACCAAACGCACCACCAATCTCGGCGCCGCTg                                                 <  7:217506/149‑1 (MQ=255)
                                                                                                     gATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCCTTCTCGCCACCAAACGCACCACCAATCTCGGCGCCGCTgg                                                >  8:49025/1‑149 (MQ=255)
                                                                                                     gATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCCTTCTCGCCACCAAACGCACCACCAATCTCGGCGCCGCTgg                                                >  8:309694/1‑149 (MQ=255)
                                                                                                          ccTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGCCGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCCTTCTCGCCACCAAATGCACCACCAATCTCGGCGCCGCTGGTAcc                                            >  6:239636/1‑148 (MQ=255)
                                                                                                          ccTGCGCCAGCGGCAGCTCGCGCGAGTAGTTAACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCAGAGCCCGACTCACGGCCACCGCCGGTTTCCTTCTCGCCACCAAACGCACCACCAATCTCGGCGCCGCTGGTAcc                                            >  6:224736/1‑148 (MQ=255)
                                                                                                               gcCAGCGGCAGCTCGCGCGAGTAGTTCACGGTCTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCCTTCTCGCCACCAAACGCACCACCAATCTCGGCGCCGCTGGTACCGATGtt                                      <  4:137399/149‑1 (MQ=255)
                                                                                                                ccAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTACCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTATCCTTCTCGCCACCAAACGCCCCACCCATATCGGCGCCGCTGGTACCGATGtt                                      >  7:44933/1‑148 (MQ=255)
                                                                                                                       cAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCCTTCTCGCCACCAAACGCACCACCAATCTCGGCGCCGCTGGTACCGATGTTGACGTTgg                              <  1:351051/149‑1 (MQ=255)
                                                                                                                        aGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCAGACTCACGGCCACCGCCGGTTTCCTTCTCGCCACCAAACGCACCACCAATCTCGGCGCCGCTGGTACCGATGTTGACGTTGGc                             >  8:196284/1‑149 (MQ=255)
                                                                                                                           tCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCCTTCTCGCCACCAAACGCACCACCAATCTCGGCGCCGCTGGTACCGATGTTGACGTTGGCGAt                          >  7:2453/1‑149 (MQ=255)
                                                                                                                           tCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCCTTCTCGCCACCAAACGCACAACCAATCTCGGCGCCGCTGGTACCGATGTTGACGTTGGCGAt                          >  3:385291/1‑149 (MQ=255)
                                                                                                                             gcgcgAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCCTTCTCGCCACCAAACGCACCACCAATCTCGGCGCCGCTGGTACCGATGTTGACGTTGGCGATGc                        >  8:92960/1‑149 (MQ=255)
                                                                                                                              cgcgAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCCTTCTCGCCACCAAACGCACCACCAATCTCGGCGCCGCTGGTACCGATGTTGACGTTGGCGATGcc                       <  1:291940/149‑1 (MQ=255)
                                                                                                                                 gAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCCTTCTCGCCACCAAACGCACCACCAATCTCGGCGCCGCTGGTACCGATGTTGACGTTGGCGATGCCg                      >  1:93983/1‑147 (MQ=255)
                                                                                                                                   gtagtTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCCTTCTCGCCACCAACCGCACCACCAATCTCGGCGCCGCTGGTACCGATGTTGACGTTGGCGATGACGCAGt                  >  2:378872/1‑149 (MQ=255)
                                                                                                                                            ggTGTTGGTCTGGCGACGCATGTAGCCTTTCAAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCCTTCTCGCCACCAAACGCACCACCAATCTCGGCGCAGCTGGTACCGATGTTGACGTTGTCGATGCCGCAGTCGCTGccc          >  1:43130/1‑148 (MQ=255)
                                                                                                                                                   gTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCCTTCTCGCCACCAAACGCACCACCAATCTCGGCGCCGCTGGTACCGATGTTGACGTTGGCGATGCCGCAGTCGCTGCCCGAGGCGct  >  6:142239/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCCTTCTCGCCACCAAACGCACCACCAATCTCGGCGCCGCTGGTACCGATGTTGACGTTGGCGATGCCGCAGTCGCTGCCCGAGGCGCT  >  NC_002947/6002357‑6002652

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: ‑