Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,462,687 G→C 7.0% H557Q (CAC→CAG 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,6870GC 7.0% 179.6 / 3.2 57H557Q (CAC→CAGPP_3078ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (24/29);  new base C (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 = 9.59e-01

AGCGGTTCCACCAGACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCGGCGACAGGGTGACCGGGTAGCCGGTGACGATCATGTCGTAGTCGCGCGCCATC  >  NC_002947/3462545‑3462834
                                                                                                                                               |                                                                                                                                                   
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:74898/149‑1 (MQ=255)
                                            gggATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTc                                                                                                   <  2:204560/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                                                                                                                                         >  8:304040/1‑59 (MQ=255)
                                                                                                gcagcaTGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGcc                                                                                                                                         <  7:304040/59‑1 (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)
                                                                                                                                           tCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCGGCGACAGGGTGACCGGGTAGCCGGTGACGATCATGTCGTAGTCgcgcgc      >  3:187423/1‑148 (MQ=255)
                                                                                                                                           tCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCGGCGACAGGGTGACCGGGTAGCCGGTGACGATCATGTCGTAGTCGCGCGcc     <  4:146962/149‑1 (MQ=255)
                                                                                                                                               gtggtCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCGGCGACAGGGTGACCGGGTAGCCGGTGACGATCATGTCGTAGTCGCGCGCCATc  >  3:191759/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
AGCGGTTCCACCAGACGGTGGAGCTGCCCGGTGGGTAATAGTTGGGGATCCAGTAGTAATTCCATTGCAGCACGCGGTCCAGGGCATGGGCGTGGCGCAGCATGTCGGCCTGGGTGTCGGCGCGGACCAGGCCGTCGATCAGGTGGTCTACGGCCGGGTCCTGCAAAACCATCAGGTTGTTCGAGCCGGGGTCGTGAGCGGCCGCCGAACCAAAGTAGTTGTACAACTCGGCGCCCGGCGACAGGGTGACCGGGTAGCCGGTGACGATCATGTCGTAGTCGCGCGCCATC  >  NC_002947/3462545‑3462834

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: