Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 274,785 T→G 5.9% D388A (GAC→GCC)  PP_0222 ← DszA family monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947274,7850TG 5.9% 452.4 / 13.4 135D388A (GAC→GCC) PP_0222DszA family monooxygenase
Reads supporting (aligned to +/- strand):  ref base T (77/50);  new base G (4/4);  total (81/54)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 136 total aligned reads displayed.

GCGGGTGGTCGGCCTGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGCGCGTGGCTGA  >  NC_002947/274636‑274933
                                                                                                                                                     |                                                                                                                                                    
gCGGGTGGTCGGCCTGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGt                                                                                                                                                      <  1:141402/150‑1 (MQ=255)
 cGGGTGGTCGGCCTGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTc                                                                                                                                                     <  2:790921/150‑1 (MQ=255)
 cGGGTGGTCGGCCTGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTc                                                                                                                                                     >  1:3955921/1‑150 (MQ=255)
         cGGCCTGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCgg                                                                                                                                             >  1:1701575/1‑150 (MQ=255)
             cTGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTc                                                                                                                                         >  1:2757823/1‑150 (MQ=255)
             cTGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTc                                                                                                                                         >  1:2757534/1‑150 (MQ=255)
              tGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCa                                                                                                                                        >  1:4221398/1‑150 (MQ=255)
              tGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCa                                                                                                                                        >  1:4246272/1‑150 (MQ=255)
              tGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCa                                                                                                                                        <  2:4814971/150‑1 (MQ=255)
               ggAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCAt                                                                                                                                       >  2:1696166/1‑150 (MQ=255)
               ggAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCAt                                                                                                                                       >  2:1696189/1‑150 (MQ=255)
               ggAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGGCCTGGACGGTTTCaa                                                                                                                                       <  2:2055969/150‑2 (MQ=255)
               ggAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGGCCTGGACGGTTTCaa                                                                                                                                       <  2:4699607/150‑2 (MQ=255)
                gAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATc                                                                                                                                      <  1:982199/150‑1 (MQ=255)
                gAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATc                                                                                                                                      >  1:326136/1‑150 (MQ=255)
                gAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGGCCTGGACGGTTTCaac                                                                                                                                      <  2:3817725/150‑3 (MQ=255)
                 aaCAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGATGAAACCGGCCTGGACGGTTTCaacc                                                                                                                                     <  1:3438924/150‑4 (MQ=255)
                       ttttCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCAc                                                                                                                               <  2:1785953/150‑1 (MQ=255)
                        tttCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGTGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCAcg                                                                                                                              >  2:672128/1‑150 (MQ=255)
                         ttCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCAcgc                                                                                                                             <  1:370691/150‑1 (MQ=255)
                         ttCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCAcgc                                                                                                                             <  2:252710/150‑1 (MQ=255)
                         ttCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCAcgc                                                                                                                             <  1:3322928/150‑1 (MQ=255)
                           cgcgCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGa                                                                                                                           <  1:1634712/150‑1 (MQ=255)
                             cgcAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATc                                                                                                                         >  1:865467/1‑150 (MQ=255)
                               cAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCaa                                                                                                                       <  2:2259523/150‑1 (MQ=255)
                               cAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCaa                                                                                                                       <  1:3910173/150‑1 (MQ=255)
                                      acTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCACCTGCTCg                                                                                                                >  2:2210058/2‑150 (MQ=255)
                                        ttCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGc                                                                                                              >  1:4632538/1‑150 (MQ=255)
                                        ttCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGc                                                                                                              >  1:4000537/1‑150 (MQ=255)
                                        ttCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGc                                                                                                              >  1:4266262/1‑150 (MQ=255)
                                        ttCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGc                                                                                                              >  2:4477636/1‑150 (MQ=255)
                                        ttCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGc                                                                                                              <  2:1559658/150‑1 (MQ=255)
                                        ttCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGc                                                                                                              >  1:2341595/1‑150 (MQ=255)
                                        ttCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGc                                                                                                              >  1:2363513/1‑150 (MQ=255)
                                           ggTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCAc                                                                                                           <  2:1295264/150‑1 (MQ=255)
                                           accgaaGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCAc                                                                                                           >  1:1775673/6‑150 (MQ=255)
                                             tgtAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCt                                                                                                         <  2:1326739/150‑1 (MQ=255)
                                                 ggCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTc                                                                                                     >  2:2425392/1‑150 (MQ=255)
                                                   cAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAg                                                                                                   <  2:376548/150‑1 (MQ=255)
                                                     gTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGgc                                                                                                 <  2:3184518/150‑1 (MQ=255)
                                                      tCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGgcg                                                                                                >  2:1628811/1‑150 (MQ=255)
                                                        ttGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGgcgcg                                                                                              >  1:5003642/1‑150 (MQ=255)
                                                        ttGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGgcgcg                                                                                              >  2:4431024/1‑150 (MQ=255)
                                                         tGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAACTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGgcgcgc                                                                                             <  2:1699055/150‑1 (MQ=255)
                                                          gTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGcc                                                                                            >  2:3540190/1‑150 (MQ=255)
                                                          gTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGcc                                                                                            >  1:4916523/1‑150 (MQ=255)
                                                          gTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGcc                                                                                            <  2:1612845/150‑1 (MQ=255)
                                                          gTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGcc                                                                                            >  1:2072643/1‑150 (MQ=255)
                                                               ggCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCa                                                                                       >  1:3416150/1‑150 (MQ=255)
                                                                        gctgAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGt                                                                              <  2:655674/150‑1 (MQ=255)
                                                                         ctgAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTa                                                                             <  1:1377299/150‑1 (MQ=255)
                                                                         ctgAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTa                                                                             >  1:302287/1‑150 (MQ=255)
                                                                          tgAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAg                                                                            >  1:304552/1‑150 (MQ=255)
                                                                                       tGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAgg                                                               >  1:2418662/1‑150 (MQ=255)
                                                                                          ccAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGt                                                            >  1:41315/1‑150 (MQ=255)
                                                                                          ccAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGt                                                            >  1:41503/1‑150 (MQ=255)
                                                                                           cAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTg                                                           <  1:2872644/150‑1 (MQ=255)
                                                                                              cTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTg                                                        >  2:158688/1‑150 (MQ=255)
                                                                                               tCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGt                                                       >  2:3665201/1‑150 (MQ=255)
                                                                                                cGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGtt                                                      <  2:932012/150‑1 (MQ=255)
                                                                                                cGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGtt                                                      <  1:4477636/150‑1 (MQ=255)
                                                                                                  aTGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGc                                                    >  1:1776229/1‑150 (MQ=255)
                                                                                                      aGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACa                                                >  1:2490220/1‑150 (MQ=255)
                                                                                                         ccTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGgc                                             <  2:1363798/150‑1 (MQ=255)
                                                                                                         ccTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGgc                                             <  1:578369/150‑1 (MQ=255)
                                                                                                         ccTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGgc                                             <  2:1364118/150‑1 (MQ=255)
                                                                                                          cTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGgcg                                            >  2:3557331/1‑150 (MQ=255)
                                                                                                           tCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGgcgc                                           >  2:4484217/1‑150 (MQ=255)
                                                                                                           tCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGgcgc                                           >  1:2106061/1‑150 (MQ=255)
                                                                                                             gAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACa                                                >  1:3679648/1‑143 (MQ=255)
                                                                                                             gAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACa                                                <  2:3679648/143‑1 (MQ=255)
                                                                                                             gAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCcg                                         >  1:3141096/1‑150 (MQ=255)
                                                                                                             gAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCcg                                         >  2:2946655/1‑150 (MQ=255)
                                                                                                             gAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCcg                                         >  2:3598563/1‑150 (MQ=255)
                                                                                                                  tctcCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCggg                                    >  2:3912873/1‑150 (MQ=255)
                                                                                                                    tcCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTc                                  <  1:2373073/150‑1 (MQ=255)
                                                                                                                    tcCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTc                                  >  2:3557569/1‑150 (MQ=255)
                                                                                                                    tcCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTc                                  >  1:3953016/1‑150 (MQ=255)
                                                                                                                            gACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGtt                          <  1:3572272/150‑1 (MQ=255)
                                                                                                                            gACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGtt                          <  2:2473533/150‑1 (MQ=255)
                                                                                                                             aCGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTg                         >  1:1601845/1‑150 (MQ=255)
                                                                                                                               ggTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTc                       >  2:2987056/1‑150 (MQ=255)
                                                                                                                                 tGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCtt                     <  2:90198/150‑1 (MQ=255)
                                                                                                                                 tGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCtt                     <  2:2757823/150‑1 (MQ=255)
                                                                                                                                 tGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCtt                     <  2:2757534/150‑1 (MQ=255)
                                                                                                                                 tGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCtt                     >  2:2592067/1‑150 (MQ=255)
                                                                                                                                  gCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTg                    <  2:338383/150‑1 (MQ=255)
                                                                                                                                  gCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTg                    <  2:4221398/150‑1 (MQ=255)
                                                                                                                                       tCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACt               >  1:3029441/1‑150 (MQ=255)
                                                                                                                                       gatcgatGAAACCGGCCTGGACGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACt               >  2:2779477/7‑150 (MQ=255)
                                                                                                                                         aGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTgg             >  1:111365/1‑150 (MQ=255)
                                                                                                                                           ggtGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGgcg           >  1:2118486/3‑150 (MQ=255)
                                                                                                                                            ttGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGgcgc          >  2:420329/1‑150 (MQ=255)
                                                                                                                                            atGAAACCGGCCTGGACGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGgcgc          >  2:1084784/2‑150 (MQ=255)
                                                                                                                                            atGAAACCGGCCTGGACGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGgcgc          >  1:3244420/2‑150 (MQ=255)
                                                                                                                                               aaaCCGGCCTGGACGGTTTCAACGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGCGCGTg       >  2:3246097/1‑150 (MQ=255)
                                                                                                                                                aaCCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGCGCGTgg      >  1:3372113/1‑150 (MQ=255)
                                                                                                                                                  ccGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGCGCGTGGCt    <  1:2438431/150‑1 (MQ=255)
                                                                                                                                                   cGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGCGCGTGGCTg   >  1:416641/1‑150 (MQ=255)
                                                                                                                                                    gTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGCGCGTGGCTGa  <  2:4916523/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCGGGTGGTCGGCCTGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGCGCGTGGCTGA  >  NC_002947/274636‑274933

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