Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,494,843 G→A 5.9% A381A (GCC→GCT cobJ ← precorrin‑3B C17‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,494,8430GA 5.9% 1114.4 / 36.3 352A381A (GCC→GCTcobJprecorrin‑3B C17‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (162/169);  new base A (14/7);  total (176/176)
Fisher's exact test for biased strand distribution p-value = 1.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 357 total aligned reads displayed.

GCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGA  >  NC_002947/5494694‑5494991
                                                                                                                                                     |                                                                                                                                                    
  cAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGgc                                                                                                                                                    <  2:5903402/150‑1 (MQ=255)
  cAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGgc                                                                                                                                                    >  1:1002056/1‑150 (MQ=255)
  cAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGgc                                                                                                                                                    >  2:184437/1‑150 (MQ=255)
         gCGGAGACCCCCGGCAAGATCTGCCGGTCACCGCGATGCCCGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGt                                                                                                                                             <  1:6462881/150‑1 (MQ=255)
            gATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCg                                                                                                                                          >  2:728668/1‑150 (MQ=255)
            gATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGCCCTTGGGCCGCCAGTTCg                                                                                                                                          >  1:2701167/1‑150 (MQ=255)
              tACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGTGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGaa                                                                                                                                        <  1:1183712/150‑1 (MQ=255)
              tACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGaa                                                                                                                                        <  2:2328285/150‑1 (MQ=255)
               aCCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCCGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAg                                                                                                                                       <  2:7653521/150‑1 (MQ=255)
               aCCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAg                                                                                                                                       <  2:3183338/150‑1 (MQ=255)
                 ccccGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGc                                                                                                                                     <  2:5574648/150‑1 (MQ=255)
                 ccccGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGc                                                                                                                                     >  2:6447818/1‑150 (MQ=255)
                  cccGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCa                                                                                                                                    >  1:6116735/1‑150 (MQ=255)
                        aaGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCgg                                                                                                                              <  2:459803/150‑1 (MQ=255)
                        aaGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCgg                                                                                                                              >  2:7030098/1‑150 (MQ=255)
                          gATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGgc                                                                                                                            >  2:3934188/1‑150 (MQ=255)
                          gATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGgc                                                                                                                            <  2:7855354/150‑1 (MQ=255)
                          gATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGgc                                                                                                                            >  1:5603507/1‑150 (MQ=255)
                             cTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGATCGAAGGCATGTCGGgcgcc                                                                                                                         >  2:5050136/1‑149 (MQ=255)
                              tGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGgcgcgc                                                                                                                        >  2:1060065/1‑150 (MQ=255)
                              tGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGgcgcgc                                                                                                                        >  1:838220/1‑150 (MQ=255)
                                cAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTg                                                                                                                      <  2:3972630/150‑1 (MQ=255)
                                  ggTCAACGCGATGCCTGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCa                                                                                                                    <  1:6602351/150‑1 (MQ=255)
                                     cAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCAtct                                                                                                                 >  2:5005513/1‑150 (MQ=255)
                                       aCGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATctct                                                                                                               >  1:5352989/1‑150 (MQ=255)
                                        cgcgATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCtt                                                                                                              >  1:7850228/1‑150 (MQ=255)
                                            aTGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTcgcg                                                                                                          <  1:3612987/150‑1 (MQ=255)
                                               cctggTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGtt                                                                                                       <  2:42151/147‑1 (MQ=255)
                                                   gTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCg                                                                                                   <  1:4962359/150‑1 (MQ=255)
                                                   gTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCg                                                                                                   <  1:7390447/150‑1 (MQ=255)
                                                    tGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCgg                                                                                                  <  2:1462966/150‑1 (MQ=255)
                                                     ggggTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCTCCGCCACCGAACGGCCTTGGGCCGCCAGTTCGCAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgt                                                                                                 <  2:6800209/150‑1 (MQ=255)
                                                      gggTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtg                                                                                                <  2:7507277/150‑1 (MQ=255)
                                                         tCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgcag                                                                                             >  2:3045614/1‑150 (MQ=255)
                                                             tGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgcagtgca                                                                                         >  2:7371513/1‑150 (MQ=255)
                                                             tGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgcagtgca                                                                                         >  1:36530/1‑150 (MQ=255)
                                                              ggATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGtgcagtgcag                                                                                        >  2:8082301/1‑150 (MQ=255)
                                                              cgATACGTGCAATGCTTCGAGCACCCAAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCCCCACCGAACGGCCTTCGGCCGCCAGATCTAAGGCATGTCGGCCGCGCTCCATCACATCGCGGTTGACGGagcagtgcag                                                                                        >  2:440260/2‑150 (MQ=255)
                                                               gATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTgcagc                                                                                       <  1:8307257/150‑1 (MQ=255)
                                                                 ttCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCAc                                                                                     >  1:7666840/1‑150 (MQ=255)
                                                                   cGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCt                                                                                   <  1:3216564/150‑1 (MQ=255)
                                                                    gTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTg                                                                                  <  2:5741921/150‑1 (MQ=255)
                                                                        aaTGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTcg                                                                              >  1:4915572/1‑150 (MQ=255)
                                                                           gCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTcgccg                                                                           >  1:7240566/1‑150 (MQ=255)
                                                                            cTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTcgccgc                                                                          <  2:3332982/150‑1 (MQ=255)
                                                                            cTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTcgccgc                                                                          <  2:7796300/150‑1 (MQ=255)
                                                                             ttCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCcgcg                                                                         <  1:7449855/150‑1 (MQ=255)
                                                                             ttCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCcgcg                                                                         >  2:3118032/1‑150 (MQ=255)
                                                                                gAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGaa                                                                      >  1:9066089/1‑150 (MQ=255)
                                                                                gAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGaa                                                                      <  2:4500972/150‑1 (MQ=255)
                                                                                 aGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAAc                                                                     >  1:2286445/1‑150 (MQ=255)
                                                                                   cACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACgg                                                                   <  1:2872477/150‑1 (MQ=255)
                                                                                      cGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGgcc                                                                >  2:6288129/1‑150 (MQ=255)
                                                                                       gCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGTAGCACCTGGTCGCCGCGGAACGGgcca                                                               <  2:8835117/150‑1 (MQ=255)
                                                                                       gCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGgcca                                                               <  2:3451064/150‑1 (MQ=255)
                                                                                        cAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGgccag                                                              >  2:658025/1‑150 (MQ=255)
                                                                                            ggCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCAt                                                          <  2:1634746/150‑1 (MQ=255)
                                                                                              ccATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCTAGCCATgc                                                        <  1:5265474/150‑1 (MQ=255)
                                                                                              ccATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATgc                                                        <  1:6994957/150‑1 (MQ=255)
                                                                                               cATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATgcg                                                       >  2:2768455/1‑150 (MQ=255)
                                                                                                   gCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACg                                                   >  1:7104777/1‑150 (MQ=255)
                                                                                                     gAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTa                                                 >  2:5003288/1‑150 (MQ=255)
                                                                                                     gAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTa                                                 >  1:869531/1‑150 (MQ=255)
                                                                                                     aaataCGCCCGGATCGCCTGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATa                                                 >  2:2544747/5‑150 (MQ=255)
                                                                                                      aaaaCGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTaa                                                <  1:1619660/150‑1 (MQ=255)
                                                                                                        taCGCCCGGATCGCCTGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGt                                              <  1:4427159/149‑1 (MQ=255)
                                                                                                         aCGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGtt                                             <  1:6740138/150‑1 (MQ=255)
                                                                                                             ccGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGt                                         >  2:4128882/1‑150 (MQ=255)
                                                                                                             ccGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGt                                         >  2:6177682/1‑150 (MQ=255)
                                                                                                               gggTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTaa                                       <  1:3934188/150‑1 (MQ=255)
                                                                                                                   cGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGa                                   <  2:282845/150‑1 (MQ=255)
                                                                                                                    gCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAg                                  >  1:8994380/1‑150 (MQ=255)
                                                                                                                     ccTGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGTTTCATAGCCAAGc                                 >  2:2712344/1‑150 (MQ=255)
                                                                                                                     ccTGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGTTTCATAGCCAAGc                                 >  1:1971787/1‑150 (MQ=255)
                                                                                                                       tGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGcac                               >  1:632291/1‑150 (MQ=255)
                                                                                                                           gacACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCt                           >  2:4906447/1‑150 (MQ=255)
                                                                                                                             caccacCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTc                         <  1:1529563/150‑1 (MQ=255)
                                                                                                                               ccaccgccaccgAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCgg                       <  2:2193355/150‑1 (MQ=255)
                                                                                                                                caccgccactgAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGc                      <  1:6291892/150‑1 (MQ=255)
                                                                                                                                caccgccaccgAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGc                      >  2:7159801/1‑150 (MQ=255)
                                                                                                                                caccgccaccgAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGc                      <  1:3054988/150‑1 (MQ=255)
                                                                                                                                caccgccaccgAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGc                      >  1:2026112/1‑150 (MQ=255)
                                                                                                                                 accgccaccgAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGcc                     >  2:4950260/1‑150 (MQ=255)
                                                                                                                                 accgccaccgAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGcc                     >  2:2941807/1‑150 (MQ=255)
                                                                                                                                   cgccaccgAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTg                   >  1:7146367/1‑150 (MQ=255)
                                                                                                                                    gccaccgAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTgc                  >  2:8263779/1‑150 (MQ=255)
                                                                                                                                     ccaccgAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTgcg                 >  2:4792945/1‑150 (MQ=255)
                                                                                                                                        ccgAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCaa              <  1:8619116/150‑1 (MQ=255)
                                                                                                                                        catcaCGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGTTTc                                           <  2:8444519/117‑1 (MQ=255)
                                                                                                                                             cGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCg         <  1:6006167/150‑1 (MQ=255)
                                                                                                                                               gCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCg                                                                                                   <  1:2527109/58‑1 (MQ=255)
                                                                                                                                               gCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCg                                                                                                   >  2:2527109/1‑58 (MQ=255)
                                                                                                                                               gCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGc       <  1:4277383/150‑1 (MQ=255)
                                                                                                                                                ccTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGcc      <  1:3512499/150‑1 (MQ=255)
                                                                                                                                                ccTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGcc      <  2:838220/150‑1 (MQ=255)
                                                                                                                                                 cTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCt     <  2:7534799/150‑1 (MQ=255)
                                                                                                                                                  ttGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCTTCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCtt    <  1:7531317/150‑1 (MQ=255)
                                                                                                                                                  ttGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCtt    <  2:5352989/150‑1 (MQ=255)
                                                                                                                                                  ttGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCtt    <  2:4720740/150‑1 (MQ=255)
                                                                                                                                                    gggCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGa  >  2:5148610/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGA  >  NC_002947/5494694‑5494991

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: