Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,486,594 G→A 11.8% R123Q (CGG→CAG) ‡ yhdZ → amino acid ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,486,5940GA11.8% 593.7 / 55.8 220R123Q (CGG→CAG) ‡yhdZamino acid ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (100/94);  new base A (11/15);  total (111/109)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 220 total aligned reads displayed.

ATGGCGTGGAACTGACCAACGACCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCT  >  NC_002947/1486447‑1486741
                                                                                                                                                   |                                                                                                                                                   
aTGGCGTGGAACTGACCAACGACCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGgc                                                                                                                                                   <  2:4476194/150‑1 (MQ=255)
  ggCGTGGAACTGACCAACGACCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGgcgc                                                                                                                                                 <  1:314553/150‑1 (MQ=255)
    cGTGGAACTGACCAACGACCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCaa                                                                                                                                               >  1:2782054/1‑150 (MQ=255)
     gTGGAACTGACCAACGACCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAg                                                                                                                                              >  1:2562245/1‑150 (MQ=255)
     gTGGAACTGACCAACGACCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAg                                                                                                                                              >  2:483757/1‑150 (MQ=255)
       ggAACTGACCAACGACCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGc                                                                                                                                            >  1:3953945/1‑150 (MQ=255)
          aCTGACCAACGACCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGa                                                                                                                                         <  2:3917126/150‑1 (MQ=255)
          aCTGACCAACGACCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGa                                                                                                                                         <  1:1926276/150‑1 (MQ=255)
              aCCAACGACCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGaa                                                                                                                                     >  2:2405532/1‑150 (MQ=255)
                cAACGACCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAAt                                                                                                                                   >  2:2253272/1‑150 (MQ=255)
                cAACGACCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAAt                                                                                                                                   >  1:3077343/1‑150 (MQ=255)
                 aaCGACCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAAtt                                                                                                                                  >  1:1424046/1‑150 (MQ=255)
                 aaCGACCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAAtt                                                                                                                                  >  1:4989395/1‑150 (MQ=255)
                 aaCGACCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAAtt                                                                                                                                  >  2:1988299/1‑150 (MQ=255)
                       cTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATg                                                                                                                            >  2:2726966/1‑150 (MQ=255)
                        tCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGc                                                                                                                           <  2:1850467/150‑1 (MQ=255)
                          aaGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCAc                                                                                                                         >  2:5292354/1‑150 (MQ=255)
                           aGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACt                                                                                                                        <  1:2405532/150‑1 (MQ=255)
                                 tCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTgg                                                                                                                  <  1:2111677/150‑1 (MQ=255)
                                  cGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGa                                                                                                                 <  1:527523/150‑1 (MQ=255)
                                  cGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGa                                                                                                                 <  1:617529/150‑1 (MQ=255)
                                      gCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCgt                                                                                                             <  1:1011138/150‑1 (MQ=255)
                                        cATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCgtgt                                                                                                           >  1:5351942/1‑150 (MQ=255)
                                        cATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCgtgt                                                                                                           <  2:2754926/150‑1 (MQ=255)
                                         aTCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGtgtg                                                                                                          >  2:1346611/1‑150 (MQ=255)
                                          tCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGTCATGCACTACCTGGAGCGTGTGc                                                                                                         >  2:5250487/1‑150 (MQ=255)
                                               cgcgAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGAtt                                                                                                    >  2:4106283/1‑150 (MQ=255)
                                                 cgAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTcc                                                                                                  <  2:1183998/150‑1 (MQ=255)
                                                  gAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCg                                                                                                 <  1:3922773/150‑1 (MQ=255)
                                                    ggTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGa                                                                                               <  1:1730323/150‑1 (MQ=255)
                                                       gggCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCa                                                                                            >  1:1651751/1‑150 (MQ=255)
                                                         gCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAgg                                                                                          <  2:5289484/150‑1 (MQ=255)
                                                          cATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAAtt                                                                                                                                  >  2:2355995/1‑109 (MQ=255)
                                                                gtTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACaa                                                                                   <  1:1605028/150‑1 (MQ=255)
                                                                gtTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACaa                                                                                   <  1:5184369/150‑1 (MQ=255)
                                                                  tCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGc                                                                                                                           >  2:5028100/1‑108 (MQ=255)
                                                                  tCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGc                                                                                                                           <  1:5028100/108‑1 (MQ=255)
                                                                    cAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTAc                                                                               >  2:4581696/1‑150 (MQ=255)
                                                                     aGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACaa                                                                                   <  2:1643770/145‑1 (MQ=255)
                                                                     aGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTAcc                                                                              >  1:3384160/1‑150 (MQ=255)
                                                                     aGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTAcc                                                                              >  2:1701512/1‑150 (MQ=255)
                                                                     aGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTAcc                                                                              >  1:3613275/1‑150 (MQ=255)
                                                                      gCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTAccc                                                                             <  1:1972264/150‑1 (MQ=255)
                                                                          ttCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGc                                                                                                                                            <  1:474733/83‑1 (MQ=255)
                                                                          ttCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGc                                                                         >  1:3604464/1‑150 (MQ=255)
                                                                          ttCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGc                                                                         <  1:1164506/150‑1 (MQ=255)
                                                                            cAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCa                                                                       >  2:4256493/1‑150 (MQ=255)
                                                                              aCCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGa                                                                                               >  1:2687330/1‑124 (MQ=255)
                                                                              aCCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGa                                                                                               <  2:2687330/124‑1 (MQ=255)
                                                                               ccTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACAGAAGTACCCTGGCCAGtt                                                                    <  1:1795249/150‑1 (MQ=255)
                                                                               ccTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGtt                                                                    <  1:3662842/150‑1 (MQ=255)
                                                                                     cccGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTgg                                                              <  2:2558926/150‑1 (MQ=255)
                                                                                      ccGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTggc                                                             <  2:2784418/150‑1 (MQ=255)
                                                                                       cGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTggcg                                                            >  2:4688578/1‑150 (MQ=255)
                                                                                         cACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTTGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTggcggc                                                          >  2:1664247/1‑150 (MQ=255)
                                                                                         cACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTggcggc                                                          <  2:3446345/150‑1 (MQ=255)
                                                                                           ccTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGAt                                                                                                     <  1:4699114/105‑1 (MQ=255)
                                                                                             tGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCcagc                                                      <  1:4486797/150‑1 (MQ=255)
                                                                                             tGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCcagc                                                      >  2:4155036/1‑150 (MQ=255)
                                                                                              gAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACaa                                                                                   >  1:3847947/1‑120 (MQ=255)
                                                                                                   tCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCagcagcagc                                                <  2:394710/150‑1 (MQ=255)
                                                                                                   tCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCagcagcagc                                                >  1:1861787/1‑150 (MQ=255)
                                                                                                        gagaACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGtgtg                                           >  2:3445797/1‑150 (MQ=255)
                                                                                                        gagaACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGtgtg                                           <  2:1442773/150‑1 (MQ=255)
                                                                                                         aaaaCTGCACCCTGGCACCGATTTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAGGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTgg                                          <  1:521162/148‑1 (MQ=255)
                                                                                                          gaACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGc                                         >  2:4160728/1‑150 (MQ=255)
                                                                                                            aCTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCGCACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCa                                       <  2:3609508/150‑1 (MQ=255)
                                                                                                             cTGCACCCTGGCACCGATTTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGAt                                      <  1:4841993/150‑1 (MQ=255)
                                                                                                              tGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCAtt                                     >  1:1562365/1‑150 (MQ=255)
                                                                                                              tGCACCCTGGCACCGATTTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGAtt                                     >  2:532042/1‑150 (MQ=255)
                                                                                                              tGCACCCTGGCACCGATTTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATAGCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGAtt                                     <  1:2861328/150‑1 (MQ=255)
                                                                                                                cACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCg                                                                                                        <  2:3882119/81‑1 (MQ=255)
                                                                                                                cACCCTGGCACCGATCTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTATCAGGCGGCCAGCAGCAGCGTGTGGCGATTGc                                   <  2:3200571/150‑1 (MQ=255)
                                                                                                                 aCCCTGGCACCGATCTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGcc                                  <  1:5233219/150‑1 (MQ=255)
                                                                                                                  cccTGGCACCGATCTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGccc                                 >  1:4624755/1‑150 (MQ=255)
                                                                                                                   ccTGGCACCGATCTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCcg                                >  2:615671/1‑150 (MQ=255)
                                                                                                                   ccTGGCACCGATCTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCcg                                >  2:1617093/1‑150 (MQ=255)
                                                                                                                     tGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGc                                                                                         <  2:485061/91‑1 (MQ=255)
                                                                                                                     tGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGc                                                                                         <  2:484871/91‑1 (MQ=255)
                                                                                                                     tGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGc                                                                                         >  1:484886/1‑91 (MQ=255)
                                                                                                                     tGGCACCGATCTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCcgtg                              <  2:3894621/150‑3 (MQ=255)
                                                                                                                         gcCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCgg                                                                                                >  1:3511830/1‑80 (MQ=255)
                                                                                                                         gcCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGtt                          <  1:1259623/150‑1 (MQ=255)
                                                                                                                           cGATCTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCcgtgcgcttt                        <  1:4149224/150‑5 (MQ=255)
                                                                                                                            gATTTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCcgtgcgctttg                       <  1:1086267/150‑6 (MQ=255)
                                                                                                                            gATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTtgtg                       >  2:3370643/1‑150 (MQ=255)
                                                                                                                             atttgGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGc                      >  2:3210377/4‑150 (MQ=255)
                                                                                                                             aTGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCTCGCGTTGTGc                      <  2:4733426/150‑1 (MQ=255)
                                                                                                                               gtgGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCAt                    <  2:2801310/150‑1 (MQ=255)
                                                                                                                               gtgGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCAt                    <  1:2021056/150‑1 (MQ=255)
                                                                                                                                   gtacgcAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGaac                >  2:408399/4‑149 (MQ=255)
                                                                                                                                   gTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCcgcgc                             <  1:2593383/137‑1 (MQ=255)
                                                                                                                                        cAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGaa           <  1:921075/150‑1 (MQ=255)
                                                                                                                                           gATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGAt        >  1:4363884/1‑150 (MQ=255)
                                                                                                                                            aTGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATc       <  1:174611/150‑1 (MQ=255)
                                                                                                                                              gCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCAt     <  1:2374593/150‑1 (MQ=255)
                                                                                                                                               ccACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATg    <  1:4155036/150‑1 (MQ=255)
                                                                                                                                                 aCGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGAtt                                                                                                    >  2:4995177/1‑52 (MQ=255)
                                                                                                                                                 aCGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTc                                                                 <  1:459504/87‑1 (MQ=255)
                                                                                                                                                 aCGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTc                                                                 >  2:459504/1‑87 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATGGCGTGGAACTGACCAACGACCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCT  >  NC_002947/1486447‑1486741

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