Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,705,470 A→G 6.4% intergenic (‑162/+15) kdpF ← / ← PP_4162 potassium‑transporting ATPase subunit F/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,705,4700AG 6.4% 268.5 / 6.2 93intergenic (‑162/+15)kdpF/PP_4162potassium‑transporting ATPase subunit F/membrane protein
Reads supporting (aligned to +/- strand):  ref base A (44/43);  new base G (4/2);  total (48/45)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

CCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAACACGCCCAG  >  NC_002947/4705321‑4705616
                                                                                                                                                     |                                                                                                                                                  
ccGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTTTCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaaa                                                                                                                                                    >  2:1132746/1‑150 (MQ=255)
 cGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAg                                                                                                                                                   >  2:2111452/1‑150 (MQ=255)
       atggtttGGGCCGCTTGGGAGAAATTGCCCGGGGGAGAGGTGTAAGGGGGGGAGAGCGGGGGGGGGGGGTGGTTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGt                                                                                                                                             <  1:457343/145‑1 (MQ=255)
         ggCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTa                                                                                                                                           >  1:911418/1‑150 (MQ=255)
          gCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGGGTAAAAATCAGCGTTGTGGTGTGACTGTAAAGGCCTCTTCGCGGTCTTGGACGCGAAAAGGCCGTTag                                                                                                                                          >  1:814598/1‑149 (MQ=255)
            ttGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAg                                                                                                                                        >  1:1525004/1‑150 (MQ=255)
             ggggCGGTTTGGGGGAAATGGTCCGGGGGAGAGGGGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAgg                                                                                                                                       <  2:1525004/149‑1 (MQ=255)
                 ccGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTa                                                                                                                                   <  2:911418/150‑1 (MQ=255)
                   gCTTTGGGGAAATTGTAGGAGGGAAAGGGGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTACAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTAtt                                                                                                                                 <  2:1794844/150‑1 (MQ=255)
                     tttGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGAGACGTTACAGGCGTAttgt                                                                                                                               >  2:1334961/1‑150 (MQ=255)
                          ggAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGACGGCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTgg                                                                                                                          >  1:145294/1‑150 (MQ=255)
                             aaTTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTggcgg                                                                                                                       >  2:240402/1‑150 (MQ=255)
                              aTTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATACAGAATGCGTCAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTGCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTggcggc                                                                                                                      >  2:262665/1‑150 (MQ=255)
                                 gTCCGAGGGAGAGGCGTAAGGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGgctgc                                                                                                                   <  2:395339/150‑1 (MQ=255)
                                    cGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCtgctgc                                                                                                                >  1:1550257/1‑150 (MQ=255)
                                      ggggAGAGGGGTAAAGGGGCGAGATGGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTgcgc                                                                                                              <  2:145294/149‑1 (MQ=255)
                                         gagaGGGGGAAGGGGGCGCGACCGGCTGGGGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGcgccgc                                                                                                           <  2:2291227/150‑1 (MQ=255)
                                         gagaGGCGTAAGGGGGAGGGTGCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGcgccgc                                                                                                           <  2:1995666/150‑1 (MQ=255)
                                            agGCGTAAAGGGGCGCGATCGGGGGGGGGGGGTGGGTAGAAAGAATGCGTAACAATCAGCGCGGAGGTCTGCCTGTAACGGCATCTTCGCGGGGTTGCGTGCGGAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGtt                                                                                                        >  2:1955272/1‑150 (MQ=255)
                                                     ggggCGAGATGGGGGGGGGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtcca                                                                                               <  1:1983009/150‑1 (MQ=255)
                                                     ggggCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTggcg                                                                                                                        >  1:1512192/1‑125 (MQ=255)
                                                     ggggCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTggcg                                                                                                                        <  2:1512192/125‑1 (MQ=255)
                                                     ggggCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtcca                                                                                               >  1:1499316/1‑150 (MQ=255)
                                                      gggCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCATTTCCTGATTAGCGGGGAAAGGGGCGTTACAGGATTACTGATGTGGTGGTCGCCTTGCCGCGTTACGGCtccag                                                                                              >  1:2156090/1‑150 (MQ=255)
                                                        gtgaaagcGGGGGGTGGGGGGGAGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtccagtc                                                                                            <  1:1657474/143‑1 (MQ=255)
                                                        gCGAGAGGGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtccagtc                                                                                            <  1:2015094/150‑1 (MQ=255)
                                                          gagaTCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTgg                                                                                                                          <  2:335813/118‑1 (MQ=255)
                                                          gagaTCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTgg                                                                                                                          >  1:335813/1‑118 (MQ=255)
                                                           agaTCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCtgctgcac                                                                                                              >  1:954797/1‑127 (MQ=255)
                                                           agaTCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCtgctgcac                                                                                                              <  2:954797/129‑3 (MQ=255)
                                                             aTCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCcagca                                                                                       >  1:226135/1‑150 (MQ=255)
                                                              tCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCcagcag                                                                                      >  2:2226810/1‑150 (MQ=255)
                                                              ggggggggTGGGGGTGGAGATAAAGAATGCGTAAAAATCCGCGTTGTGGTCTCCCTGTAACCGCCTCTTCGCGCGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGCGGCGGCTGCTGCGCTGCGTTACGGCTCCAGTCcagcag                                                                                      <  2:1517673/148‑1 (MQ=255)
                                                                 gggggtgggggtGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGACTGCTGGACCGCGTTACGGCTCCAGTCcagcagcag                                                                                   >  1:516158/1‑150 (MQ=255)
                                                                    ggtgggggtgggTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCt                                                                                >  1:1280995/1‑150 (MQ=255)
                                                                     gtgggggtgggTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTg                                                                               <  2:226135/150‑1 (MQ=255)
                                                                       gggggtgggTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTa                                                                             <  2:1499316/150‑1 (MQ=255)
                                                                           gtgggtATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGc                                                                                                                                      <  2:1793551/89‑1 (MQ=255)
                                                                           gtgggtATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGc                                                                                                                                      >  1:1793551/1‑89 (MQ=255)
                                                                            tgggtATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCt                                                                                                                                     <  2:2133494/89‑1 (MQ=255)
                                                                            tgggtATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCt                                                                                                                                     >  1:2133494/1‑89 (MQ=255)
                                                                               gtATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCg                                                                     >  1:1237946/1‑150 (MQ=255)
                                                                               gtATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCg                                                                     >  2:1209980/1‑150 (MQ=255)
                                                                                tataAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGc                                                                    >  1:1305973/1‑150 (MQ=255)
                                                                                 ataAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTggc                                                                                                                         >  2:1622015/1‑96 (MQ=255)
                                                                                 ataAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTggc                                                                                                                         <  1:1622015/96‑1 (MQ=255)
                                                                                 ataAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTATGCTACCGcc                                                                   >  1:1524671/1‑150 (MQ=255)
                                                                                      gAATGCGTAACAATCAGCGTTGTGGTATGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCATTACAGGCTTAGTGTTGTGACGGCTGCTGCGCCGCGTTACGCACCCGGTCCAACAACAAGCCGTAGCCTACCGCCAGcca                                                              >  2:537857/1‑148 (MQ=255)
                                                                                      gAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCaa                                                              >  1:1945393/1‑150 (MQ=255)
                                                                                            gTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGcgcc                                                                                                             >  2:631457/1‑97 (MQ=255)
                                                                                            gTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGcgcc                                                                                                             <  1:631457/97‑1 (MQ=255)
                                                                                            gTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtc                                                                                                 <  1:1898974/109‑1 (MQ=255)
                                                                                            gTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtc                                                                                                 >  2:1898974/1‑109 (MQ=255)
                                                                                            gTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCt                                                                                                  <  2:1521932/108‑1 (MQ=255)
                                                                                            gTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCt                                                                                                  >  1:1521932/1‑108 (MQ=255)
                                                                                             tAAAAATCAGCGTTGTGGTCTGCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGAGTCCCTCCGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGc                                                       <  1:1334961/150‑1 (MQ=39)
                                                                                             tAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTAtt                                                                                                                                 <  2:2286617/76‑1 (MQ=255)
                                                                                             tAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTAtt                                                                                                                                 >  1:2286617/1‑76 (MQ=255)
                                                                                                 aaTCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGa                                                   <  1:1209980/150‑1 (MQ=255)
                                                                                                        gTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTGCGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCa                                                                  <  1:817732/128‑1 (MQ=255)
                                                                                                        gTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTGCGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCa                                                                  >  2:817732/1‑128 (MQ=255)
                                                                                                        gTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAg                                            <  1:2111452/150‑1 (MQ=255)
                                                                                                         ttGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAATGTCGGCCCGATGAAAAgg                                           >  2:837593/1‑150 (MQ=255)
                                                                                                          tgtgGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGc                                          <  1:240402/150‑1 (MQ=255)
                                                                                                             ggTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCaa                                       <  2:1757126/150‑1 (MQ=255)
                                                                                                             ggTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCaa                                       <  1:262665/150‑1 (MQ=255)
                                                                                                              gtctttcTGTAACGGCCTCGTCGCGGGCTTGCCGGCGAAAGGGCCGTTACAGGCTTATTGTTGTGCCGGAGATCGGGCCGCGTTACGGCTCCGGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAAt                                      <  1:519927/144‑1 (MQ=255)
                                                                                                              gTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTAttgttg                                                                                                                             >  2:577837/1‑63 (MQ=255)
                                                                                                              gTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTAttgttg                                                                                                                             <  1:577837/63‑1 (MQ=255)
                                                                                                                 tGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGTCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGaa                                   <  1:2226810/150‑1 (MQ=255)
                                                                                                                 tGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGaa                                   <  1:471569/150‑1 (MQ=255)
                                                                                                                  gCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAg                                  >  1:455177/1‑150 (MQ=255)
                                                                                                                    cTGTAACGGCCCTTTCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGcc                                >  2:239285/1‑150 (MQ=37)
                                                                                                                     tGTAACGTCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCa                               <  1:420619/150‑1 (MQ=255)
                                                                                                                     tGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACGGTCTTGTTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCa                               >  2:450960/1‑150 (MQ=255)
                                                                                                                      gTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCaa                              >  2:1603336/1‑150 (MQ=255)
                                                                                                                       tagcgaccTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCaaa                             >  2:1163030/7‑150 (MQ=255)
                                                                                                                       tAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTa                                                                             >  2:772830/1‑102 (MQ=255)
                                                                                                                       tAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTa                                                                             <  1:772830/102‑1 (MQ=255)
                                                                                                                          cGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgc                          >  2:1117731/1‑150 (MQ=255)
                                                                                                                            gccctttCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGa                                                   >  2:1802293/6‑123 (MQ=37)
                                                                                                                            gccctttCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGa                                                   <  1:1802293/118‑1 (MQ=37)
                                                                                                                            gccctttCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGGAGCCAAAGgcaa                        >  1:1662490/6‑150 (MQ=38)
                                                                                                                            gCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaa                        <  2:1280995/150‑1 (MQ=255)
                                                                                                                                 ttCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGcaagcaag                   <  1:1163030/150‑1 (MQ=255)
                                                                                                                                 ttCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAACCAAAGgcaa                        <  1:1578671/145‑1 (MQ=255)
                                                                                                                                 ttCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAACCAAAGgcaa                        >  2:1578671/1‑145 (MQ=255)
                                                                                                                                      gggCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAgccgcc              <  2:2156090/150‑1 (MQ=255)
                                                                                                                                         cTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGaa           <  1:1124253/150‑1 (MQ=255)
                                                                                                                                          aagCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtcca                                                                                               <  2:724079/63‑1 (MQ=37)
                                                                                                                                          aagCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtcca                                                                                               >  1:724079/3‑65 (MQ=37)
                                                                                                                                              ccGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGaa                                   >  1:2205989/1‑121 (MQ=255)
                                                                                                                                              ccGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGaa                                   <  2:2205989/121‑1 (MQ=255)
                                                                                                                                                  gAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAACACGCCcag  >  1:1954540/1‑150 (MQ=255)
                                                                                                                                                   aaaGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAacac        >  1:1865065/1‑143 (MQ=255)
                                                                                                                                                   aaaGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAacac        <  2:1865065/143‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAACACGCCCAG  >  NC_002947/4705321‑4705616

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