Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,462,690 G→A 7.0% D556D (GAC→GAT)  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,6900GA 7.0% 164.9 / 2.1 57D556D (GAC→GAT) PP_3078ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (24/29);  new base A (2/2);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.78e-01

AGCGGTTCCACCAGACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCGGCGACAGGGTGACCGGGTAGCCGGTGACGATCATGTCGTAGTCGCGCGCCATCAAC  >  NC_002947/3462545‑3462837
                                                                                                                                                  |                                                                                                                                                   
aGCGGTTCCACCAGACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCt                                                                                                                                                  <  6:149326/149‑1 (MQ=255)
   ggTTCCACCAGACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACg                                                                                                                                               <  4:351192/149‑1 (MQ=255)
   ggTTCCACCAGACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATAGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACAAGGCCGTCGATCAGGTGGTCTCCg                                                                                                                                               >  6:141874/1‑149 (MQ=255)
     ttCCACCAGACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGc                                                                                                                                             >  3:93591/1‑149 (MQ=255)
        cacaaGAAGGTGGAGATGCCCGGTGGGTAATAGTTGGGGATCAAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCgg                                                                                                                                          <  1:297576/149‑1 (MQ=255)
         accaGACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCggg                                                                                                                                         >  4:56984/1‑149 (MQ=255)
           caGACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTc                                                                                                                                       >  4:208033/1‑149 (MQ=255)
           caGACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTc                                                                                                                                       <  7:227476/149‑1 (MQ=255)
           caGACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTAGACCACCTGATCGACGGcctg                                                                                                                                          <  1:118801/146‑3 (MQ=255)
           caGACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTAGACCACCTGATCGACGGcctg                                                                                                                                          >  2:118801/1‑144 (MQ=255)
           caGACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGAGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCCGGCCGTAGA‑CACCTGATCGACGGCCTGGTcc                                                                                                                                      <  4:75156/149‑1 (MQ=255)
             gACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGGCGGC‑CGGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCt                                                                                                                                     <  4:258359/149‑1 (MQ=255)
                     aGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGGCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAAcc                                                                                                                             <  8:305275/149‑1 (MQ=255)
                          cccGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGCCCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAg                                                                                                                        >  4:42443/1‑149 (MQ=255)
                             ggTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGgtt                                                                                                                     >  2:81709/1‑149 (MQ=255)
                             ggTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGgtt                                                                                                                     <  5:57634/149‑1 (MQ=255)
                                 ggTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTc                                                                                                                 <  3:208033/149‑1 (MQ=255)
                                 ggTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTc                                                                                                                 <  1:378650/149‑1 (MQ=255)
                                   tAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGTGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTc                                                                                                                 >  1:281270/1‑147 (MQ=255)
                                   tAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCTAAACCATCAGGTTGTTc                                                                                                                 >  7:316916/1‑147 (MQ=255)
                                    aaTAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAg                                                                                                              >  2:279353/1‑149 (MQ=255)
                                       aGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCg                                                                                                           >  3:59337/1‑149 (MQ=255)
                                        gTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATTAGGTTGTTCGAGCCgg                                                                                                          <  7:310189/149‑1 (MQ=255)
                                        gTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATTAGGTTGTTCGAGCCgg                                                                                                          >  8:310189/1‑149 (MQ=255)
                                            gggATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTc                                                                                                      <  2:204560/149‑1 (MQ=255)
                                            gggATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTc                                                                                                      <  2:74898/149‑1 (MQ=255)
                                                  cAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGc                                                                                                <  7:288409/149‑1 (MQ=255)
                                                   ggtagtaATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCg                                                                                               >  2:352853/2‑149 (MQ=255)
                                                   agtagtAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCg                                                                                               >  6:214169/1‑149 (MQ=255)
                                                    gtagtaATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCgg                                                                                              >  1:384881/1‑149 (MQ=255)
                                                          aTTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCAGTGGGTGTCGGCCCGGACCAGGCCGACCAACAGGTGGAATAGGCCCGGGACATGCAAACCCATCAGGTTGTGCGAGCCGGGGACGTGAGCGgccgcc                                                                                         >  1:95832/1‑149 (MQ=255)
                                                             ccATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGACGTGAGCGGCCGCCGAAc                                                                                     >  1:359406/1‑149 (MQ=255)
                                                                  gcagcaCGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAag                                                                                >  5:81355/1‑149 (MQ=255)
                                                                       aggcgGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑GTGGGTGTCGGCGCGGACCAGGCGGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTTAGCGGCCGCCGAACCAAAGTAGtt                                                                           <  5:141874/147‑1 (MQ=255)
                                                                         gcgGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTg                                                                          >  6:141768/1‑148 (MQ=255)
                                                                           ggTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTAc                                                                       <  3:56984/149‑1 (MQ=255)
                                                                                 gggCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCg                                                                 >  3:280126/1‑149 (MQ=255)
                                                                                 gggCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCg                                                                 <  2:384881/149‑1 (MQ=255)
                                                                                  ggCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCgg                                                                <  4:59337/149‑1 (MQ=255)
                                                                                         gCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCgg                                                         <  1:279353/149‑1 (MQ=255)
                                                                                         gCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCgg                                                         >  4:275680/1‑149 (MQ=255)
                                                                                          cGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGCCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCGGc                                                        <  2:359406/149‑1 (MQ=255)
                                                                                              gcgcAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCGGCGACa                                                    <  4:93591/149‑1 (MQ=255)
                                                                                                gcagcaTGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGcc                                                                                                                                            >  8:304040/1‑59 (MQ=255)
                                                                                                gcagcaTGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGcc                                                                                                                                            <  7:304040/59‑1 (MQ=255)
                                                                                                gcagcaTGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCGGCGACAgg                                                  <  2:281270/149‑1 (MQ=255)
                                                                                                     aTGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCGGCGACAGGGTGAc                                             <  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)
                                                                                                                                           tCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCGGCGACAGGGTGACCGGGTAGCCGGTGACGATCATGTCGTAGTCgcgcgc         >  3:187423/1‑148 (MQ=255)
                                                                                                                                           tCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCGGCGACAGGGTGACCGGGTAGCCGGTGACGATCATGTCGTAGTCGCGCGcc        <  4:146962/149‑1 (MQ=255)
                                                                                                                                               gtggtCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCGGCGACAGGGTGACCGGGTAGCCGGTGACGATCATGTCGTAGTCGCGCGCCATc     >  3:191759/1‑148 (MQ=255)
                                                                                                                                                  gtCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCGGCGACAGGGTGACCGGGTAGCCGGTGACGATCATGTCGTAGTCGCGCGCCATCAAc  <  5:214169/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AGCGGTTCCACCAGACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGC‑CTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCGGCGACAGGGTGACCGGGTAGCCGGTGACGATCATGTCGTAGTCGCGCGCCATCAAC  >  NC_002947/3462545‑3462837

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