Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 447,159:1 +C 5.4% intergenic (+39/‑263) 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,1591.C 5.4% 381.5 / 15.4 92intergenic (+39/‑263)PP_0368/PP_0369acyl‑CoA dehydrogenase/diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base . (50/37);  new base C (3/2);  total (53/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

ACTTCCAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAACCCGCTCCTAC  >  NC_002947/447014‑447303
                                                                                                                                                  |                                                                                                                                                
aCTTCCAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAgg                                                                                                                                               >  3:257742/1‑149 (MQ=255)
  ttCCAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAg                                                                                                                                                >  5:50420/1‑146 (MQ=255)
  ttCCAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAg                                                                                                                                                <  6:50420/146‑1 (MQ=255)
    ccAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCgg                                                                                                                                            <  3:34478/148‑1 (MQ=255)
       gcgcATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGttt                                                                                                                                        >  5:98290/1‑149 (MQ=255)
          cATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGtttttt                                                                                                                                     <  7:39320/149‑1 (MQ=255)
           aTCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTa                                                                                                                                    <  1:110027/149‑1 (MQ=255)
           aTCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTa                                                                                                                                    >  8:153665/1‑149 (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                                                                                                                                      >  4:7092/1‑141 (MQ=38)
                gCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTCCCGGTCAGGGCAGGCGGGttttg                                                                                                                                      <  3:7092/142‑2 (MQ=38)
                gCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTCCCGGTCAGGGCAGGCGGGttttg                                                                                                                                      <  7:289487/142‑2 (MQ=38)
                gCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTCCCGGTCAGGGCAGGCGGGttttg                                                                                                                                      >  8:289487/1‑141 (MQ=38)
                   aCGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCa                                                                                                                            <  5:379700/149‑1 (MQ=255)
                     gtacccgtacccATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAgg                                                                                                                          >  3:60107/1‑149 (MQ=255)
                         ccgtacccATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGtt                                                                                                                        >  5:306652/1‑147 (MQ=255)
                         ccgtacccATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGtt                                                                                                                        >  5:311761/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                                                                                                         <  5:307824/149‑1 (MQ=255)
                                      gCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCg                                                                                                         >  1:243272/1‑149 (MQ=255)
                                               cTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGATCgg                                                                                                <  5:260195/149‑1 (MQ=255)
                                               cTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCgg                                                                                                <  8:198967/149‑1 (MQ=255)
                                               cTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCgg                                                                                                <  7:153665/149‑1 (MQ=255)
                                                 gTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGt                                                                                              >  8:48334/1‑149 (MQ=255)
                                                           aacaacCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGa                                                                                    <  2:243272/149‑1 (MQ=255)
                                                                 cTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATTCTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTc                                                                                                                             <  8:226318/102‑1 (MQ=255)
                                                                 cTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATTCTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTc                                                                                                                             >  7:226318/1‑102 (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                                                                         <  7:83346/136‑1 (MQ=255)
                                                                                   gagCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAggccggcc                                                                         >  8:83346/1‑136 (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                                                              >  5:18663/1‑126 (MQ=255)
                                                                                                        tGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCaa                                                              <  6:18663/126‑1 (MQ=255)
                                                                                                        tGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCAc                                            >  8:90252/1‑144 (MQ=255)
                                                                                                        tGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCAc                                            <  7:90252/144‑1 (MQ=255)
                                                                                                          aTACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGc                                      <  6:306652/148‑1 (MQ=255)
                                                                                                          aTACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGc                                      <  6:311761/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:354174/1‑149 (MQ=255)
                                                                                                                         ccAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTc                      >  6:70425/1‑149 (MQ=255)
                                                                                                                         ccAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTc                      >  8:122189/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                  >  2:141741/1‑147 (MQ=255)
                                                                                                                               cccGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCgg                  <  1:141741/147‑1 (MQ=255)
                                                                                                                                ccGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGt                >  7:257403/1‑148 (MQ=255)
                                                                                                                                ccGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGt                >  6:207538/1‑148 (MQ=255)
                                                                                                                                ccGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGt                >  6:158693/1‑148 (MQ=255)
                                                                                                                                ccGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGt                >  2:343478/1‑148 (MQ=255)
                                                                                                                                   ccTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAAc            <  4:173897/149‑1 (MQ=255)
                                                                                                                                    cTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAAcc           <  4:164417/149‑1 (MQ=255)
                                                                                                                                       cccTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAACCCGc        <  8:201049/149‑1 (MQ=255)
                                                                                                                                          tGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCAcc                                           <  6:116725/111‑1 (MQ=255)
                                                                                                                                          tGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAACCCGCTc      >  3:199590/1‑148 (MQ=255)
                                                                                                                                          agACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCAcc                                           >  5:116725/2‑111 (MQ=255)
                                                                                                                                               gggAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAACCCGCTCCTAc  >  2:124897/1‑147 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
ACTTCCAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAACCCGCTCCTAC  >  NC_002947/447014‑447303

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: