Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,318,991 T→A 8.5% H128L (CAC→CTC)  ilvI ← acetohydroxybutanoate synthase/acetolactate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,318,9910TA 8.5% 178.6 / 3.8 62H128L (CAC→CTC) ilvIacetohydroxybutanoate synthase/acetolactate synthase
Reads supporting (aligned to +/- strand):  ref base T (31/25);  new base A (2/3);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTATTCGAACTTCTCAGCCGGGTTGGTCATATCTTTTGGAATGTCGACCACGACTGGACCTGGACGACCGGATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGATGACCATCGGAATCGAGTCCATGTAGGCGGTGGCAATACCGGTAATGGCATTGGTC  >  NC_002947/5318852‑5319137
                                                                                                                                           |                                                                                                                                                  
gTATTCGAACTTCTCAGACGGGTTGGTAATATCTTTTGGAATGTCGACCACGACTGGACCTGGACGACCGGATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGCGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGa                                                                                                                                           >  5:296023/1‑149 (MQ=255)
     cGAACTTCTCACCCGGGAAGGTCATATCTTTTGGAATGACGACCACGACTGGACCTGGACGACCGGATTGCGCCAGGTAGAAAGCTTTTTTCAGAAATTCAGGGATCACGCTAGGATTCTAGATCGTGAAGCACAGCTTCATGATc                                                                                                                                         >  2:53182/1‑146 (MQ=255)
      gAACTTCTCAGCCGGGTTGGTCATATCTTTTGGAATGTCGACCACGACTGGACCTGGACGACCGGATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCGTGAAGCACAGCTTCATGATcaag                                                                                                                                      <  1:202389/148‑4 (MQ=255)
      gAACTTCTCAGCCGGGTTGGTCATATCTTTTGGAATGTCGACCACGACTGGACCTGGACGACCGGATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGcc                                                                                                                                     >  8:40576/1‑149 (MQ=255)
       aaCTTCTCAGCCGGGTTGGTCATATCTTTTGGAATGTCGACCACGACTGGACCTGGACGACCGGATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCGTGAAGCACAGCTTCATGATcaagaa                                                                                                                                    <  1:53182/149‑6 (MQ=255)
        aCTTCTCAGCCGGGTTGGTCATATCTTTTGGAATGTCGACCACGACTGGACCTGGACGACCGGATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCGTGAAGCACAGCTTCATGATcaagaat                                                                                                                                   <  7:311992/149‑7 (MQ=255)
            ctcAGCCGGGTTGGTCATATCTTTTGGAATGTCGACCACGACTGGACCTGGACGACCGGATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCcgc                                                                                                                                   >  4:231340/1‑145 (MQ=255)
               aGCCGGGTTGGTCATATCTTTTGGAATGTCGACCACGACTGGACCTGGACGACCGGATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCgaga                                                                                                                               <  3:148741/146‑1 (MQ=255)
               aGCCGGGTTGGTCATATCTTTTGGAATGTCGACCACGACTGGACCTGGACGACCGGATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATAc                                                                                                                            <  7:40576/149‑1 (MQ=255)
                          tCATATCTTTTGGAATGTCGACCACGACTGGACCTGGACGACCGGATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCg                                                                                                                 <  1:327709/149‑1 (MQ=255)
                          tCATATCTTTTGGAATGTCGACCACGACTGGACCTGGACGACCGGATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCg                                                                                                                 <  4:331151/149‑1 (MQ=255)
                              atCTTTTGGAATGTCGACCACGACTGGACCTGGACGACCGGATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACg                                                                                                                                            <  6:165215/118‑1 (MQ=255)
                              atCTTTTGGAATGTCGACCACGACTGGACCTGGACGACCGGATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACg                                                                                                                                            >  5:165215/1‑118 (MQ=255)
                                 ttttGGAATGTCGACCACGACTGGACCTGGACGACCGGATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCt                                                                                                          >  1:137332/1‑149 (MQ=255)
                                  tttGGAATGTCGACCACGACTGGACCTGGACGACCGGATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTg                                                                                                         >  2:250911/1‑149 (MQ=255)
                                      gAATGTCGACCACGACTGGACCTGGACGACCGGATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATc                                                                                                                       <  1:124905/131‑1 (MQ=255)
                                      gAATGTCGACCACGACTGGACCTGGACGACCGGATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATc                                                                                                                       >  2:124905/1‑131 (MQ=255)
                                          gTCGACCACGACTGGACCTGGACGACCGGATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCAt                                                                                                 >  4:183333/1‑149 (MQ=255)
                                             gACCACGACTGGACCTGGACGACCGGATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATa                                                                                                                             >  7:314558/1‑118 (MQ=255)
                                             gACCACGACTGGACCTGGACGACCGGATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATa                                                                                                                             <  8:314558/118‑1 (MQ=255)
                                             gACCACGACTGGACCTGGACGACCGGATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCg                                                                                               <  7:112218/148‑1 (MQ=255)
                                                cACGACTGGACCTCGACGACCGGATCGCGCCAGGTAGACAGCTTTTTTCCGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGCGCTTCACGATCGGCCGCGAGATACCAATCATGTCGCTTTCCTGGAAGGCATCGGTGc                                                                                           <  6:296023/149‑1 (MQ=255)
                                                    aCTGGACCTGGACGACCGGATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCccac                                                                                       <  3:183333/149‑1 (MQ=255)
                                                             ggacGACCGGATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGc                                                                               >  2:268758/1‑148 (MQ=255)
                                                              gacgacCGGATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGc                                                                               >  2:180400/1‑147 (MQ=255)
                                                              gacgacCGGATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTa                                                                             <  7:24211/149‑1 (MQ=255)
                                                                   ccGGATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAggacc                                                                        >  8:115987/1‑146 (MQ=255)
                                                                      gATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCt                                                                                                          <  6:32824/112‑1 (MQ=255)
                                                                      gATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCt                                                                                                          >  5:32824/1‑112 (MQ=255)
                                                                           cgcCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCg                                                                 >  4:141131/1‑148 (MQ=255)
                                                                                 gTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGa                                                          >  7:125233/1‑149 (MQ=255)
                                                                                   aGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGAt                                                         <  8:125233/148‑1 (MQ=255)
                                                                                    gAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGATGa                                                       <  8:33031/149‑1 (MQ=255)
                                                                                    gAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGATGa                                                       >  7:33031/1‑149 (MQ=255)
                                                                                     aaaGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGATGAc                                                      >  3:34741/1‑149 (MQ=255)
                                                                                        gCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACa                                                             >  4:265880/1‑139 (MQ=255)
                                                                                        gCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACa                                                             <  3:265880/139‑1 (MQ=255)
                                                                                         cTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGATGACCATc                                                  <  3:141131/149‑1 (MQ=255)
                                                                                          ttttttCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGATGACCATCg                                                 >  2:316289/1‑149 (MQ=255)
                                                                                          ttttttCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGATGACCATCg                                                 >  5:74208/1‑149 (MQ=255)
                                                                                              ttCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGATGACCATCGGaa                                              <  1:122337/148‑1 (MQ=255)
                                                                                               tCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGATGACCATCg                                                 <  7:108823/144‑1 (MQ=255)
                                                                                               tCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGATGACCATCg                                                 >  8:108823/1‑144 (MQ=255)
                                                                                                  gAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGATGACCATCGGAATCGAg                                         >  3:67272/1‑149 (MQ=255)
                                                                                                      ttCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGGACCTGGCCGGACAGGATGACCATCGGAATc                                            <  1:29695/142‑1 (MQ=255)
                                                                                                      ttCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGATGACCATCGGAATc                                            >  2:29695/1‑142 (MQ=255)
                                                                                                        cAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGATGACCATCGGAATCGAGTCCATg                                   >  3:50973/1‑149 (MQ=255)
                                                                                                         aGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGATGACCATCGGAATCGAGTCCATg                                   >  5:255124/1‑148 (MQ=255)
                                                                                                          gggATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGATGACCATCGGAATCGAGTCCATGTa                                 <  1:268758/149‑1 (MQ=255)
                                                                                                           ggATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGATGACCATCGGAATCGAGTCCATGTAg                                <  1:250911/149‑1 (MQ=255)
                                                                                                                  ggtcggATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATg                                                                                   >  7:171879/1‑91 (MQ=255)
                                                                                                                  ggtcggATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATg                                                                                   <  8:171879/91‑1 (MQ=255)
                                                                                                                  ggtcggATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGATGACCATCGGAATCGAGTCCATGTAGGCgg                            <  1:180400/146‑1 (MQ=255)
                                                                                                                  ggtcggATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGATGACCATCGGAATCGAGTCCATGTAGGCGGTgg                         >  6:63791/1‑149 (MQ=255)
                                                                                                                  ggtcggATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGATGACCATCGGAATCGAGTCCATGTAGGCGGTg                          <  1:316289/148‑1 (MQ=255)
                                                                                                                   gtcggATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGATGACCATCGGAATCGAGTCCATGTAGGCGGTgg                         >  4:9494/1‑148 (MQ=255)
                                                                                                                     cggATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGATGACCATCGGAATCGAGTCCATGTAGGCGGTGGCa                       >  7:211355/1‑148 (MQ=255)
                                                                                                                        aTTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGATGACCATCGGAATCGAGTCCATGTAGGCGGTGGCAATAc                   >  5:335206/1‑149 (MQ=255)
                                                                                                                            ttGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGATGACCATCGGAATCGAGTCCATGTAGGCGGTGGCAATACCgg                >  7:217747/1‑148 (MQ=255)
                                                                                                                              gATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGATGACCATCGGAATCGAGTCCATGTAGGCGGTGGCAATACCgg                <  8:217747/146‑1 (MQ=255)
                                                                                                                               aTCGTGAAGCACAGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGATGACCATCGGAATCGAGTCCATGTAGGCGGTGGCAATACCGGTAAt            >  3:98508/1‑149 (MQ=255)
                                                                                                                               aTCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGATGACCATCGGAATCGAGTCCATGTAGGCGGTGGCAATACCGGTaa             <  4:303653/148‑1 (MQ=255)
                                                                                                                                         tgtgCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGATGACCATCGGAATCGAGTCCATGTAGGCGGTGGCAATACCGGTAATGGCATTGGTc  >  3:60788/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                  
GTATTCGAACTTCTCAGCCGGGTTGGTCATATCTTTTGGAATGTCGACCACGACTGGACCTGGACGACCGGATTGCGCCAGGTAGAAAGCTTTTTTCAGAACTTCAGGGATCTCGGTCGGATTCTTGATCATGAAGCTGTGCTTCACGATCGGCCGCGAGATACCAATCATGTCGGTTTCCTGGAAGGCATCGGTGCCCACCATGGTGCTAGGCACCTGGCCGGACAGGATGACCATCGGAATCGAGTCCATGTAGGCGGTGGCAATACCGGTAATGGCATTGGTC  >  NC_002947/5318852‑5319137

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