Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,545,946 T→G 6.9% I349M (ATT→ATG)  PP_2238 → metallopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,545,9460TG 6.9% 275.0 / 8.8 87I349M (ATT→ATG) PP_2238metallopeptidase
Reads supporting (aligned to +/- strand):  ref base T (43/38);  new base G (3/3);  total (46/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

CCTGCGGGGGCGCGGCTGCTGGTGGACCCGACCCGGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAATTTCGTGTCAC  >  NC_002947/2545800‑2546091
                                                                                                                                                  |                                                                                                                                                 
ccTGCGGGGGCGCGGCTGCTGGTGGACCCGACCCGGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCATCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTc                                                                                                                                                  <  1:75288/148‑1 (MQ=255)
ccTGCGGGGGCGCGGCTGCTGGTGGACCCGACCCGGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCg                                                                                                                                                 >  8:7062/1‑149 (MQ=255)
       gggCGCGGCTGCTGGGGGACCCGACCCGGGTCACCTGCGGCCGGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGCGAATATGCACCAGGTc                                                                                                                                          <  8:242219/149‑1 (MQ=255)
           gcgGCTGCTGGTGGACCCGACCCGGGTCACCTGCGGCCTGCTGGACAACCAGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGTATATTCGCCAGGTCATg                                                                                                                                       >  1:258232/1‑148 (MQ=255)
             ggCTGCTGGTGGACCCGACCCGGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAg                                                                                                                                    >  8:39771/1‑149 (MQ=255)
              gctgctGGTGGACCCGACCCGGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGc                                                                                                                                   <  6:178563/149‑1 (MQ=255)
                  ctgGTGGACCCGACCCGGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAg                                                                                                                                 <  8:101516/147‑1 (MQ=255)
                   tggtggACCCGACCCGGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGAc                                                                                                                              >  6:327801/1‑149 (MQ=255)
                    ggtggACCCGACCCGGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACg                                                                                                                             >  6:70196/1‑149 (MQ=255)
                          cccgacccgGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAg                                                                                                                       >  8:157739/1‑149 (MQ=255)
                          cccgacccgGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACAACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAg                                                                                                                       >  8:30391/1‑149 (MQ=255)
                           ccgacccgGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGc                                                                                                                      >  6:26642/1‑149 (MQ=255)
                           ccgacccgGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGc                                                                                                                      >  3:51516/1‑149 (MQ=255)
                            cgacccgGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGcc                                                                                                                     >  1:78662/1‑149 (MQ=255)
                             gacccgGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCt                                                                                                                    >  8:40568/1‑149 (MQ=255)
                              acccgGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCtt                                                                                                                   <  7:30391/149‑1 (MQ=255)
                              acccgGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCtt                                                                                                                   >  6:314759/1‑149 (MQ=255)
                                  gggTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGc                                                                                                               <  5:314759/149‑1 (MQ=255)
                                    gTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGc                                                                                                               >  6:343602/1‑147 (MQ=255)
                                         cTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAACCCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGtttt                                                                                                        >  5:180640/1‑149 (MQ=255)
                                            cGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTc                                                                                                      <  5:343602/148‑1 (MQ=255)
                                                 tgctgGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGt                                                                                                           >  1:93650/1‑138 (MQ=255)
                                                 tgctgGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGt                                                                                                           <  2:93650/138‑1 (MQ=255)
                                                   ctgGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTgg                                                                                                <  7:157739/147‑1 (MQ=255)
                                                     ggACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTc                                                                                             >  7:240207/1‑148 (MQ=255)
                                                          aCCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGc                                                                                        >  6:182830/1‑148 (MQ=255)
                                                              ggcggcAGAGGTGGTTCTTGTCGAGGTGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGaa                                                                                     >  8:165016/1‑147 (MQ=255)
                                                              ggcggcAGAGGTGGTTCTTGTCGAGGTGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGaa                                                                                     <  7:165016/147‑1 (MQ=255)
                                                                cggcAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGCCTCGAAccag                                                                                 <  4:312369/149‑3 (MQ=255)
                                                                cggcAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGCCTCGAAccag                                                                                 <  5:26642/149‑3 (MQ=255)
                                                                      agGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCgg                                                                           <  5:70196/149‑1 (MQ=255)
                                                                       ggtggtTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCggg                                                                          <  5:327801/149‑1 (MQ=255)
                                                                       ggtggtTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCgg                                                                           <  5:182830/148‑1 (MQ=255)
                                                                         tggtTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGa                                                                         <  1:232469/148‑1 (MQ=255)
                                                                         tggtTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAg                                                                        <  7:40568/149‑1 (MQ=255)
                                                                          ggtTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAgg                                                                       >  7:237831/1‑149 (MQ=255)
                                                                           gtTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGTGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGt                                                                      >  6:186746/1‑149 (MQ=255)
                                                                              cTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGtcat                                                                   >  2:187145/1‑149 (MQ=255)
                                                                                tGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGtcatca                                                                 >  5:251747/1‑149 (MQ=255)
                                                                                 gTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGTGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGtcatc                                                                  <  5:186746/147‑1 (MQ=255)
                                                                                 gTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCAc                                                                >  2:77655/1‑149 (MQ=255)
                                                                                  tCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCAcc                                                               <  1:187145/149‑1 (MQ=255)
                                                                                  tCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCAcc                                                               >  1:260944/1‑149 (MQ=255)
                                                                                  tCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCATGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCAcc                                                               >  8:262902/1‑149 (MQ=255)
                                                                                     aGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAg                                                            >  2:180347/1‑149 (MQ=255)
                                                                                      ggggCTAAACACCACCACACTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGCAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCTAGGCGAACCTCGGCCGGGAGGTCATCACCGAGc                                                           >  3:193189/1‑149 (MQ=255)
                                                                                           tAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACg                                                      >  2:26476/1‑149 (MQ=255)
                                                                                              aCCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGtt                                                   >  1:303901/1‑149 (MQ=255)
                                                                                              aCCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGtt                                                   <  2:303901/149‑1 (MQ=255)
                                                                                                 ccaccaCGCTTAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGAc                                                <  4:51516/149‑1 (MQ=255)
                                                                                                 ccaccaCGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGAc                                                <  6:251747/149‑1 (MQ=255)
                                                                                                 ccaccaCGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGAc                                                <  1:77655/149‑1 (MQ=255)
                                                                                                 ccaccaCGCTGAGCAAGTCCTACAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCAACGAGCTGACGGTTGAc                                                >  8:220694/1‑149 (MQ=255)
                                                                                                   accacGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGa                                              <  7:262902/149‑1 (MQ=255)
                                                                                                    ccacGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAg                                             >  2:305434/1‑149 (MQ=255)
                                                                                                       cGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTCGAACCAGGCGaaa                                                                                    >  5:326990/1‑106 (MQ=255)
                                                                                                       cGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTCGAACCAGGCGaaa                                                                                    <  6:326990/107‑2 (MQ=255)
                                                                                                        gCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAAt                                         >  2:116059/1‑149 (MQ=255)
                                                                                                        gCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAAt                                         <  7:7062/149‑1 (MQ=255)
                                                                                                         cTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAAtt                                        <  7:39771/149‑1 (MQ=255)
                                                                                                           gAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGtc                                                                     >  5:82982/1‑118 (MQ=255)
                                                                                                           gAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGtc                                                                     <  6:82982/118‑1 (MQ=255)
                                                                                                            aGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGa                                      <  1:180347/148‑1 (MQ=255)
                                                                                                            aGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGa                                      <  8:240207/148‑1 (MQ=255)
                                                                                                            aGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGa                                      <  2:78662/148‑1 (MQ=255)
                                                                                                              cAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTg                                   <  1:305434/149‑1 (MQ=255)
                                                                                                                        aaGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGc                         <  1:26476/149‑1 (MQ=255)
                                                                                                                                cgacgacCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGc                 <  2:260944/149‑1 (MQ=255)
                                                                                                                                ccatgacCTGGCGAATATGCACCAGGTCTTGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGc                 >  1:192324/5‑149 (MQ=255)
                                                                                                                                 gacgacCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGcc                <  7:220694/149‑1 (MQ=255)
                                                                                                                                 gacgacCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGcc                >  6:58505/1‑149 (MQ=255)
                                                                                                                                 gacgacCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGcc                >  2:300469/1‑149 (MQ=255)
                                                                                                                                 gacgacCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGACCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGcc                >  3:336193/1‑149 (MQ=255)
                                                                                                                                       cTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAAttt          >  7:303373/1‑149 (MQ=255)
                                                                                                                                        tGGAGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGGTGACGGTTGACGAGCAATTGAGTGCCGCTCGTGCTCGACGGCCTAAttac         >  6:989/1‑147 (MQ=255)
                                                                                                                                         ggcgaatatGCACCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGc                         <  7:91336/127‑1 (MQ=255)
                                                                                                                                         ggcgaatatGCACCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGc                         <  7:86608/127‑1 (MQ=255)
                                                                                                                                         ggcgaatatGCACCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGc                         >  8:86608/6‑132 (MQ=255)
                                                                                                                                         ggcgaatatGCACCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGc                         >  8:91336/6‑132 (MQ=255)
                                                                                                                                         ggTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAATTTCg        >  5:44466/1‑149 (MQ=255)
                                                                                                                                          gTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAATTTc         <  8:303373/147‑1 (MQ=255)
                                                                                                                                          gTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAATTTCgt       <  6:44466/149‑1 (MQ=255)
                                                                                                                                             cATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGc                 <  2:24306/136‑1 (MQ=255)
                                                                                                                                             cATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGc                 >  1:24306/1‑136 (MQ=255)
                                                                                                                                               tatTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAATTTCGTGTc    >  5:55883/1‑147 (MQ=255)
                                                                                                                                               tatTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAATTTCGTGTc    <  6:55883/147‑1 (MQ=255)
                                                                                                                                               tatTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAATTTCGTGTCAc  >  6:228805/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CCTGCGGGGGCGCGGCTGCTGGTGGACCCGACCCGGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAATTTCGTGTCAC  >  NC_002947/2545800‑2546091

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