Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 447,148 Δ1 bp 5.6% intergenic (+28/‑274) PP_0368 → / → PP_0369 acyl‑CoA dehydrogenase/diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947447,1480G. 5.6% 291.7 / 9.5 89intergenic (+28/‑274)PP_0368/PP_0369acyl‑CoA dehydrogenase/diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (48/36);  new base . (3/2);  total (51/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

CGCTTCTACTTCCAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAACC  >  NC_002947/447007‑447294
                                                                                                                                             |                                                                                                                                                   
cgcTTCTACTTCCAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGACGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGAc                                                                                                                                              <  6:42452/149‑1 (MQ=255)
       aCTTCCAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAgg                                                                                                                                      >  3:257742/1‑149 (MQ=255)
         ttCCAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAg                                                                                                                                       <  6:50420/146‑1 (MQ=255)
         ttCCAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAg                                                                                                                                       >  5:50420/1‑146 (MQ=255)
           ccAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCgg                                                                                                                                   <  3:34478/148‑1 (MQ=255)
              gcgcATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGttt                                                                                                                               >  5:98290/1‑149 (MQ=255)
                 cATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGtttttt                                                                                                                            <  7:39320/149‑1 (MQ=255)
                  aTCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTa                                                                                                                           >  8:153665/1‑149 (MQ=255)
                  aTCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTa                                                                                                                           <  1:110027/149‑1 (MQ=255)
                   tCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAGCCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTAc                                                                                                                          <  7:78584/149‑1 (MQ=255)
                    ccTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACg                                                                                                                         >  3:213963/1‑149 (MQ=255)
                      tGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGtttttt                                                                                                                            <  1:180620/144‑1 (MQ=255)
                      tGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGtttttt                                                                                                                            >  2:180620/1‑144 (MQ=255)
                       gCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTCCCGGTCAGGGCAGGCGGGttttg                                                                                                                             <  7:289487/142‑2 (MQ=38)
                       gCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTCCCGGTCAGGGCAGGCGGGttttg                                                                                                                             >  8:289487/1‑141 (MQ=38)
                       gCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTCCCGGTCAGGGCAGGCGGGttttg                                                                                                                             >  4:7092/1‑141 (MQ=38)
                       gCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTCCCGGTCAGGGCAGGCGGGttttg                                                                                                                             <  3:7092/142‑2 (MQ=38)
                          aCGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCa                                                                                                                   <  5:379700/149‑1 (MQ=255)
                            gtacccgtacccATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAgg                                                                                                                 >  3:60107/1‑149 (MQ=255)
                                ccgtacccATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGtt                                                                                                               >  5:311761/1‑147 (MQ=255)
                                ccgtacccATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGtt                                                                                                               >  5:306652/1‑147 (MQ=255)
                                  gtacccATGTGGTGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGc                                                                                                           >  6:337654/1‑149 (MQ=255)
                                  gtacccATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAg                                                                                                                                       <  7:226533/121‑1 (MQ=255)
                                  gtacccATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAg                                                                                                                                       >  8:226533/1‑121 (MQ=255)
                                     cccATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTg                                                                                                        >  8:97961/1‑149 (MQ=255)
                                      ccATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCa                                                                                                                   >  2:68649/1‑137 (MQ=255)
                                      ccATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCa                                                                                                                   <  1:68649/137‑1 (MQ=255)
                                             gCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCg                                                                                                >  1:243272/1‑149 (MQ=255)
                                             gCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCg                                                                                                <  5:307824/149‑1 (MQ=255)
                                                      cTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGATCgg                                                                                       <  5:260195/149‑1 (MQ=255)
                                                      cTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCgg                                                                                       <  7:153665/149‑1 (MQ=255)
                                                      cTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCgg                                                                                       <  8:198967/149‑1 (MQ=255)
                                                        gTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGt                                                                                     >  8:48334/1‑149 (MQ=255)
                                                                  aacaacCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGa                                                                           <  2:243272/149‑1 (MQ=255)
                                                                        cTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATTCTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTc                                                                                                                    >  7:226318/1‑102 (MQ=255)
                                                                        cTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATTCTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTc                                                                                                                    <  8:226318/102‑1 (MQ=255)
                                                                        cTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAATACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAggc                                                                     <  6:332659/149‑1 (MQ=255)
                                                                          gATGGCGATGGACGAAGAGCACTTCGGCATGTCGTACTGATACTGGCGCTGTACCCAAAACCCGACTGCCATGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCAGCGAAGAggccg                                                                   >  4:138852/1‑149 (MQ=255)
                                                                           aTGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAggccg                                                                   <  7:48334/148‑1 (MQ=255)
                                                                                gATGGACGAAGAGCACTTCGGCCTGTAGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCctg                                                             >  4:124753/1‑149 (MQ=255)
                                                                                acaggccGAAGTGCTCTTCGTCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCctg                                                             >  2:159004/7‑149 (MQ=39)
                                                                                   ggACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAggccagcagggccg                                                          >  8:64571/1‑142 (MQ=255)
                                                                                   ggACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCctgctg                                                          >  7:201049/1‑149 (MQ=255)
                                                                                    gACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCGGGCCctgctg                                                          <  5:337614/148‑1 (MQ=255)
                                                                                    gACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTgg                                                         <  4:60107/149‑1 (MQ=255)
                                                                                     aCGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGc                                                        <  8:153904/149‑1 (MQ=255)
                                                                                        aaGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCaa                                                     <  7:97961/149‑1 (MQ=255)
                                                                                         agagCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAAt                                                    >  2:261939/1‑149 (MQ=255)
                                                                                          gagCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAggccggcc                                                                >  8:83346/1‑136 (MQ=255)
                                                                                          gagCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAggccggcc                                                                <  7:83346/136‑1 (MQ=255)
                                                                                          gagCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAAtc                                                   <  4:213963/149‑1 (MQ=255)
                                                                                             cACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTcc                                                <  4:257742/149‑1 (MQ=255)
                                                                                                  cGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCCGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAg                                           >  8:63750/1‑149 (MQ=255)
                                                                                                    gCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGcc                                         >  3:164417/1‑149 (MQ=255)
                                                                                                            tACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCg                                 <  6:104130/149‑1 (MQ=255)
                                                                                                             aCTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTCCTGCTGGCCAATCTCCATCAGCCGTTCACCGt                                >  8:25011/1‑149 (MQ=255)
                                                                                                             aCTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCg                                 <  3:124753/148‑1 (MQ=255)
                                                                                                              cTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGtt                               >  6:112301/1‑149 (MQ=255)
                                                                                                               tGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCaa                                                     <  6:18663/126‑1 (MQ=255)
                                                                                                               tGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCaa                                                     >  5:18663/1‑126 (MQ=255)
                                                                                                               tGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCAc                                   >  8:90252/1‑144 (MQ=255)
                                                                                                               tGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCAc                                   <  7:90252/144‑1 (MQ=255)
                                                                                                                 aTACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGc                             <  6:311761/148‑1 (MQ=255)
                                                                                                                 aTACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGc                             <  6:306652/148‑1 (MQ=255)
                                                                                                                 aTACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGc                             >  7:219600/1‑148 (MQ=255)
                                                                                                                  tACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGc                             >  2:341815/1‑147 (MQ=255)
                                                                                                                   aCTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGGGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAAtc                                                   >  3:328033/1‑124 (MQ=255)
                                                                                                                   aCTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAAtc                                                   <  4:328033/124‑1 (MQ=255)
                                                                                                                         gcTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGGTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTAc                    >  6:44265/1‑149 (MQ=255)
                                                                                                                           tGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCg                  >  6:145482/1‑149 (MQ=255)
                                                                                                                           tGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCg                  >  7:17206/1‑149 (MQ=255)
                                                                                                                               cccAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTAtt              >  3:173897/1‑149 (MQ=255)
                                                                                                                               cccAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCACTCGTAGGAGCGGGTTTAGCCTCGAAGAGGCCGGCCGTGCTGGCCAATCTCCATCAGCAGTTCAACGTTGCCCGCGCTACCGTAtt              >  4:28111/1‑149 (MQ=255)
                                                                                                                                ccAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTc             >  7:332026/1‑149 (MQ=255)
                                                                                                                                ccAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTc             >  8:122189/1‑149 (MQ=255)
                                                                                                                                ccAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTc             >  6:70425/1‑149 (MQ=255)
                                                                                                                                ccAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTc             >  8:354174/1‑149 (MQ=255)
                                                                                                                                 cAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTTCGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGGCAATCTCCATGAGGGGTTCACCGTTGCCAGCGCTACCGTATTcg            >  3:345049/1‑149 (MQ=38)
                                                                                                                                 cAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTcg            >  3:345048/1‑149 (MQ=255)
                                                                                                                                 aaaaaCCCGCCTCCCGGTCAGGGCAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCGGCACTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTc             >  2:193864/2‑148 (MQ=25)
                                                                                                                                     aCCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCggg        <  8:219600/149‑1 (MQ=255)
                                                                                                                                      cccGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCgg         <  1:141741/147‑1 (MQ=255)
                                                                                                                                      cccGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCgg         >  2:141741/1‑147 (MQ=255)
                                                                                                                                       ccGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGt       >  2:343478/1‑148 (MQ=255)
                                                                                                                                       ccGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGt       >  7:257403/1‑148 (MQ=255)
                                                                                                                                       ccGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGt       >  6:158693/1‑148 (MQ=255)
                                                                                                                                       ccGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGt       >  6:207538/1‑148 (MQ=255)
                                                                                                                                          ccTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAAc   <  4:173897/149‑1 (MQ=255)
                                                                                                                                           cTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAAcc  <  4:164417/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
CGCTTCTACTTCCAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAACC  >  NC_002947/447007‑447294

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