Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,062,900 C→G 6.3% G480A (GGC→GCC)  PP_3581 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,062,9000CG 6.3% 201.0 / 3.4 63G480A (GGC→GCC) PP_3581membrane protein
Reads supporting (aligned to +/- strand):  ref base C (30/29);  new base G (2/2);  total (32/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.25e-01

TCACCGCCATGCGCCAGCCATGGTGCCCCTGGTCCAGCCACTGACGCATCTGACGGCAGGCCTCAGTGAGTACCCACTCACCAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCG  >  NC_002947/4062757‑4063047
                                                                                                                                               |                                                                                                                                                   
tCACCGCCATGCGCCAGCCATGGTGCCCCTGGTCCAGCCACTGACGCATCTGACGGCAGGCCTCAGTGAGTACCCACTCACCAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCctgc                                                                                                                                                >  1:10598/1‑149 (MQ=255)
     gCCATGCGCCAGCCATGGTGCCCCTGGTCCAGCCACTGACGCATCTGACGGCAGGCCTCAGTGAGTACCCACTCACCAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCgg                                                                                                                                           <  7:211919/149‑1 (MQ=255)
     gCCATGCGCCAGCCATGGTGCCCCTGGTCCAGCCACTGACGCATCTGACGGCAGGCCTCAGTGAGTACCCACTCACCAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCgg                                                                                                                                           <  4:195997/149‑1 (MQ=255)
       cATGCGCCAGCCATGGTGCCCCTGGTCCAGCCACTGACGCATCTGACGGCAGGCCTCAGTGAGTACCCACTCACCAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGt                                                                                                                                         >  1:113007/1‑149 (MQ=255)
            gccagccaTGGTGCCCCTGGTCCAGCCACTGACGCATCTGACGGCAGGCCTCAGTGAGTACCCACTCACCAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTcc                                                                                                                                    >  4:70278/1‑149 (MQ=255)
              cagccaTGGTGCCCCTGGTCCAGCCACTGACGCATCTGACGGCAGGCCTCAGTGAGTACCCACTCACTAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCa                                                                                                                                  >  8:23107/1‑149 (MQ=255)
                       tGCCCCTGGTCCAGCCACTGACGCATCTGACGGCAGGCCTCAGTGAGTACCCACTCACCAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCaa                                                                                                                         >  7:23866/1‑149 (MQ=255)
                                tCCAGCCACTGACGCATCTGACGGCAGGCCTCAGTGAGTACCCACTCACCAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCgg                                                                                                                 >  3:251725/1‑148 (MQ=255)
                                       aCTGACGCATCTGACGGCAGGCCTCAGTGAGTACCCACTCACCAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATc                                                                                                         <  8:23866/149‑1 (MQ=255)
                                           aCGCATCTGACGGCAGGCCTCAGTGAGTACCCACTCACCAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGtt                                                                                                     <  3:108379/149‑1 (MQ=255)
                                               aTCTGACGGCAGGCCTCAGTGAGTACCCACTCACCAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTc                                                                                                                   <  3:180922/131‑1 (MQ=255)
                                               aTCTGACGGCAGGCCTCAGTGAGTACCCACTCACCAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTc                                                                                                                   >  4:180922/1‑131 (MQ=255)
                                                tCTGACGGCAGGCCTCAGTGAGTACCCACTCACCAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTcc                                                                                                                                    >  5:69425/1‑113 (MQ=255)
                                                tCTGACGGCAGGCCTCAGTGAGTACCCACTCACCAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTcc                                                                                                                                    <  6:69425/113‑1 (MQ=255)
                                                  tGACGGCAGGCCTCAGTGAGTACCCACTCACCAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCTCAGCAACGCTTCGGCACCGATCGGTTGGCAGGc                                                                                              >  7:240235/1‑149 (MQ=255)
                                                       gCAGGCCTCAGTGAGTACCCACTCACCAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGc                                                                                                                                       <  3:226993/103‑1 (MQ=255)
                                                       gCAGGCCTCAGTGAGTACCCACTCACCAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGc                                                                                                                                       >  4:226993/1‑103 (MQ=255)
                                                           gCCTCAGTGAGTACCCACTCACCAATGGGGATGACCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTc                                                                                     >  2:88430/1‑149 (MQ=255)
                                                            ccTCAGTGAGTACCCACTCACCAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTc                                                                                     >  2:15142/1‑148 (MQ=255)
                                                              tCAGTGAGTACCCACTCACCAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGaa                                                                                  >  5:167056/1‑149 (MQ=255)
                                                                  ggAGTACCCACTCACCAATGGGGATGAGCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTg                                                                              <  7:51389/148‑1 (MQ=255)
                                                                     gTACCCACTCACCAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGc                                                                           <  4:251725/149‑1 (MQ=255)
                                                                       aCCCACTCACCAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTg                                                                         >  7:77610/1‑149 (MQ=255)
                                                                        cccACTCACCAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTgg                                                                        >  3:191487/1‑149 (MQ=255)
                                                                         ccACTCACCAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGt                                                                       <  5:136635/149‑1 (MQ=255)
                                                                            cTCACCAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAAt                                                                    <  2:113007/149‑1 (MQ=255)
                                                                            cTCACCAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAAt                                                                    >  1:99381/1‑149 (MQ=255)
                                                                               aCCAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTgg                                                                                          <  7:35208/124‑1 (MQ=255)
                                                                               aCCAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTgg                                                                                          >  8:35208/1‑124 (MQ=255)
                                                                                  aaTGGGGATGATCAGCCCGGTCTTTTCCGCTAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGc                                                                                              >  3:74579/1‑117 (MQ=255)
                                                                                  aaTGGGGATGATCAGCCCGGTCTTTTCCGCTAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGc                                                                                              <  4:74579/117‑1 (MQ=255)
                                                                                      gggATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGaa                                                          >  5:25183/1‑149 (MQ=255)
                                                                                       ggATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAAc                                                         >  1:226416/1‑149 (MQ=255)
                                                                                        gatgatCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACt                                                        >  1:238426/1‑149 (MQ=255)
                                                                                        gatgatAAGACCGGTCTTTTCCTCCAGGCCGATTAAGCGGTCATGCAGCAGCAGGCACTGCTCCGCCTGCTCCCAGCGCAGCAACGCAACCGCACAGAACACTTTGCAGGACTCGGCCTCCAACTAGGGCAGGTAATGCAACCGGAACt                                                        >  8:85098/1‑149 (MQ=255)
                                                                                          tgatCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAAc                                                         >  4:58102/1‑146 (MQ=255)
                                                                                          tgatCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAAc                                                         <  3:58102/146‑1 (MQ=255)
                                                                                          tgatCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTgg                                                      <  4:191487/149‑1 (MQ=255)
                                                                                           gatCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGgc                                                     >  8:18497/1‑149 (MQ=255)
                                                                                                gCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCgctgc                                                <  7:85098/149‑1 (MQ=255)
                                                                                                      tCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAgc                                          <  3:70278/149‑1 (MQ=255)
                                                                                                      tCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAgc                                          <  1:88430/149‑1 (MQ=255)
                                                                                                         tttCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCgcagc                                                                                                                           >  2:240692/1‑65 (MQ=255)
                                                                                                         tttCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCgcagc                                                                                                                           <  1:240692/65‑1 (MQ=255)
                                                                                                          ttCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTGCAgcgc                                        >  4:236350/1‑147 (MQ=255)
                                                                                                           tCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCa                                      <  7:18497/148‑1 (MQ=255)
                                                                                                               ccAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGGAGCAGCGCCAgttcc                                 <  2:238426/149‑5 (MQ=255)
                                                                                                                 aGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCaa                               <  6:167056/149‑1 (MQ=255)
                                                                                                                        tGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTgc                        <  5:119934/149‑1 (MQ=255)
                                                                                                                          aaGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTgca                       <  8:77610/148‑1 (MQ=255)
                                                                                                                           aGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTgcagc                     <  6:25183/149‑1 (MQ=255)
                                                                                                                                 cAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGc               <  1:60524/149‑1 (MQ=255)
                                                                                                                                  aGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTCGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGc               >  4:59504/1‑148 (MQ=255)
                                                                                                                                    gcagcagcaGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAAc            <  3:236350/149‑1 (MQ=255)
                                                                                                                                     ccgcagcagGGCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAgc                                          >  4:20988/3‑118 (MQ=255)
                                                                                                                                     ccgcagcagGGCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAgc                                          <  3:20988/116‑1 (MQ=255)
                                                                                                                                      agcagcagGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAActg          >  8:156302/1‑149 (MQ=255)
                                                                                                                                           cagGGCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTg                                                                         >  1:85348/1‑81 (MQ=255)
                                                                                                                                           cagGGCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTg                                                                         <  2:85348/81‑1 (MQ=255)
                                                                                                                                           cagGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGgc     >  3:99525/1‑149 (MQ=255)
                                                                                                                                           cagGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGgc     <  3:59504/149‑1 (MQ=255)
                                                                                                                                              gCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGgcgcg  >  8:171092/1‑149 (MQ=255)
                                                                                                                                              gCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGgcgcg  <  4:99525/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TCACCGCCATGCGCCAGCCATGGTGCCCCTGGTCCAGCCACTGACGCATCTGACGGCAGGCCTCAGTGAGTACCCACTCACCAATGGGGATGATCAGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCG  >  NC_002947/4062757‑4063047

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑