Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,462,680 C→A 6.7% D560Y (GAC→TAC)  PP_3078 ← ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,462,6800CA 6.7% 189.0 / 3.8 60D560Y (GAC→TAC) PP_3078ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (25/31);  new base A (2/2);  total (27/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 = 1.35e-01

GGCAGGCCGAAGCGGTTCCACCAGACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCGGCGACAGGGTGACCGGGTAGCCGGTGACGATCATGTCGTAGTCG  >  NC_002947/3462535‑3462825
                                                                                                                                                  |                                                                                                                                                 
ggCAGGCCGAAGCGGTTCCACCAGACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGa                                                                                                                                                 <  7:208724/148‑1 (MQ=255)
          aGCGGTTCCACCAGACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCt                                                                                                                                      <  6:149326/149‑1 (MQ=255)
             ggTTCCACCAGACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACg                                                                                                                                   <  4:351192/149‑1 (MQ=255)
             ggTTCCACCAGACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATAGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACAAGGCCGTCGATCAGGTGGTCTCCg                                                                                                                                   >  6:141874/1‑149 (MQ=255)
               ttCCACCAGACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGc                                                                                                                                 >  3:93591/1‑149 (MQ=255)
                  cacaaGAAGGTGGAGATGCCCGGTGGGTAATAGTTGGGGATCAAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCgg                                                                                                                              <  1:297576/149‑1 (MQ=255)
                   accaGACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCggg                                                                                                                             >  4:56984/1‑149 (MQ=255)
                     caGACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTc                                                                                                                           >  4:208033/1‑149 (MQ=255)
                     caGACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTc                                                                                                                           <  7:227476/149‑1 (MQ=255)
                     caGACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTAGACCACCTGATCGACGGcctg                                                                                                                              <  1:118801/146‑3 (MQ=255)
                     caGACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTAGACCACCTGATCGACGGcctg                                                                                                                              >  2:118801/1‑144 (MQ=255)
                     caGACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGAGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCCGGCCGTAGACACCTGATCGACGGCCTGGTcc                                                                                                                          <  4:75156/149‑1 (MQ=255)
                       gACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGGCGGCCGGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCt                                                                                                                         <  4:258359/149‑1 (MQ=255)
                               aGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGGCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAAcc                                                                                                                 <  8:305275/149‑1 (MQ=255)
                                    cccGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGCCCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAg                                                                                                            >  4:42443/1‑149 (MQ=255)
                                       ggTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGgtt                                                                                                         <  5:57634/149‑1 (MQ=255)
                                       ggTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGgtt                                                                                                         >  2:81709/1‑149 (MQ=255)
                                           ggTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTc                                                                                                     <  1:378650/149‑1 (MQ=255)
                                           ggTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTc                                                                                                     <  3:208033/149‑1 (MQ=255)
                                             tAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGTGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTc                                                                                                     >  1:281270/1‑147 (MQ=255)
                                             tAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAgg                                                                                                                                            <  6:200123/108‑1 (MQ=255)
                                             tAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAgg                                                                                                                                            >  5:200123/1‑108 (MQ=255)
                                             tAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCTAAACCATCAGGTTGTTc                                                                                                     >  7:316916/1‑147 (MQ=255)
                                              aaTAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAg                                                                                                  >  2:279353/1‑149 (MQ=255)
                                                 aGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCg                                                                                               >  3:59337/1‑149 (MQ=255)
                                                  gTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATTAGGTTGTTCGAGCCgg                                                                                              >  8:310189/1‑149 (MQ=255)
                                                  gTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATTAGGTTGTTCGAGCCgg                                                                                              <  7:310189/149‑1 (MQ=255)
                                                      gggATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTc                                                                                          <  2:204560/149‑1 (MQ=255)
                                                      gggATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTc                                                                                          <  2:74898/149‑1 (MQ=255)
                                                            cAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGc                                                                                    <  7:288409/149‑1 (MQ=255)
                                                             ggtagtaATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCg                                                                                   >  2:352853/2‑149 (MQ=255)
                                                             agtagtAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAggtg                                                                                                                                          <  4:130226/94‑1 (MQ=255)
                                                             agtagtAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAggtg                                                                                                                                          >  3:130226/1‑94 (MQ=255)
                                                             agtagtAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCg                                                                                   >  6:214169/1‑149 (MQ=255)
                                                              gtagtaATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCgg                                                                                  >  1:384881/1‑149 (MQ=255)
                                                                    aTTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCAGTGGGTGTCGGCCCGGACCAGGCCGACCAACAGGTGGAATAGGCCCGGGACATGCAAACCCATCAGGTTGTGCGAGCCGGGGACGTGAGCGgccgcc                                                                             >  1:95832/1‑149 (MQ=255)
                                                                       ccATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGACGTGAGCGGCCGCCGAAc                                                                         >  1:359406/1‑149 (MQ=255)
                                                                            gcagcaCGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAag                                                                    >  5:81355/1‑149 (MQ=255)
                                                                                 aggcgGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCGTGGGTGTCGGCGCGGACCAGGCGGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTTAGCGGCCGCCGAACCAAAGTAGtt                                                               <  5:141874/147‑1 (MQ=255)
                                                                                   gcgGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTg                                                              >  6:141768/1‑148 (MQ=255)
                                                                                     ggTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTAc                                                           <  3:56984/149‑1 (MQ=255)
                                                                                           gggCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCg                                                     <  2:384881/149‑1 (MQ=255)
                                                                                           gggCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCg                                                     >  3:280126/1‑149 (MQ=255)
                                                                                            ggCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCgg                                                    <  4:59337/149‑1 (MQ=255)
                                                                                                   gCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCgg                                             <  1:279353/149‑1 (MQ=255)
                                                                                                   gCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCgg                                             >  4:275680/1‑149 (MQ=255)
                                                                                                    cGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGCCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCGGc                                            <  2:359406/149‑1 (MQ=255)
                                                                                                        gcgcAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCGGCGACa                                        <  4:93591/149‑1 (MQ=255)
                                                                                                          gcagcaTGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGcc                                                                                                                                <  7:304040/59‑1 (MQ=255)
                                                                                                          gcagcaTGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGcc                                                                                                                                >  8:304040/1‑59 (MQ=255)
                                                                                                          gcagcaTGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCGGCGACAgg                                      <  2:281270/149‑1 (MQ=255)
                                                                                                               aTGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCGGCGACAGGGTGAc                                 <  1:352853/149‑1 (MQ=255)
                                                                                                                           ggTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTTGTTGTACAACTCGGCGCCCGGCGACAGGGTGACCGGGTAGCCgg                      <  1:81709/149‑1 (MQ=255)
                                                                                                                           ggTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCTGCGCCCGGCGACAGGGTGACCGGGTAGCCgg                      <  6:81355/149‑1 (MQ=255)
                                                                                                                                 ggCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCGGCGACAGGGTGACCGGGTAGCCGGTGACg                 <  2:228161/148‑1 (MQ=255)
                                                                                                                                  gcgcGGACCAGGCCGTCGATCAGGTGGGCTACGGCCGGGTCCTGCAAACCCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCGGCGACAGGGGGACCGGGTAGCCGGTGACGAt               <  7:280549/149‑1 (MQ=255)
                                                                                                                                         cccggCCGTAGACCACCTGATCGACGGCCTGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCGGCGACAGGGTGACCGGGTAGCCGGTGACGATCATGTCg        >  6:151132/4‑149 (MQ=38)
                                                                                                                                          cAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCGGCGACAGGGTGACCGGGTAGCCGGTGACGATCATGTCGt       <  3:275680/149‑1 (MQ=255)
                                                                                                                                               cGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCGGCGACAGGGTGACCGGGTAGCCGGTGACGATCATGTCGTAGTcg  >  3:146962/1‑149 (MQ=255)
                                                                                                                                               cGTCGATCAGGTGGTCTACGGCCGCGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCGGCGACAGGGTGACCGGGTAGCCGGTGACGATCATGTCGTAGTcg  >  8:123831/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GGCAGGCCGAAGCGGTTCCACCAGACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCGGCGACAGGGTGACCGGGTAGCCGGTGACGATCATGTCGTAGTCG  >  NC_002947/3462535‑3462825

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: