Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,821,953 T→A 11.8% L173Q (CTG→CAG)  metW → methionine biosynthesis protein MetW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,821,9530TA11.8% 259.3 / 19.7 91L173Q (CTG→CAG) metWmethionine biosynthesis protein MetW
Reads supporting (aligned to +/- strand):  ref base T (39/41);  new base A (3/8);  total (42/49)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCCGCTGGTACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATCGCACCATGCGT  >  NC_002947/5821804‑5822102
                                                                                                                                                     |                                                                                                                                                      
cTGCCGCTGGTACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCt                                                                                                                                                        <  2:1293346/150‑1 (MQ=255)
   ccGCTGGTACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGc                                                                                                                                                     <  2:574573/150‑1 (MQ=255)
      cTGGTACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGt                                                                                                                                                  <  1:705711/150‑1 (MQ=255)
         gtgccgggCGACCAAAGGCCGCATGCCGGTATCGGACTTCATCCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGCGCACGAACGCCGGGCCAATGTACTCGACCGCCAGGCGGTCGa                                                                                                                                               <  1:2130283/144‑1 (MQ=255)
         gTACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGCACATCCACTTCTGCACCTTAGCCGAATTCGACGCCCTTGTCCACGCAAGCCGGGCCAAGGTACTCGACCGCCAGGCGGTCGa                                                                                                                                               <  2:1091664/150‑1 (MQ=255)
         gTACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGa                                                                                                                                               >  2:1940943/1‑150 (MQ=255)
         gTACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCAGGCGGTCGa                                                                                                                                               <  1:1671577/150‑1 (MQ=255)
         gTACCCGGCCGCCAAAGTCCGCATGCCGGTATCGGACTCCATGCCGTATACCTGGTACAACACGTCGAACATCCACCTCTGCCCATTCGCCGACTTCTAACTCCCTGTGCTCGAACCCCGGGTCAAGGTACTCGCCCGCTAGTCGGTCGa                                                                                                                                               <  2:30629/150‑1 (MQ=255)
            ccTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCGGAGCCAAGTTACTCGACCGCCTGGCGGTCGACCa                                                                                                                                            >  2:737460/1‑150 (MQ=255)
            ccTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCAGGCGGTCGagta                                                                                                                                            <  2:1361694/150‑4 (MQ=255)
              tGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACt                                                                                                                                          >  1:1096782/1‑150 (MQ=255)
               ggCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACtt                                                                                                                                         <  1:1553712/150‑1 (MQ=255)
                     cAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCg                                                                                                                                   >  1:852207/1‑150 (MQ=255)
                      aaaGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGtt                                                                                                                                                 >  2:1801237/1‑134 (MQ=255)
                      aaaGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGtt                                                                                                                                                 <  1:1801237/135‑2 (MQ=255)
                      aaaGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGc                                                                                                                                  >  1:1656508/1‑150 (MQ=255)
                               aTGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCAGGCGGTCGag                                                                                                                                              >  2:182663/1‑128 (MQ=255)
                               aTGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCAGGCGGTCGag                                                                                                                                              <  1:182663/129‑2 (MQ=255)
                                cgacGGTATCTGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACgggtggg                                                                                                                        >  2:1726613/4‑150 (MQ=255)
                                    ggTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGCGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGCTGAAGAGCTGATCCGCGACGggtggggtgg                                                                                                                   >  1:2332799/1‑150 (MQ=255)
                                     gTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGc                                                                                                                  <  1:1940943/150‑1 (MQ=255)
                                        tCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCgg                                                                                                               >  1:2277298/1‑150 (MQ=255)
                                            aCTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTAt                                                                                                           <  2:1092695/150‑1 (MQ=255)
                                             cTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACGGCCTGGCGGTCGTCCGCGTGCACCGCAACGGGTGGGGTGGCCGGCTATg                                                                                                          <  1:1624667/150‑1 (MQ=255)
                                              ttCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATgg                                                                                                         >  1:1980768/1‑150 (MQ=255)
                                              ttCATGCCGTATACCTGGTACAACACGCCGAACAGCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGACCGCCGGGCCACGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATgg                                                                                                         <  2:362810/150‑1 (MQ=255)
                                                  tGCCGTATACCTGGTACAACACGACGTCCGCCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTa                                                                                                     >  1:1269045/1‑150 (MQ=255)
                                                          aCCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTa                                                                                             <  2:967001/150‑1 (MQ=255)
                                                           ccTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACgggtggg                                                                                                                        <  1:2320553/123‑1 (MQ=255)
                                                           ccTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACgggtggg                                                                                                                        >  2:2320553/1‑123 (MQ=255)
                                                           ccTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAg                                                                                            >  2:2401383/1‑150 (MQ=255)
                                                           ccTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAg                                                                                            >  2:458582/1‑150 (MQ=255)
                                                             tGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGt                                                                                          <  2:451336/150‑1 (MQ=255)
                                                             tGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGt                                                                                          <  1:1726613/150‑1 (MQ=255)
                                                              ggTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTg                                                                                         <  2:952103/150‑1 (MQ=255)
                                                               gTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTga                                                                                        <  1:458582/150‑1 (MQ=255)
                                                                   aacaCGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACgggt                                                                                                                           <  2:1044195/112‑1 (MQ=255)
                                                                   aacaCGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACgggt                                                                                                                           >  1:1044195/1‑112 (MQ=255)
                                                                   aacaCGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGTGTGGcggtca                                                                                                             >  1:806284/1‑121 (MQ=255)
                                                                   aacaCGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGTGTGGcggtca                                                                                                             <  2:806284/126‑6 (MQ=255)
                                                                   aacaCGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGt                                                                                          <  2:2214316/144‑1 (MQ=255)
                                                                   aacaCGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGt                                                                                          >  1:2214316/1‑144 (MQ=255)
                                                                   aacaCGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATc                                                                                    >  2:1811047/1‑150 (MQ=255)
                                                                   aacaCGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATc                                                                                    >  1:1454887/1‑150 (MQ=255)
                                                                      acGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGc                                                                                 >  2:971018/1‑150 (MQ=255)
                                                                         ccGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAg                                                                                            <  1:792428/136‑1 (MQ=255)
                                                                         ccGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAg                                                                                            >  2:792428/1‑136 (MQ=255)
                                                                            aaCATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCCAc                                                                           <  1:2401383/150‑1 (MQ=255)
                                                                             aCATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCAGGCGGTCGa                                                                                                                                               <  1:1253975/82‑1 (MQ=255)
                                                                             aCATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCAGGCGGTCGa                                                                                                                                               >  2:1253975/1‑82 (MQ=255)
                                                                                tCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACgggtggggtg                                                                                                                    <  2:533038/105‑1 (MQ=255)
                                                                                tCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACgggtggggtg                                                                                                                    >  1:533038/1‑105 (MQ=255)
                                                                                  cACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTc                                                                     <  2:1096782/150‑1 (MQ=255)
                                                                                   aCTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGACCGCAGATTCTGATGATGCTACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTca                                                                    <  2:1269045/150‑1 (MQ=255)
                                                                                   aCTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTca                                                                    <  2:1454887/150‑1 (MQ=255)
                                                                                      tCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTcagct                                                                 <  2:1656508/150‑2 (MQ=255)
                                                                                      tCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTcagca                                                                 <  1:971018/150‑1 (MQ=255)
                                                                                            ccTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCgg                                                           <  2:2277298/150‑1 (MQ=255)
                                                                                             cTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCt                                                                                                      >  1:1107297/1‑106 (MQ=255)
                                                                                             cTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCt                                                                                                      <  2:1107297/106‑1 (MQ=255)
                                                                                              ttCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGc                                                         >  1:647834/1‑150 (MQ=255)
                                                                                                cGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCt                                                       <  1:700109/150‑1 (MQ=255)
                                                                                                  ccGACTGCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGc                                                     <  2:233185/150‑1 (MQ=255)
                                                                                                   cGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCa                                                    >  2:1697038/1‑150 (MQ=255)
                                                                                                       ttCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGc                                                                                 >  2:841629/1‑117 (MQ=255)
                                                                                                       ttCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGc                                                                                 <  1:841629/117‑1 (MQ=255)
                                                                                                        tCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGc                                               <  2:2183281/150‑1 (MQ=255)
                                                                                                          gaagaaCTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGc                                               >  2:589088/1‑148 (MQ=255)
                                                                                                          gaagaaCTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGc                                               <  1:589088/148‑1 (MQ=255)
                                                                                                          gaagaaCTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCAt                                             >  1:475555/1‑150 (MQ=255)
                                                                                                           aagaaCTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCgg                                                           <  1:1033546/135‑1 (MQ=255)
                                                                                                           aagaaCTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCgg                                                           >  2:1033546/1‑135 (MQ=255)
                                                                                                           aagaaCTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATc                                            <  2:852207/150‑1 (MQ=255)
                                                                                                              aaCTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGc                                         <  2:2406152/150‑1 (MQ=255)
                                                                                                                cTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCTTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTc                                       <  1:737460/150‑1 (MQ=255)
                                                                                                                 ttGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGctgct                                                 <  2:2353794/139‑5 (MQ=255)
                                                                                                                 ttGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGctgct                                                 >  1:2353794/1‑135 (MQ=255)
                                                                                                                    tgCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCgg                                   <  2:1980768/150‑1 (MQ=255)
                                                                                                                     gCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTcgcgtt                                  <  2:643519/150‑3 (MQ=255)
                                                                                                                       aCGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTCTAGTTCAAAAACCTAGGTGAGGTCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCt                                >  1:1621977/1‑150 (MQ=255)
                                                                                                                       aCGAACGCCGGGCCAAGGTACTCGACCGACTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTGCTAGGTGAGATCGGCATCTACCGCGGCAGCCGTCTGGGCCTACCAGAGCATCAGCTCGCGGTCt                                >  1:947619/1‑150 (MQ=255)
                                                                                                                             gCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGcc                          >  2:277369/1‑150 (MQ=255)
                                                                                                                             gCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGcc                          >  1:2233592/1‑150 (MQ=255)
                                                                                                                             gCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGGGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCCTCAGCAGTCCGGGCCTGGAAGAGCACCAGCTCGCGGTCTGACGcc                          >  1:976890/1‑150 (MQ=255)
                                                                                                                                 ggCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCAcc                      >  1:121285/1‑150 (MQ=255)
                                                                                                                                     aaGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAaga                                                 >  2:281025/1‑119 (MQ=255)
                                                                                                                                     aaGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAaga                                                 <  1:281025/119‑1 (MQ=255)
                                                                                                                                     aaGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCggag                  >  1:1183330/1‑150 (MQ=255)
                                                                                                                                              gACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATCGCAc         <  1:277369/150‑1 (MQ=255)
                                                                                                                                               aCCGCCAGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTa                                                                                                     <  1:2377857/57‑1 (MQ=255)
                                                                                                                                               aCCGCCAGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTa                                                                                                     >  2:2377857/1‑57 (MQ=255)
                                                                                                                                                   ccTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATCGCACCATGc    >  1:259359/1‑150 (MQ=255)
                                                                                                                                                     tGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATCGCACCATGCGt  >  1:1424209/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CTGCCGCTGGTACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATCGCACCATGCGT  >  NC_002947/5821804‑5822102

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