Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,406,056 C→G 8.3% A184A (GCG→GCC PP_4750 ← amino acid ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,406,0560CG 8.3% 199.9 / 7.3 67A184A (GCG→GCCPP_4750amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (26/34);  new base G (2/3);  total (29/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 = 3.59e-01

AGCACAAAGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGA  >  NC_002947/5405908‑5406198
                                                                                                                                                    |                                                                                                                                              
aGCACAAAGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTCCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCAccgccgcc                                                                                                                                                >  5:235591/1‑149 (MQ=255)
  cacaAAGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACcgccgccgc                                                                                                                                              >  3:132568/1‑149 (MQ=255)
   acaAAGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCcgcg                                                                                                                                             <  5:197527/149‑1 (MQ=255)
   acaAAGCACTTCCACAACTCCGGCAACCAAAGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCCACAGCCAACGAGAGAACAAGGCAACTGGCACGCTGAGCACCAGGTACACCAccgcagccgag                                                                                                                                             >  6:82152/1‑147 (MQ=255)
     aaaGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGgt                                                                                                                                           >  8:209062/1‑149 (MQ=255)
         cccTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCGGCGGCGgtggtg                                                                                                                                       <  2:84274/149‑4 (MQ=255)
         cccTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAg                                                                                                                                       <  7:78329/149‑1 (MQ=255)
          ccTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGa                                                                                                                                      >  6:27508/1‑149 (MQ=255)
          ccTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGa                                                                                                                                      >  2:56393/1‑149 (MQ=255)
           cTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGgt                                                                                                                                           >  3:3570/1‑143 (MQ=255)
           cTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGgt                                                                                                                                           <  4:3570/143‑1 (MQ=255)
              ccACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGc                                                                                                                                   <  3:108674/148‑1 (MQ=255)
                    cTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCGGCGGCGgtggtgtacc                                                                                                                                   <  6:126917/142‑8 (MQ=255)
                    cTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCGGCGGCGgtggtgtacc                                                                                                                                   >  5:126917/1‑135 (MQ=255)
                         gCAACGATTTCATGCACTTCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGAAAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTa                                                                                                                                        <  4:230154/132‑1 (MQ=255)
                          cAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTCCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGaaa                                                                                                                      <  6:235591/149‑1 (MQ=255)
                                    aTCCACATCTTCAGCAGTTCTCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACGAGGTACACCACCGACGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCtgt                                                                                                            <  5:82152/149‑1 (MQ=255)
                                       cacaTCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGa                                                                                                          >  6:126673/1‑148 (MQ=255)
                                       cacaTCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGa                                                                                                          >  2:166475/1‑148 (MQ=255)
                                             gTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTg                                                                                                    <  7:184138/148‑1 (MQ=255)
                                                        gCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGaaa                                                                                        <  1:166475/149‑1 (MQ=255)
                                                          tCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGc                                                                                      >  5:239419/1‑149 (MQ=255)
                                                            aTCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGc                                                                                      <  8:35736/147‑1 (MQ=255)
                                                            aTCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGc                                                                                      <  4:29012/147‑1 (MQ=255)
                                                               gggCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGAt                                                                                  <  6:210993/148‑1 (MQ=255)
                                                                  cAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCTTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCa                                                                              >  6:200175/1‑149 (MQ=255)
                                                                     tCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTATGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCAcc                                                                            >  8:123770/1‑148 (MQ=255)
                                                                        tctGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGAt                                                                                  <  7:109517/139‑1 (MQ=255)
                                                                        tctGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGAt                                                                                  >  8:109517/1‑139 (MQ=255)
                                                                                 aCGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAAtcgtc                                                               >  5:189189/1‑149 (MQ=255)
                                                                                     ggCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAgg                                                           >  8:247130/1‑149 (MQ=255)
                                                                                        cGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCtgtg                                                                                                           <  4:113929/98‑1 (MQ=255)
                                                                                        cGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCtgtg                                                                                                           >  3:113929/1‑98 (MQ=255)
                                                                                               cAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGcc                                                   >  1:123319/1‑147 (MQ=255)
                                                                                                aaCGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAg                                                >  2:58188/1‑149 (MQ=255)
                                                                                                 aCGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAAc                                                                         <  8:190622/123‑1 (MQ=255)
                                                                                                 aCGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAAc                                                                         >  7:190622/1‑123 (MQ=255)
                                                                                                   gagagaACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTcc                                             >  5:116686/1‑149 (MQ=255)
                                                                                                     gagaACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCtt                                           <  6:239419/149‑1 (MQ=255)
                                                                                                     gagaACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCtt                                           <  7:209062/149‑1 (MQ=255)
                                                                                                     gagaACAAGGCCACTGGCACGCTGAGCAACTGGTACACCACCGCCGCAGCTGTGTAGACGCTCAGCCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGACTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCt                                            >  1:22231/1‑148 (MQ=255)
                                                                                                      agaACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCTGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTg                                          <  1:180495/149‑1 (MQ=255)
                                                                                                      agaACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCtt                                           <  7:123770/148‑1 (MQ=255)
                                                                                                        aaCAAGGCCTCTGGCACGCTGAGCACCAGGTAGACGACGGCCGCCGCGGTGTAGACTCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAATGGTTAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAg                                        <  2:22231/149‑1 (MQ=255)
                                                                                                         aCAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCTTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGc                                       <  5:200175/149‑1 (MQ=255)
                                                                                                         aCAACGCCCCTGGCACGCCGAGCACCAGGTACCCCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGc                                       <  4:166750/149‑1 (MQ=255)
                                                                                                          cAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCa                                      <  1:58188/149‑1 (MQ=255)
                                                                                                           aaGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAg                                     >  1:89/1‑149 (MQ=255)
                                                                                                            aGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAgg                                    <  5:126673/149‑1 (MQ=255)
                                                                                                             ggCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGc                                   >  5:71091/1‑149 (MQ=255)
                                                                                                                 aCTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAAt                               >  4:27803/1‑149 (MQ=255)
                                                                                                                  cTGGCACGCTGAGCACCAGGTACACCACCGCCTCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATg                              <  6:189189/149‑1 (MQ=255)
                                                                                                                     gCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCt                                                                                     >  4:239125/1‑91 (MQ=255)
                                                                                                                      cACCCTCAGCACCACGAACACCACCGCCGCTGCGTTGTAGACGCTCAGGCCCTGAAAGGATTGCTGTGAAACCTGGTCAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGt                          <  8:47430/149‑1 (MQ=255)
                                                                                                                          cTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGt                      <  5:27508/149‑1 (MQ=255)
                                                                                                                               cACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCaa                                                                    >  7:130228/1‑98 (MQ=255)
                                                                                                                               cACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCaa                                                                    <  8:130228/98‑1 (MQ=255)
                                                                                                                               cACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCt                  >  6:40016/1‑148 (MQ=255)
                                                                                                                               cACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTg                 <  1:56393/149‑1 (MQ=255)
                                                                                                                               cACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTg                 >  3:195170/1‑149 (MQ=255)
                                                                                                                                 ccAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCa                                      <  8:140731/126‑1 (MQ=255)
                                                                                                                                 ccAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCa                                      >  7:140731/1‑126 (MQ=255)
                                                                                                                                     gTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTgg           <  3:122502/149‑1 (MQ=255)
                                                                                                                                     gTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTgg           <  3:27803/149‑1 (MQ=255)
                                                                                                                                     cgtctacaccGCGGCGGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTgg           >  8:85332/7‑149 (MQ=255)
                                                                                                                                     cgtctacaccGCGGCGGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTgg           <  7:85332/143‑1 (MQ=255)
                                                                                                                                         accaccGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAAt       <  4:132568/149‑1 (MQ=255)
                                                                                                                                           caccGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGa                                                 <  3:227657/105‑1 (MQ=255)
                                                                                                                                           caccGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGa                                                 >  4:227657/1‑105 (MQ=255)
                                                                                                                                           caccGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATgg     <  7:247130/149‑1 (MQ=255)
                                                                                                                                            accGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGc    <  5:40016/149‑1 (MQ=255)
                                                                                                                                             ccgccgccgCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCg   <  2:123319/149‑1 (MQ=255)
                                                                                                                                              cgccgccgcGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGa  <  6:116686/149‑1 (MQ=255)
                                                                                                                                               gccgccgcGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCg                           >  6:80817/1‑123 (MQ=255)
                                                                                                                                               gccgccgcGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCg                           <  5:80817/123‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
AGCACAAAGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGA  >  NC_002947/5405908‑5406198

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: