Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,707,364 C→G 9.4% intergenic (‑20/+90) PP_4163 ← / ← PP_4164 hypothetical protein/alpha/beta hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,707,3640CG 9.4% 319.5 / 19.5 104intergenic (‑20/+90)PP_4163/PP_4164hypothetical protein/alpha/beta hydrolase
Reads supporting (aligned to +/- strand):  ref base C (44/50);  new base G (4/6);  total (48/56)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

Only 100 of 108 total aligned reads displayed.

GCTCCGTCGCCGGGCGTTGATCGACCACCTTGGCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATC  >  NC_002947/4707215‑4707513
                                                                                                                                                     |                                                                                                                                                     
gCTCCGTCGCCGGGCGTTGATCGACCACCTTGGCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCACTc                                                                                                                                                       >  1:873859/1‑150 (MQ=255)
   ccGTCGCCGGGCGTTGATCGACCACCTTGGCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAgc                                                                                                                                                    <  2:439038/150‑1 (MQ=255)
    cGTCGCCGGGCGTTGATCGACCACCTTGGCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAgcg                                                                                                                                                   <  2:2255223/150‑1 (MQ=255)
     gTCGCCGGGCGTTGATCGACCACCTTGGCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAgcgc                                                                                                                                                  <  2:629099/150‑1 (MQ=255)
     gTCGCCGGGCGTTGATCGACCACCTTGGCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAgcgc                                                                                                                                                  >  2:2243565/1‑150 (MQ=255)
          cGGGCGTTGATCGACCACCTTGGCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTGAGCGCCgcct                                                                                                                                             <  2:595173/150‑4 (MQ=255)
           gggCGTTGATCGACCACCTTGGCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTGAGCGCCgcctg                                                                                                                                            <  2:2073419/150‑5 (MQ=255)
           gggCGTTGATCGACCACCTTGGCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTGAGCGCCgcctg                                                                                                                                            <  2:202759/150‑5 (MQ=255)
           gggCGTTGATCGACCACCTTGGCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAAtt                                                                                                                                            <  1:810048/150‑1 (MQ=255)
               gTTGATCGACCACCTTGGCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAg                                                                                                                                        <  2:2285787/150‑1 (MQ=255)
                  gATCGACCACCTTGGCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGTCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTc                                                                                                                                     >  1:974418/1‑150 (MQ=255)
                     cGATCACCTGGGCGTCGTGCGGGACCGGTAGGCTGCGGGACAACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAGGGCAGTGCGAATGGAAACAGGCACGCGCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAg                                                                                                                                  <  2:514701/150‑1 (MQ=255)
                      gACCACCTTGGCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTccggg                                                                                                                                 >  2:1974293/1‑147 (MQ=255)
                      gACCACCTTGGCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAgg                                                                                                                                 >  2:377071/1‑150 (MQ=255)
                          accTTGGCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGAt                                                                                                                               >  2:2335800/1‑148 (MQ=255)
                          accTTGGCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGAt                                                                                                                               <  1:2335800/148‑1 (MQ=255)
                           ccTTGGCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTc                                                                                                                            >  1:619831/1‑150 (MQ=255)
                                   tCCTGCGGGACCTGTAGGCTGTCGGCCTATGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGtt                                                                                                                    <  2:1075029/150‑1 (MQ=255)
                                       gCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCGGGATGTCTGGCAGTTGCAt                                                                                                                <  1:1974293/150‑1 (MQ=255)
                                       gCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCAt                                                                                                                <  2:2050444/150‑1 (MQ=255)
                                       gCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCAt                                                                                                                <  2:822127/150‑1 (MQ=255)
                                        cGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATg                                                                                                               <  1:1181396/150‑1 (MQ=255)
                                           gACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGtt                                                                                                            <  2:1526216/150‑1 (MQ=255)
                                           gACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGtt                                                                                                            <  2:615091/150‑1 (MQ=255)
                                                gTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGtctc                                                                                                       >  1:2322818/1‑150 (MQ=255)
                                                gTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGtctc                                                                                                       <  1:1183440/150‑1 (MQ=255)
                                                gTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGtctc                                                                                                       >  1:830057/1‑150 (MQ=255)
                                                 tAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGaa                                                                                                                                         >  1:1279731/1‑115 (MQ=255)
                                                 tAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGaa                                                                                                                                         <  2:1279731/115‑1 (MQ=255)
                                                 tAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTcc                                                                                                      <  1:534935/150‑1 (MQ=255)
                                                      tGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGtt                                                                                                 <  1:236716/150‑1 (MQ=255)
                                                      tGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGtt                                                                                                 >  1:2325957/1‑150 (MQ=255)
                                                      tGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGtt                                                                                                 <  2:313459/150‑1 (MQ=255)
                                                           gCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAAcc                                                                                            >  2:1947635/1‑150 (MQ=255)
                                                            cccACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCGCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAAccc                                                                                           <  1:2243565/150‑1 (MQ=255)
                                                             ccACGACACGGGGGTGGCGAGCACCAGGCAGGCAGTCAAGGCAGTGCGAATGGAAATAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCg                                                                                          >  1:1213311/1‑150 (MQ=255)
                                                             ccACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCg                                                                                          >  2:1596700/1‑150 (MQ=255)
                                                              cACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGc                                                                                         >  1:2264053/1‑150 (MQ=255)
                                                               agtacgcGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGcc                                                                                        >  2:2211659/7‑150 (MQ=255)
                                                                 gACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCaa                                                                                      <  2:830057/150‑1 (MQ=255)
                                                                   cacGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGc                                                                                    >  1:567843/1‑150 (MQ=255)
                                                                     cGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCt                                                                                  >  2:1710207/1‑150 (MQ=255)
                                                                      gggggTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTc                                                                                 <  2:171549/150‑1 (MQ=255)
                                                                        gggTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGa                                                                               >  2:1253089/1‑150 (MQ=255)
                                                                        gggTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCACTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGa                                                                               <  2:873859/150‑1 (MQ=255)
                                                                          gTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTcc                                                                                                      >  2:1255388/1‑125 (MQ=255)
                                                                          gTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTcc                                                                                                      <  1:1255388/125‑1 (MQ=255)
                                                                          gTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGaca                                                                             <  1:1710207/150‑1 (MQ=255)
                                                                           tGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGacac                                                                            <  2:1693634/150‑1 (MQ=255)
                                                                           tGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCAGAGCTGCCAGTCGAATTGCAGCTTGAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGacac                                                                            <  1:1498873/150‑1 (MQ=255)
                                                                           tGGCGAGAACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAGACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCTTGTTCATCCCGCCAAGCGTCGacac                                                                            <  2:400608/150‑1 (MQ=255)
                                                                             gCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCAGAATTGAAGCTCCAGGGTGTCTATCATTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACtt                                                                          <  1:1253089/150‑1 (MQ=255)
                                                                                  cACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCa                                                                                                                       >  1:2345608/1‑100 (MQ=255)
                                                                                  cACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCa                                                                                                                       <  2:2345608/100‑1 (MQ=255)
                                                                                     caggcaggCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTg                                                                  >  2:390829/1‑150 (MQ=255)
                                                                                     caggcaggCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTg                                                                  >  1:498067/1‑150 (MQ=255)
                                                                                      aggcaggcAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAg                                                                                                                                  <  2:833999/85‑1 (MQ=255)
                                                                                      aggcaggcAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAg                                                                                                                                  >  1:833999/1‑85 (MQ=255)
                                                                                       ggcaggcaGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGa                                                                                                                                          <  2:638158/76‑1 (MQ=255)
                                                                                       ggcaggcaGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGa                                                                                                                                          >  1:638158/1‑76 (MQ=255)
                                                                                        gcaggcagCCAAGTCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATTGTTGTTTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCaa                                                               <  1:300922/150‑1 (MQ=255)
                                                                                        gcaggcagCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTc                                                                                 >  2:973736/1‑132 (MQ=255)
                                                                                        gcaggcagCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTc                                                                                 <  1:973736/132‑1 (MQ=255)
                                                                                          aggcagCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGt                                                             >  2:985797/1‑150 (MQ=255)
                                                                                            gcagCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGtt                                                                                                                    >  1:32029/1‑93 (MQ=255)
                                                                                            gcagCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGtt                                                                                                                    <  2:32029/93‑1 (MQ=255)
                                                                                            gcagCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCa                                                           >  1:2003218/1‑150 (MQ=255)
                                                                                             cagCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCaa                                                          <  1:569699/150‑1 (MQ=255)
                                                                                                ccAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCa                                                       >  1:1094824/1‑150 (MQ=255)
                                                                                                ccAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCa                                                       >  2:1408812/1‑150 (MQ=255)
                                                                                                  aaGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCaaa                                                     <  2:2325957/150‑1 (MQ=255)
                                                                                                  aaGGCAGGGCGAATGGAAACAGGCCCGCTCATTGGCGCACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCaaa                                                     <  2:1599336/150‑1 (MQ=255)
                                                                                                     gCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGcc                                                  >  2:536765/1‑150 (MQ=255)
                                                                                                         tGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAaag                                              >  2:1458320/1‑150 (MQ=255)
                                                                                                             aaTGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGc                                          >  2:768272/1‑150 (MQ=255)
                                                                                                               tGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGc                                          <  2:1903620/148‑1 (MQ=255)
                                                                                                               tGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGc                                          >  1:1903620/1‑148 (MQ=255)
                                                                                                               tGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGa                                        >  2:1331811/1‑150 (MQ=255)
                                                                                                                 gAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGaa                                                                                                                                              >  1:754491/1‑46 (MQ=255)
                                                                                                                 gAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGaa                                                                                                                                              <  2:754491/46‑1 (MQ=255)
                                                                                                                     cAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATa                                  <  1:1331811/150‑1 (MQ=255)
                                                                                                                     cAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATa                                  <  2:607482/150‑1 (MQ=255)
                                                                                                                       ggCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACa                                <  1:390829/150‑1 (MQ=255)
                                                                                                                          aCGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACa                                <  2:708847/147‑1 (MQ=255)
                                                                                                                          aCGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACa                                >  1:708847/1‑147 (MQ=255)
                                                                                                                                aTTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGcaac                       <  1:768272/150‑1 (MQ=255)
                                                                                                                                 ttGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAGAGCCAAAGAAGCGATTGATACAGTCTGcaaca                      >  2:1107954/1‑150 (MQ=255)
                                                                                                                                  tGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGcaacaa                     <  2:1661325/150‑1 (MQ=255)
                                                                                                                                   ggCGAACCAGGCGGCGCGCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGGCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAg                    <  2:467408/150‑1 (MQ=255)
                                                                                                                                    gCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCaa                                                          <  2:727474/111‑1 (MQ=255)
                                                                                                                                    gCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCaa                                                          >  1:727474/1‑111 (MQ=255)
                                                                                                                                    gCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAgg                   >  1:248485/1‑150 (MQ=255)
                                                                                                                                       aaCCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCt                <  1:987411/150‑1 (MQ=255)
                                                                                                                                          aattcggcgCTGAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCaa                                                               >  1:2238058/5‑100 (MQ=255)
                                                                                                                                          aattcggcgCTGAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCaa                                                               <  2:2238058/96‑1 (MQ=255)
                                                                                                                                            ttcggcgCTGAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCaaa                                                     >  1:437004/3‑108 (MQ=255)
                                                                                                                                            ttcggcgCTGAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCaaa                                                     <  2:437004/106‑1 (MQ=255)
                                                                                                                                            ttcggcgCTGAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGcc           >  1:1914890/3‑150 (MQ=255)
                                                                                                                                            ttcggcgCTGAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGcc           >  1:929449/3‑150 (MQ=255)
                                                                                                                                                     cAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATc  <  2:415185/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTCCGTCGCCGGGCGTTGATCGACCACCTTGGCGTCCTGCGGGACCGGTAGGCTGCCGGCCCACGACACGGGGGTGGCGAGCACCAGGCAGGCAGCCAAGGCAGTGCGAATGGAAACAGGCACGCTCATTGGCGAACCAGGCGGCGCTCAGCGCCGAATTGAAGCTCCAGGATGTCTGGCAGTTGCATGGTTGTCTCCCTGTTCAACCCGCCAAGCCTCGACACTTGTCCGGTGCAAGTCAAGCAAAGCCAAAGAAGCGATTGATACAGTCTGCAACAAGGGCTGCGCCCTGCGCATC  >  NC_002947/4707215‑4707513

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: