Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,091,750 T→G 5.7% L40R (CTG→CGG)  PP_5343 → putative Transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,091,7500TG 5.7% 372.9 / 11.7 124L40R (CTG→CGG) PP_5343putative Transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (60/57);  new base G (4/3);  total (64/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.95e-01

Only 100 of 126 total aligned reads displayed.

GGATGGAACCATTCGACAAGGATCCCGCCATGTACAACAGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGCCTATGTCAGCCC  >  NC_002947/6091603‑6091891
                                                                                                                                                   |                                                                                                                                             
ggATGGAACCATTCGACAAGGATCCCGCCATGTACAACAGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGc                                                                                                                                             <  2:845693/150‑1 (MQ=255)
    ggAACCATTCGACAAGGATCCCGCCATGTACAACAGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGt                                                                                                                                         >  1:1371572/1‑150 (MQ=255)
      aaCCATTCGACAAGGATCCCGCCATGTACAACAGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGc                                                                                                                                             <  2:224829/144‑1 (MQ=255)
      aaCCATTCGACAAGGATCCCGCCATGTACAACAGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGc                                                                                                                                             >  1:224829/1‑144 (MQ=255)
      aaCCATTCGACAAGGATCCCGCCATGTACAACAGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGc                                                                                                                                       >  1:2174663/1‑150 (MQ=255)
          aTTCGACAAGGATCCCGCCATGTACAACAGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAAGGCGGACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGt                                                                                                                                   >  2:1770532/1‑150 (MQ=255)
               aCAAGGATCCCGCCATGTACAACAGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACCGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATa                                                                                                                              <  1:1770532/150‑1 (MQ=255)
                       cccGCCATGTACAACAGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAg                                                                                                                      <  1:1753524/150‑1 (MQ=255)
                        ccGCCATGTACAACAGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAgg                                                                                                                     >  2:1407170/1‑150 (MQ=255)
                        ccGCCATGTACAACAGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAgg                                                                                                                     <  2:1080998/150‑1 (MQ=255)
                                   aacaGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGc                                                                                                          >  2:155967/1‑150 (MQ=255)
                                    acaGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCa                                                                                                         <  2:1226026/150‑1 (MQ=255)
                                     caGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGcc                                                                                                                                            <  1:239074/114‑1 (MQ=255)
                                     caGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGcc                                                                                                                                            >  2:239074/1‑114 (MQ=255)
                                     caGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCAc                                                                                                        <  2:319284/150‑1 (MQ=255)
                                        tAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCATGGCGAACAGGGCCTTCTGGCCACGCACCTGCCgg                                                                                                                                          >  1:1738970/1‑113 (MQ=255)
                                        tAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCgg                                                                                                                                          <  2:1738970/113‑1 (MQ=255)
                                        tAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACtgt                                                                                                     >  1:221008/1‑150 (MQ=255)
                                                aTCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCa                                                                                                         >  2:1048210/1‑138 (MQ=255)
                                                aTCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCa                                                                                                         <  1:1048210/138‑1 (MQ=255)
                                                aTCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGccc                                                                                             <  1:2193789/150‑1 (MQ=255)
                                                 tCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTgtgt                                                                                                   <  1:1747350/143‑1 (MQ=255)
                                                 tCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTgtgt                                                                                                   >  2:1747350/1‑143 (MQ=255)
                                                   aGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCCCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCCCCTGCCGGTGCTGGTCGATACCGCCGAGGGCGCTTTCGCCACTGTGTATGCCCAcc                                                                                          >  2:1950840/1‑150 (MQ=255)
                                                    ggAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCt                                                                                         >  1:855490/1‑150 (MQ=255)
                                                       acacGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGa                                                                                                                       <  2:865695/117‑1 (MQ=255)
                                                       acacGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGACGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGc                                                                                      >  1:658589/1‑150 (MQ=255)
                                                        cacGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGcc                                                                                     <  1:155967/150‑1 (MQ=255)
                                                         acGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGccc                                                                                    >  1:2058827/1‑150 (MQ=255)
                                                           gACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCgg                                                                                  >  2:6896/1‑150 (MQ=255)
                                                           gACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGACGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCgg                                                                                  >  2:1767640/1‑150 (MQ=255)
                                                            aCCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCggg                                                                                 <  2:270271/150‑1 (MQ=255)
                                                               tCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGAGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGc                                                                                                                   >  1:1589375/1‑113 (MQ=255)
                                                               tCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGAGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGc                                                                                                                   <  2:1589375/113‑1 (MQ=255)
                                                                 aCCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCa                                                                                            <  2:855419/134‑1 (MQ=255)
                                                                 aCCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCa                                                                                            >  1:855419/1‑134 (MQ=255)
                                                                    cGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGc                                                                         <  2:1924077/150‑1 (MQ=255)
                                                                       ttGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCaa                                                                      <  1:1714695/150‑1 (MQ=255)
                                                                        tGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAAt                                                                     <  2:1760176/150‑1 (MQ=255)
                                                                         gCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATg                                                                    <  2:2058827/150‑1 (MQ=255)
                                                                             cAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACtgt                                                                                                     <  2:1768630/113‑1 (MQ=255)
                                                                             cAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACtgt                                                                                                     >  1:1768630/1‑113 (MQ=255)
                                                                             cAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGc                                                                         <  1:1521634/141‑1 (MQ=255)
                                                                             cAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCACGTCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCaa                                                                >  2:2269198/1‑150 (MQ=255)
                                                                                cacaTGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGAc                                                             <  2:221008/150‑1 (MQ=255)
                                                                                cacaTGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGAc                                                             <  1:1407170/150‑1 (MQ=255)
                                                                                     gCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTgca                                                        <  1:838436/150‑1 (MQ=255)
                                                                                      cTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTgcag                                                       <  2:800940/150‑1 (MQ=255)
                                                                                       tAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGAc                                                             <  1:1650892/143‑1 (MQ=255)
                                                                                       tAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTgcagc                                                      <  1:217897/150‑1 (MQ=255)
                                                                                         gACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCg                                                                                                            >  2:478610/1‑94 (MQ=255)
                                                                                         gACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCAc                                                    >  2:780741/1‑150 (MQ=255)
                                                                                         gACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCAc                                                    >  1:1404865/1‑150 (MQ=255)
                                                                                            acCCGCCTGGCCGTTCTGGTCAGCTAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATg                                                                    <  1:1372963/131‑1 (MQ=255)
                                                                                            acCCGCCTGGCCGTTCTGGTCAGCTAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATg                                                                    >  2:1372963/1‑131 (MQ=255)
                                                                                            acCCGCCTGGCCGTTCTGGTCAGCTAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATg                                                                    >  2:1373313/1‑131 (MQ=255)
                                                                                            acCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCAcggc                                                 >  2:2035860/1‑150 (MQ=255)
                                                                                            acCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCAcggc                                                 <  2:2174663/150‑1 (MQ=255)
                                                                                              ccGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGa                                                                                                                       >  1:1888012/1‑78 (MQ=255)
                                                                                              ccGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGa                                                                                                                       <  2:1888012/78‑1 (MQ=255)
                                                                                              ccGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAAcc                                                                           >  2:2174852/1‑122 (MQ=255)
                                                                                              ccGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCAcggcgg                                               >  2:1827363/1‑150 (MQ=255)
                                                                                               cGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACggcggc                                              >  2:1565288/1‑150 (MQ=255)
                                                                                               cGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACggcggc                                              <  1:1565288/150‑1 (MQ=255)
                                                                                                           ttGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCActgctg                                  <  1:1082856/149‑1 (MQ=255)
                                                                                                           cTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATTGCAAGACCTGCAGCACGGCGGCGAAGCActgctg                                  <  1:1769868/150‑1 (MQ=255)
                                                                                                            tGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTgg                                 >  2:1345834/1‑150 (MQ=255)
                                                                                                             ggTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGt                                <  1:1037052/150‑1 (MQ=255)
                                                                                                              gTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCa                                        <  1:1941429/141‑1 (MQ=255)
                                                                                                                  gCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTcc                           <  2:170266/150‑1 (MQ=255)
                                                                                                                   ccAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATg                                                                    >  1:2272594/1‑108 (MQ=255)
                                                                                                                   ccAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATg                                                                    <  2:2272594/108‑1 (MQ=255)
                                                                                                                    cAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCTCCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTcccc                         >  2:1830348/1‑150 (MQ=255)
                                                                                                                    cAGGACGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACTGCGGCGAAGCACTGCTGGTATTcccc                         <  1:1767640/150‑1 (MQ=255)
                                                                                                                       ggCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGgc                      >  2:499517/1‑150 (MQ=255)
                                                                                                                         cGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGgcgc                    <  1:499517/150‑1 (MQ=255)
                                                                                                                               gggCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGc              <  1:833949/150‑1 (MQ=255)
                                                                                                                                ggCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGcc             >  1:115182/1‑150 (MQ=255)
                                                                                                                                ggCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGcc             >  2:429341/1‑150 (MQ=255)
                                                                                                                                  ccTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACtgt                                                                                                     <  1:519418/60‑1 (MQ=255)
                                                                                                                                  ccTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACtgt                                                                                                     >  2:519418/1‑60 (MQ=255)
                                                                                                                                   cTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGCCTAt          <  1:1368248/150‑1 (MQ=255)
                                                                                                                                   cTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGCCTAt          <  2:115182/150‑1 (MQ=255)
                                                                                                                                      cTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGCCTATGTc       >  1:457707/1‑150 (MQ=255)
                                                                                                                                       tGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAggg                                                                                                                    >  1:1394556/1‑40 (MQ=255)
                                                                                                                                       tGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAggg                                                                                                                    <  2:1394556/40‑1 (MQ=255)
                                                                                                                                         gCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCActgc                                    >  1:2046555/1‑118 (MQ=255)
                                                                                                                                         gCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCAc                                       <  1:768932/115‑1 (MQ=255)
                                                                                                                                         gCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCAc                                       >  2:768932/1‑115 (MQ=255)
                                                                                                                                          accagCACCGGCAGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCActgctg                                  <  2:2390817/115‑1 (MQ=255)
                                                                                                                                          accagCACCGGCAGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCActgctg                                  >  1:2390817/5‑119 (MQ=255)
                                                                                                                                           ccagCACCGGCAGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGCCTAt          >  1:2021004/4‑142 (MQ=255)
                                                                                                                                           ccagCACCGGCAGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGCCTAt          <  2:2021004/139‑1 (MQ=255)
                                                                                                                                           ccagCACCGGCAGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGCCTATGTCAGccc  >  2:1656948/4‑150 (MQ=255)
                                                                                                                                           cACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGCCTATGTCAGccc  >  2:1254750/1‑150 (MQ=255)
                                                                                                                                           cACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGCCTATGTCAGccc  >  2:1584407/1‑150 (MQ=255)
                                                                                                                                              gCACCGGCAGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTa                               >  2:2081632/1‑118 (MQ=255)
                                                                                                                                              gCACCGGCAGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTa                               <  1:2081632/118‑1 (MQ=255)
                                                                                                                                                   tGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTgcagc                                                      <  2:1288207/90‑1 (MQ=255)
                                                                                                                                                   tGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTgcagc                                                      >  1:1288207/1‑90 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
GGATGGAACCATTCGACAAGGATCCCGCCATGTACAACAGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGCCTATGTCAGCCC  >  NC_002947/6091603‑6091891

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