Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,264,296 C→A 7.3% intergenic (‑50/+34) zinU ← / ← PP_4640 P‑loop guanosine triphosphatase‑dependent zinc chaperone/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,264,2960CA 7.3% 211.5 / 5.7 68intergenic (‑50/+34)zinU/PP_4640P‑loop guanosine triphosphatase‑dependent zinc chaperone/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/30);  new base A (2/3);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

GAGCTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGTTCGACA  >  NC_002947/5264151‑5264444
                                                                                                                                                   |                                                                                                                                                     
gAGCTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCa                                                                                                                                                    <  8:119365/149‑1 (MQ=255)
    ttCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCt                                                                                                                                                >  7:97171/1‑149 (MQ=255)
         gCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCTGGGTAACCGGCTCCcac                                                                                                                                           <  6:103573/149‑1 (MQ=255)
           cATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAg                                                                                                                                         >  8:115213/1‑149 (MQ=255)
           cATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAg                                                                                                                                         >  8:120123/1‑149 (MQ=255)
           cATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAg                                                                                                                                         >  8:247418/1‑149 (MQ=255)
                   gCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTacac                                                                                                                                >  8:87971/1‑149 (MQ=255)
                         agaaTGTACTTGAGCAAGGTGGTCTTGGCGGCAGCGAGGAAGCCTGTGAGGACGGTTTCGGGAATTGGCATCTGCAAGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTAGCCAGGc                                                                                                                                                 <  6:263509/124‑1 (MQ=255)
                         aGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGACTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCAGAAAATGACTCTGGGAGACGGTTACCAAgac                                                                                                                                                 >  5:263509/1‑125 (MQ=255)
                           cATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCa                                                                                                                        <  3:329420/149‑1 (MQ=255)
                           cATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCa                                                                                                                        <  8:97171/149‑1 (MQ=255)
                            aTGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCTGGGTAACCGGCTCCcaca                                                                                                                                          <  1:275668/131‑1 (MQ=255)
                            aTGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCTGGGTAACCGGCTCCcaca                                                                                                                                          >  2:275668/1‑131 (MQ=255)
                            aTGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAg                                                                                                                       <  7:115213/149‑1 (MQ=255)
                            aTGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAg                                                                                                                       <  7:120123/149‑1 (MQ=255)
                                     aGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAg                                                                                                              <  7:87971/149‑1 (MQ=255)
                                          ggtggtCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATCGACTGTGGGAGCCGGTTACCCAGGCTcc                                                                                                                                              <  1:197270/114‑1 (MQ=255)
                                          ggtggtCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTcc                                                                                                                                              >  2:197270/1‑113 (MQ=255)
                                           gtggtCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTg                                                                                                        >  7:136131/1‑149 (MQ=255)
                                            tggtCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGTCATCTGCACGCAAAATACTCCTTAATCTGAAAATGTCTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTgg                                                                                                       <  1:277874/149‑1 (MQ=255)
                                            tggtCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTgg                                                                                                       >  8:149171/1‑149 (MQ=255)
                                              gtCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTggg                                                                                                      <  7:284145/148‑1 (MQ=255)
                                                      ggCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACtt                                                                                             <  1:333452/149‑1 (MQ=255)
                                                         aCCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAgg                                                                                                                      >  3:63261/1‑121 (MQ=255)
                                                         aCCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAgg                                                                                                                      <  4:63261/121‑1 (MQ=255)
                                                                    cTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTCCCACCGTAACGt                                                                                <  7:149171/148‑1 (MQ=255)
                                                                     tgtgAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGc                                                                              >  2:188852/1‑149 (MQ=255)
                                                                                    gggAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAAc                                                                                  <  6:49523/130‑1 (MQ=255)
                                                                                    gggAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAAc                                                                                  >  5:49523/1‑130 (MQ=255)
                                                                                      gAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCgg                                                              >  5:158903/1‑148 (MQ=255)
                                                                                           ggCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTgg                                                                                                       <  5:154779/102‑1 (MQ=255)
                                                                                           ggCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTgg                                                                                                       >  6:154779/1‑102 (MQ=255)
                                                                                                 tGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCg                                                  >  7:49325/1‑149 (MQ=255)
                                                                                                   cACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCg                                                  <  7:247418/147‑1 (MQ=255)
                                                                                                        aaaaTACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCtc                                            >  2:324693/1‑148 (MQ=255)
                                                                                                        aaaaTACTACAGAATCTGAAAATGACTGTGGGAGACCGTTACCCAGGCTCCCACAGGTTTACACGAACAGCAGCCTCAACAGCACTTGGGCCCGTACTTGCCAACGTAACGTCCTTCCTGACGTTCTAGTAAGAATGACTCTTAGCtca                                           >  8:177872/1‑149 (MQ=255)
                                                                                                              cTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCg                                     <  1:188852/149‑1 (MQ=255)
                                                                                                               tCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCtc                                            <  7:41969/141‑1 (MQ=255)
                                                                                                               tCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCtc                                            >  8:41969/1‑141 (MQ=255)
                                                                                                                 cTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGc                                   >  1:158285/1‑148 (MQ=255)
                                                                                                                 cTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGc                                   <  2:158285/148‑1 (MQ=255)
                                                                                                                  tGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGgaa                                                            <  6:50400/122‑1 (MQ=255)
                                                                                                                  tGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCtt                                 >  5:104959/1‑149 (MQ=255)
                                                                                                                  tGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCtt                                 >  5:159791/1‑149 (MQ=255)
                                                                                                                  tGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGATTCCTGACGTTCGCGgaa                                                            >  5:50400/1‑122 (MQ=255)
                                                                                                                     aTCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTc                              <  8:49325/149‑1 (MQ=255)
                                                                                                                      tCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACGACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTc                              >  1:257251/1‑149 (MQ=255)
                                                                                                                      tCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTcc                             >  1:124005/1‑149 (MQ=255)
                                                                                                                        tGAGAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCg                            >  5:173618/1‑148 (MQ=255)
                                                                                                                          agaaTGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCg                            <  6:173618/144‑1 (MQ=255)
                                                                                                                                  cTGTGGGAGCCTGGGTAACCGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGAt                         >  1:45147/1‑142 (MQ=255)
                                                                                                                                  cTGTGGGAGCCTGGGTAACCGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGAt                         <  2:45147/142‑1 (MQ=255)
                                                                                                                                   tgtgGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTc                 >  5:174313/1‑149 (MQ=255)
                                                                                                                                     tgGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGtttt                  <  3:69944/146‑1 (MQ=255)
                                                                                                                                     tgGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCt                <  6:104959/148‑1 (MQ=255)
                                                                                                                                     tgGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGATCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGtttt                  >  4:69944/1‑146 (MQ=255)
                                                                                                                                       ggAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCCCCGTAACGAGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCa             >  3:99039/1‑149 (MQ=255)
                                                                                                                                       ggAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCg                                                  >  6:45156/1‑112 (MQ=255)
                                                                                                                                       ggAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCg                                                  <  5:45156/112‑1 (MQ=255)
                                                                                                                                       ggAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGc              <  6:158903/148‑1 (MQ=255)
                                                                                                                                        gAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCAt            <  2:257251/149‑1 (MQ=255)
                                                                                                                                        gAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCAt            <  6:159791/149‑1 (MQ=255)
                                                                                                                                        gAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCAt            <  2:124005/149‑1 (MQ=255)
                                                                                                                                              ggTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCAt            <  8:4891/144‑1 (MQ=255)
                                                                                                                                              ggTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGtt       >  5:104693/1‑149 (MQ=255)
                                                                                                                                              ggTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGGTTTCTGCAt            >  7:4891/1‑144 (MQ=255)
                                                                                                                                               gTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGTTc      >  3:331489/1‑149 (MQ=255)
                                                                                                                                                ttACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGTTCg     <  4:331489/149‑1 (MQ=255)
                                                                                                                                                   cccATGATCCCACATGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCAACGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTTTCCGTATGCTTTTTCTGCATGTGTTCGACa  <  7:177872/149‑1 (MQ=255)
                                                                                                                                                   cccAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGTTCGAc   >  3:21325/1‑148 (MQ=255)
                                                                                                                                                   cccAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGTTCGAc   <  4:21325/148‑1 (MQ=255)
                                                                                                                                                   cccAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGTTCGACa  <  1:324693/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GAGCTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGTTCGACA  >  NC_002947/5264151‑5264444

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