Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,460,988 C→T 5.9% S35S (AGC→AGT)  PP_4797 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,460,9880CT 5.9% 206.5 / 2.4 68S35S (AGC→AGT) PP_4797hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/39);  new base T (2/2);  total (27/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

CCTACAGTGCGCCCGAAACTCACCCAACCCGGGAAGCCCACCATGAGCAAAAAGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAGCGCAGGCATGGTATGGTCGTCACCTGACCTGCAACTCTGGATCTGTGCATGCTCTTCCGTCTTATCCCCCGCTGGGA  >  NC_002947/5460842‑5461137
                                                                                                                                                  |                                                                                                                                                     
ccTACAGTGCGCCCGAAACTCACCCAACCCGGGAAGCCCACCATGAGCAAAAAGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTAc                                                                                                                                                    <  1:588686/150‑1 (MQ=255)
  tACAGTGCGCCCGAAACTCACCCAACCCGGGAAGCCCACCATGAGCAAAAAGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTAccg                                                                                                                                                  >  2:495697/1‑150 (MQ=255)
         cgcCCGAAACTCACCCAACCCGGGAAGCCCACCATGAGCAAAAAGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCCCTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGaa                                                                                                                                           <  1:292794/150‑1 (MQ=255)
         cgcCCGAAACTCACCCAACCCGGGAAGCCCACCATGAGCAAAAAGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGaa                                                                                                                                           <  1:2311874/150‑1 (MQ=255)
          gcCCGAAACTCACCCAACCCGGGAAGCCCACCATGAGCAAAAAGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAg                                                                                                                                          <  2:1333337/150‑1 (MQ=255)
           cccGAAACTCACCCAACCCGGGAAGCCCACCATGAGCAAAAAGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGc                                                                                                                                         <  1:978495/150‑1 (MQ=255)
            ccGAAACTCACCCAACCCGGGAAGCCCACCATGAGCAAAAAGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGcc                                                                                                                                        <  1:1975822/150‑1 (MQ=255)
             cGAAACTCACCCAACCCGGGAAGCCCACCATGAGCAAAAAGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCt                                                                                                                                       >  2:1222831/1‑150 (MQ=255)
              gAAACTCACCCAACCCGGGAAGCCCACCATGAGCAAAAAGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCtt                                                                                                                                      <  2:1076064/150‑1 (MQ=255)
              gAAACTCACCCAACCCGGGAAGCCCACCATGAGCAAAAAGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCtt                                                                                                                                      <  2:1610658/150‑1 (MQ=255)
              gAAACTCACCCAACCCGGGAAGCCCACCATGAGCAAAAAGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCtt                                                                                                                                      <  1:1745134/150‑1 (MQ=255)
               aaaCTCACCCAACCCGGGAAGCCCACCATGAGCAAAAAGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTc                                                                                                                                     <  1:1008091/150‑1 (MQ=255)
               aaaCTCACCCAACCCGGGAAGCCCACCATGAGCAAAAAGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTc                                                                                                                                     <  1:2337706/150‑1 (MQ=255)
                  cTCACCCAACCCGGGAAGCCCACCATGAGCAAAAAGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCaa                                                                                                                                  >  2:507064/1‑150 (MQ=255)
                   tCACCCAACCCGGGAAGCCCACCATGAGCAAAAAGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAAt                                                                                                                                 <  2:1954208/150‑1 (MQ=255)
                   tCACCCAACCCGGGAAGCCCACCATGAGCAAAAAGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAAt                                                                                                                                 >  1:1117188/1‑150 (MQ=255)
                    cACCCAACCCGGGAAGCCCACCATGAGCAAAAAGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATc                                                                                                                                <  2:2032222/150‑1 (MQ=255)
                      cccAACCCGGGAAGCCCACCATGAGCAAAAAGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCga                                                                                                                              <  2:868221/150‑1 (MQ=255)
                      cccAACCCGGGAAGCCCACCATGAGCAAAAAGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCga                                                                                                                              <  2:459768/150‑1 (MQ=255)
                        cAACCCGGGAAGCCCACCATGAGCAAAAAGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTCCCAATctctc                                                                                                                            <  1:605630/150‑5 (MQ=255)
                            ccGGGAAGCCCACCATGAGCAAAAAGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGAtt                                                                                                                        <  1:1125376/150‑1 (MQ=255)
                                 aaGCCCACCATGAGCAAAAAGCCGAAAAAGCCCGGCCCCAACAAGGCCAAGTCCATCCTTGCCCAGCCACTGTTCCGCTGCCGCCCGGAACAACCGAACAATGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAg                                                                                                                   <  2:987977/150‑1 (MQ=255)
                                      caccaTGAGCAAAAAGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTc                                                                                                              >  2:826837/1‑150 (MQ=255)
                                         caTGAGCAAAAAGCCGAACAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCGCTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTa                                                                                                           <  1:710823/150‑1 (MQ=255)
                                              gCAAAAAGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCt                                                                                                                                       >  2:1754904/1‑117 (MQ=255)
                                              gCAAAAAGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCt                                                                                                                                       <  1:1754904/117‑1 (MQ=255)
                                                aaaaaGCCGAAAAAGCACGGCCCCAAAAATGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCCGGCACAACCGAACAGGGGCAAAGGCAGCTACCGCCGCGAGGCCTTCCAATCGAGAGATTGGGGGGCTTCGTACTTTTTg                                                                                                    >  1:901684/1‑150 (MQ=255)
                                                   aaGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCCATATATCCAGTTGGAGGCTTCCTACTTTTTGGCt                                                                                                 <  1:1814450/150‑1 (MQ=255)
                                                   aaGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACttttt                                                                                                     <  1:47876/146‑1 (MQ=255)
                                                   aaGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACttttt                                                                                                     >  2:47876/1‑146 (MQ=255)
                                                        gAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATg                                                                                            >  2:693076/1‑150 (MQ=255)
                                                        gAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATg                                                                                            >  1:9807/1‑150 (MQ=255)
                                                                aCGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATc                                                                                    <  2:1117188/150‑1 (MQ=255)
                                                                      ccAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGATGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCCTGAAAGCATCGCAAgg                                                                              <  2:308594/150‑2 (MQ=255)
                                                                        aacaaGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAgcgc                                                                            >  2:1830005/1‑150 (MQ=255)
                                                                             ggCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAGCTCATGcc                                                                       <  2:1416687/150‑2 (MQ=255)
                                                                              gCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCgaga                                                                                                                            >  1:1188625/1‑96 (MQ=255)
                                                                              gCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCgaga                                                                                                                            <  2:1188625/96‑1 (MQ=255)
                                                                              gCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAGCGCAGGCat                                                                      <  1:693076/150‑1 (MQ=255)
                                                                               ccAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAGCGCAGGCatg                                                                     >  1:321394/1‑150 (MQ=255)
                                                                                        tCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAGCGCAGGCatggtatg                                                                <  2:2105318/146‑1 (MQ=255)
                                                                                        tCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAGCGCAGGCatggtatg                                                                >  1:2105318/1‑146 (MQ=255)
                                                                                          gTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGTTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAGCGCAGGCATGGTATGgtcgtc                                                          >  2:1036269/1‑150 (MQ=255)
                                                                                            tGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGTTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAGCGCAGGCATGGTATGGTCGTCac                                                        <  1:1036269/150‑1 (MQ=255)
                                                                                            tGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAGCGCAGGCATGGTATGGTCGTCac                                                        <  2:598349/150‑1 (MQ=255)
                                                                                                 aGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAGCGCAGGCATGGTATGGTCGTCacctgac                                                   >  2:998794/1‑150 (MQ=255)
                                                                                                         ttCCGCTGCCGCCAGGAACAACCGAATAAGGTCAAAGGCAGCTACCGCGGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAGCGCAGGCATGGTATGGTCGTCACCTGACCTGCAAct                                           >  1:1221032/1‑150 (MQ=255)
                                                                                                         ttCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCga                                                                                                                              <  2:1076223/67‑1 (MQ=255)
                                                                                                         ttCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCga                                                                                                                              >  1:1076223/1‑67 (MQ=255)
                                                                                                           ccGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAGCGCAGGCATGGTATGGTCGTCACCTGACCTGCAActct                                         <  1:998794/150‑1 (MQ=255)
                                                                                                              cTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGaa                                                                                                                                           <  2:17099/49‑1 (MQ=255)
                                                                                                              cTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGaa                                                                                                                                           >  1:17099/1‑49 (MQ=255)
                                                                                                              cTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAGCGCAGGCATGGTATGGTCGTCACCTGACCTGCAACTCTGGa                                      <  1:507064/150‑1 (MQ=255)
                                                                                                              cTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAGCGCAGGCATGGTATGGTCGTCACCTGACCTGCAACTCTGGa                                      <  1:1222831/150‑1 (MQ=255)
                                                                                                                 ccgccAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGGAAGCATCGCAAGCGCAGGCATGGTATGGTCGTCACCTGACCTGCAACTCTGGATCt                                   <  2:9807/150‑1 (MQ=255)
                                                                                                                            aacCGAACAAGGGCAAAGGCAGTTACCGCCGCGAAGCCTTCCAATCGATAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAGCGCAGGCATGGTATGGTCGTCACCTGACCTGCAACTCTGGATCTGTGCAt                             <  1:1028824/145‑1 (MQ=255)
                                                                                                                            aacCGAACAAGGGCAAAGGCAGTTACCGCCGCGAAGCCTTCCAATCGATAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAGCGCAGGCATGGTATGGTCGTCACCTGACCTGCAACTCTGGATCTGTGCAt                             >  2:1028824/1‑145 (MQ=255)
                                                                                                                               cGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAGCGCAGGCATGGTATGGTCGTCACCTGACCTGCAACTCTGGATCTGTGCATGCTCTTcc                     >  2:1758166/1‑150 (MQ=255)
                                                                                                                                gAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAGCGCAGGCATGGTATGGTCGTCACCTGACCTGCAACTCTGGATCTGTGCATGCTCTTCCg                    >  1:798135/1‑150 (MQ=255)
                                                                                                                                    aaGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAGCGCAg                                                                          <  2:1602590/92‑1 (MQ=255)
                                                                                                                                    aaGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAGCGCAg                                                                          >  1:1602590/1‑92 (MQ=255)
                                                                                                                                    aaGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAGCGCAGGCATGGTATGGTCGTCACCTGACCTGCAACTCTGGATCTGTGCATGCTCTTCCGTCtt                <  1:495697/150‑1 (MQ=255)
                                                                                                                                      gggCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAGCGCAGGCATGGTATGGTCGTCACCTGACCTGCAACTCTGGATCTGTGCATGCTCTTCCGt                   >  1:1059521/1‑145 (MQ=255)
                                                                                                                                      gggCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCATACTTTTTGGCTGCATGAAAGCATCGCAAGCGCAGGCATGGTATGGTCGTCACCTGACCTGCAACTCTGGATCTGTGCATGCTCTTCCGt                   <  2:1059521/145‑1 (MQ=255)
                                                                                                                                       ggCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAGCGCAGGCATGGTATGGTCGTCACCTGACCTGCAACTCTGGATCTGTGCATGCTCTTCCGTCTTATc             >  1:1437239/1‑150 (MQ=255)
                                                                                                                                          aaaGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGGGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAGCGCAGGCATGGTATGGTCGTCacctga                                                    >  1:1970222/1‑108 (MQ=255)
                                                                                                                                          aaaGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGGGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAGCGCAGGCATGGTATGGTCGTCacctga                                                    <  2:1970222/108‑1 (MQ=255)
                                                                                                                                            aGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAGCGCAGGCATGGTATGGTCGTCACCTGACCTGCAACTCTGGATCTGTGCATGCTCTTCCGTCTTATCCCCCg        >  1:990515/1‑150 (MQ=255)
                                                                                                                                                  cTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAGCGCAGGCATGGTCTGGTCGTCCCCTGACCTGCAACGCTGGATCTGTGCATGCTCTTCCGTCTTATCCCCCGCTGGGa  >  1:1687438/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CCTACAGTGCGCCCGAAACTCACCCAACCCGGGAAGCCCACCATGAGCAAAAAGCCGAAAAAGCACGGCCCCAACAAGGCCAAGTCCATCGTTGCCCAGCCACTGTTCCGCTGCCGCCAGGAACAACCGAACAAGGGCAAAGGCAGCTACCGCCGCGAAGCCTTCCAATCGAGAGATTGGGAGGCTTCCTACTTTTTGGCTGCATGAAAGCATCGCAAGCGCAGGCATGGTATGGTCGTCACCTGACCTGCAACTCTGGATCTGTGCATGCTCTTCCGTCTTATCCCCCGCTGGGA  >  NC_002947/5460842‑5461137

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: ‑