Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,902,672 T→C 9.6% V196A (GTG→GCG) ‡ PP_4310 → hydantoin racemase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,902,6720TC 9.6% 199.5 / 10.1 73V196A (GTG→GCG) ‡PP_4310hydantoin racemase
Reads supporting (aligned to +/- strand):  ref base T (33/33);  new base C (5/2);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

TGGAACTGGAAGCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTCGGCCGCTAGAGATCC  >  NC_002947/4902525‑4902820
                                                                                                                                                   |                                                                                                                                                    
tggaactggaaGCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTg                                                                                                                                                     >  7:276014/1‑149 (MQ=255)
 ggaactggaaGCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGc                                                                                                                                                    >  6:37575/1‑149 (MQ=255)
  gaactggaaGCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGcc                                                                                                                                                   <  1:7933/149‑1 (MQ=255)
     ctggaaGCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCggt                                                                                                                                                >  7:233165/1‑149 (MQ=255)
        gaaGCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGTGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCggtggt                                                                                                                                             <  8:112354/149‑1 (MQ=255)
        gaaGCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCggtggt                                                                                                                                             <  2:332409/149‑1 (MQ=255)
         aaGCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGgtggtg                                                                                                                                            >  3:338988/1‑149 (MQ=255)
            ccGACCCGCGGCGGGCGGTGGAGGCTATTGTCGAGCAGGCCGATGGCGCGGGGCGATACGAGAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGAATGGCCTGGCTGGACGAGCAGATTCGCCAACGCAGCGGCGTGCCGGTGGTGGAt                                                                                                                                         <  6:33095/149‑1 (MQ=255)
             cGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATg                                                                                                                                        <  4:61890/149‑1 (MQ=255)
             cGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATg                                                                                                                                        <  4:303028/149‑1 (MQ=255)
               aCCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCg                                                                                                                                                  >  5:234001/1‑137 (MQ=255)
               aCCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCg                                                                                                                                                  <  6:234001/137‑1 (MQ=255)
                   gCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGa                                                                                                                                  <  6:326484/149‑1 (MQ=255)
                   gCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGa                                                                                                                                  <  4:338988/149‑1 (MQ=255)
                   gCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGa                                                                                                                                  <  5:37575/149‑1 (MQ=255)
                     tGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGt                                                                                                                                <  2:289985/149‑1 (MQ=255)
                        gggCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCACGCCggtgcgttg                                                                                                                                          >  6:164325/1‑131 (MQ=255)
                        gggCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCACGCCggtgcgttg                                                                                                                                          <  5:164325/136‑6 (MQ=255)
                         ggCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGTAg                                                                                                                            <  1:68796/149‑1 (MQ=255)
                                gAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAGGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGAc                                                                                                                     <  3:280443/149‑1 (MQ=255)
                                     gATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTg                                                                                                                >  6:144395/1‑149 (MQ=255)
                                       ttGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGcc                                                                                                              >  6:276258/1‑149 (MQ=255)
                                          tCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGaa                                                                                                           <  7:9018/149‑1 (MQ=255)
                                            gAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATc                                                                                                         >  3:199687/1‑149 (MQ=255)
                                            gAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATc                                                                                                         >  5:113335/1‑149 (MQ=255)
                                            gAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCAAGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATc                                                                                                         >  1:108206/1‑149 (MQ=255)
                                                  gCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGc                                                                                                       <  3:24398/145‑1 (MQ=255)
                                                  gCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGc                                                                                                       >  4:24398/1‑145 (MQ=255)
                                                   ccGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTg                                                                                                  <  3:142763/149‑1 (MQ=255)
                                                         gcgcCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATg                                                                                            <  1:134040/149‑1 (MQ=255)
                                                            ccGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGt                                                                                         <  8:143427/149‑1 (MQ=255)
                                                              gTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGa                                                                                                                                  <  2:254813/106‑1 (MQ=255)
                                                              gTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGa                                                                                                                                  >  1:254813/1‑106 (MQ=255)
                                                               tGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTg                                                                                      <  5:144395/149‑1 (MQ=255)
                                                                            ggCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAg                                                                         >  6:208749/1‑149 (MQ=255)
                                                                             gCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGt                                                                        <  5:276258/149‑1 (MQ=255)
                                                                                gAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTgcg                                                                     >  6:174445/1‑149 (MQ=255)
                                                                                      aTCTGCCTGGGGTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCt                                                                <  8:233165/148‑1 (MQ=255)
                                                                                        cTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATg                                                             >  7:38540/1‑149 (MQ=255)
                                                                                           ccTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGc                                                            <  8:276014/147‑1 (MQ=255)
                                                                                               ggTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCa                                                                   >  7:133515/1‑136 (MQ=255)
                                                                                               ggTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCa                                                                   <  8:133515/136‑1 (MQ=255)
                                                                                                gTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCa                                                     <  8:38540/149‑1 (MQ=255)
                                                                                                    tgGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACggg                                                 <  6:113335/149‑1 (MQ=255)
                                                                                                    tgGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACggg                                                 <  4:199687/149‑1 (MQ=255)
                                                                                                      gcggcATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGc                                                >  6:191825/1‑148 (MQ=255)
                                                                                                       cggcATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCga                                              <  5:208749/149‑1 (MQ=255)
                                                                                                       cggcATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCg                                               >  6:175113/1‑148 (MQ=255)
                                                                                                                       gACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTg                              >  6:55131/1‑149 (MQ=255)
                                                                                                                       gACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTg                              >  7:255392/1‑149 (MQ=255)
                                                                                                                           aGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGcc                           >  5:106252/1‑148 (MQ=255)
                                                                                                                            gCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGa                         >  7:152672/1‑149 (MQ=255)
                                                                                                                                aTTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATgcg                     >  2:180413/1‑149 (MQ=255)
                                                                                                                                aTTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATgcg                     >  2:209234/1‑149 (MQ=255)
                                                                                                                                    gCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACgg                                                  <  5:109183/116‑1 (MQ=255)
                                                                                                                                    gCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACgg                                                  >  6:109183/1‑116 (MQ=255)
                                                                                                                                    gCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTc                 >  4:187341/1‑149 (MQ=255)
                                                                                                                                     ccAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGa                                                          >  6:4115/1‑107 (MQ=255)
                                                                                                                                     ccAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGa                                                          <  5:4115/107‑1 (MQ=255)
                                                                                                                                      tccaccACCGGCACGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTCgg               >  2:181300/6‑149 (MQ=255)
                                                                                                                                      cAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGc                            <  5:132670/136‑1 (MQ=255)
                                                                                                                                      cAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGc                            >  6:132670/1‑136 (MQ=255)
                                                                                                                                      cAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTCgg               >  3:61029/1‑149 (MQ=255)
                                                                                                                                      cAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGAAGCGCTTCgg               >  5:242312/1‑149 (MQ=255)
                                                                                                                                          ccACCGGCACGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAgtggtg                                   <  6:62544/124‑1 (MQ=255)
                                                                                                                                          ccACCGGCACGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAgtggtg                                   >  5:62544/2‑125 (MQ=255)
                                                                                                                                          ccACCGGCACGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTCGGCCGc           >  3:46108/2‑149 (MQ=255)
                                                                                                                                          ccACCGGCACGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTCGGCCGc           >  8:289937/2‑149 (MQ=255)
                                                                                                                                            aCCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGc                                                <  6:31629/110‑1 (MQ=255)
                                                                                                                                            aCCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGc                                                >  5:31629/1‑110 (MQ=255)
                                                                                                                                            aCCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTCGGCCGc           >  6:302357/1‑147 (MQ=255)
                                                                                                                                            aCCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCAATGCGCTTCGGCCGc           >  4:323858/1‑147 (MQ=255)
                                                                                                                                                gCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGAAGCGCTTCGGCCGCTAgaga     <  6:242312/149‑1 (MQ=255)
                                                                                                                                                   tGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCAATGCGCTTCGGCCGCTAGAGATcc  <  3:323858/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGGAACTGGAAGCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTCGGCCGCTAGAGATCC  >  NC_002947/4902525‑4902820

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: