Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,896,444 T→C 6.4% intergenic (+32/‑176) PP_4304 → / → sbp‑I VIC family cation transporter/sulfate ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,896,4440TC 6.4% 258.8 / 6.4 93intergenic (+32/‑176)PP_4304/sbp‑IVIC family cation transporter/sulfate ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (39/48);  new base C (3/3);  total (42/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.17e-01

CAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGC  >  NC_002947/4896295‑4896589
                                                                                                                                                     |                                                                                                                                                 
cAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAt                                                                                                                                                   >  1:204947/1‑150 (MQ=255)
cAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAt                                                                                                                                                   >  2:1785092/1‑150 (MQ=255)
cAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAt                                                                                                                                                   >  2:569495/1‑150 (MQ=255)
 aaTGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATc                                                                                                                                                  >  1:1072874/1‑150 (MQ=255)
   gggcATGCGTGGGGGGCAGCGGCAGCATGGTTGCGCGACGTGCCACGAGAAGGCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCg                                                                                                                                                <  1:94354/147‑1 (MQ=255)
     ccATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCt                                                                                                                                              <  2:880340/150‑1 (MQ=255)
      cATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCTGCCTCTTCGCGGGTAAATCCGCTc                                                                                                                                             <  2:607621/150‑1 (MQ=255)
      cATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTc                                                                                                                                             <  2:1072874/150‑1 (MQ=255)
          cGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTAc                                                                                                                                         <  1:740934/150‑1 (MQ=255)
             ggggAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAgg                                                                                                                                      >  2:977854/1‑150 (MQ=255)
                         cagcagGACTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTGCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTgc                                                                                                                          <  1:107182/150‑1 (MQ=255)
                              tGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTgcgc                                                                                                                        >  2:1880902/1‑147 (MQ=255)
                              tGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTgcgc                                                                                                                        <  1:1880902/147‑1 (MQ=255)
                              tGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGGACCGTGCGCGGTccc                                                                                                                  >  2:1584961/1‑150 (MQ=38)
                               gATTGCCCGACCTGCCACAAGAAGACCCAAGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTc                                                                                                                    <  2:204947/150‑1 (MQ=255)
                                    cccGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCg                                                                                                               >  1:1114890/1‑150 (MQ=255)
                                    cccGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACGGGGGACCGCGCACGGTCCgcgcctg                                                                                                            <  2:383602‑M1/150‑8 (MQ=255)
                                       gACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAg                                                                                                            >  2:873383/1‑150 (MQ=255)
                                       gACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAg                                                                                                            <  1:873383/150‑1 (MQ=255)
                                        aCCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATGGCGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGGACCGTGCGCGGTCCCCCGTAGGAGc                                                                                                        <  1:1584961/150‑1 (MQ=38)
                                         ccTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGa                                                                                                          >  1:700028/1‑150 (MQ=255)
                                          cTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCgcgcctgtaggag                                                                                                      >  2:158523‑M1/1‑137 (MQ=255)
                                           tGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGGACCGCGCACGGTCCgcgcctgtaggagc                                                                                                     <  1:158523‑M1/150‑15 (MQ=255)
                                              cacaCGAAGACACACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCggg                                                                                                     <  1:959682/150‑1 (MQ=255)
                                              cacaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCggg                                                                                                     <  1:578012/150‑1 (MQ=255)
                                              cacaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGGACCGCGCACGGTCCgcgcctgtaggagcgga                                                                                                  <  1:74714‑M1/150‑18 (MQ=255)
                                               acaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGt                                                                                                    >  2:1616901/1‑150 (MQ=255)
                                                caAGAAGACCCACGTGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGtt                                                                                                   <  1:272414/150‑1 (MQ=255)
                                                 aagaagACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTa                                                                                                             <  2:1277776/139‑1 (MQ=255)
                                                 aagaagACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTa                                                                                                             >  1:1277776/1‑139 (MQ=255)
                                                 aagaagACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGttt                                                                                                  >  2:1465665/1‑150 (MQ=255)
                                                  tggagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGAGGGTAAATCCGCTCCTACAGGCGCGGACAGTGCGCGGTCCACCGTAGAAGCGGGTGTa                                                                                                 >  1:1523596/4‑150 (MQ=255)
                                                   gaagcCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTAc                                                                                                <  2:97593/150‑1 (MQ=255)
                                                     agaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTAccc                                                                                              <  1:1177776/150‑1 (MQ=255)
                                                      gaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCcg                                                                                             >  2:980792/1‑150 (MQ=255)
                                                         ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCg                                                                                                       >  1:1323813/1‑137 (MQ=255)
                                                         ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCg                                                                                                       <  2:1323813/137‑1 (MQ=255)
                                                         ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGa                                                                                          >  2:1862311/1‑150 (MQ=255)
                                                           aCGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAag                                                                                        >  2:1358623/1‑150 (MQ=255)
                                                             gAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAagag                                                                                      <  1:1373806/150‑1 (MQ=255)
                                                             gAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAagag                                                                                      <  1:1616901/150‑1 (MQ=255)
                                                              aGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAgg                                                                                     <  1:1862311/150‑1 (MQ=255)
                                                                 tgggggCGGCGTCGTGTTCTCGCTGTGGGAATGGGTTGTTTCGCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGGCCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACGCGCGAAGAGGCCg                                                                                  <  2:45912/149‑1 (MQ=39)
                                                                   ggggCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGGTTACCCGCGAAGAGGCCgga                                                                                >  1:45912/1‑149 (MQ=255)
                                                                        ggcgTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCcg                                                                                             <  2:1586330/132‑1 (MQ=255)
                                                                        ggcgTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCcg                                                                                             >  1:1586330/1‑132 (MQ=255)
                                                                           gttctgttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGGGGACCGCGCACGGTCCgcgcctgtaggagcggatttacccgcgaagaggccggcgcagg                                                                        <  2:96359‑M1/150‑44 (MQ=255)
                                                                           gttctgttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGGGGACCGCGCACGGTCCgcgcctgtaggagcggatttacccgcgaagaggccggcgcagg                                                                        >  1:96359‑M1/1‑107 (MQ=255)
                                                                           gttctgttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAg                                                                        >  1:1060638/1‑150 (MQ=255)
                                                                                    tcGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCgg                                                                                                      <  2:728421/111‑1 (MQ=255)
                                                                                    tcGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCgg                                                                                                      >  1:728421/1‑111 (MQ=255)
                                                                                        tgtgGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCaaa                                                           >  1:568039/1‑150 (MQ=255)
                                                                                        ggtgGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGCCCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCaaa                                                           <  1:400908/149‑1 (MQ=255)
                                                                                         gtgGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAg                                                          >  1:879349/1‑150 (MQ=255)
                                                                                             gAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTcccc                                                                                                                 >  2:710155/1‑91 (MQ=255)
                                                                                             gAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTcccc                                                                                                                 <  1:710155/91‑1 (MQ=255)
                                                                                             gAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGca                                                                               >  2:370314/1‑125 (MQ=255)
                                                                                             gAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGca                                                                               <  1:370314/125‑1 (MQ=255)
                                                                                             gAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAAcc                                                              >  1:979721/1‑142 (MQ=255)
                                                                                             gAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAAcc                                                              <  2:979721/142‑1 (MQ=255)
                                                                                             gAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCAtt                                                      >  2:1837926/1‑150 (MQ=255)
                                                                                               aTGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCAtttt                                                    <  2:695046/150‑1 (MQ=255)
                                                                                                 gggTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAg                                                                        <  2:1213688/128‑1 (MQ=255)
                                                                                                 gggTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAg                                                                        >  1:1213688/1‑128 (MQ=255)
                                                                                                  ggTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTg                                                 <  2:312514/150‑1 (MQ=255)
                                                                                                   gttgttTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGt                                                <  1:977854/150‑1 (MQ=255)
                                                                                                    ttgttTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGCCCGTGCGCGGTCCCCCGTAGGAGCGGGATTACCCGCGAAGCGGGCCGGACACGCAGAAGCATAACCAAAGCACTTTTTGtt                                               >  2:1610839/1‑150 (MQ=255)
                                                                                                    ttgttTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGtt                                               <  1:980792/150‑1 (MQ=255)
                                                                                                       tttCCCAGAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCtt                                            >  1:1184020/1‑150 (MQ=255)
                                                                                                       tttCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGaa                                                                                         <  2:1291813/105‑1 (MQ=255)
                                                                                                       tttCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGaa                                                                                         >  1:1291813/1‑105 (MQ=255)
                                                                                                       tttCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCtt                                            <  2:1060638/150‑1 (MQ=255)
                                                                                                        ttCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCttt                                           <  2:568039/150‑1 (MQ=255)
                                                                                                           ccAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAc                                        >  1:1454443/1‑150 (MQ=255)
                                                                                                            cAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAgg                                                                           <  1:1199156/114‑1 (MQ=255)
                                                                                                            cAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAgg                                                                           >  2:1199156/1‑114 (MQ=255)
                                                                                                            cAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAcc                                       >  2:1212283/1‑150 (MQ=255)
                                                                                                             aaaGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAccc                                      <  1:1785092/150‑1 (MQ=255)
                                                                                                                gCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAccc                                      >  2:59522/1‑147 (MQ=255)
                                                                                                                gCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCCGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAccc                                      <  1:59522/147‑1 (MQ=255)
                                                                                                                    gATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                    <  1:614417/97‑1 (MQ=255)
                                                                                                                    gATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                    >  2:614417/1‑97 (MQ=255)
                                                                                                                           ccTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCtt                                            >  2:517117/1‑130 (MQ=255)
                                                                                                                           ccTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCtt                                            <  1:517117/130‑1 (MQ=255)
                                                                                                                                gcCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTtat                   <  2:1454443/150‑1 (MQ=255)
                                                                                                                                gcCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTtat                   <  1:1465665/150‑1 (MQ=255)
                                                                                                                                 ccGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTtata                  <  2:1184020/150‑1 (MQ=255)
                                                                                                                                       tctTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCt            <  1:1358623/150‑1 (MQ=255)
                                                                                                                                         ttCGCGGGTAAACCCGCTCCTACGGGGGACCGCGCACGGTCCgcgcctgtaggagcggatttacccgcgaagaggccggcacaggcagaagcataaccaaagcattttttgttctttaacccgtataaacccacccgttatagtcgctgg       >  1:96559‑M1/1‑42 (MQ=255)
                                                                                                                                            gcgGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCaa       >  1:825718/1‑150 (MQ=255)
                                                                                                                                              gggTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTg                                                 >  2:1176369/1‑106 (MQ=255)
                                                                                                                                              gggTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTg                                                 <  1:1176369/106‑1 (MQ=255)
                                                                                                                                                gTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTg   >  1:1250712/1‑150 (MQ=255)
                                                                                                                                                 tAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGc  <  2:1114890/150‑1 (MQ=255)
                                                                                                                                                    aTCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTg                                                 <  1:1673677/100‑1 (MQ=255)
                                                                                                                                                    aTCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTg                                                 >  2:1673677/1‑100 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGC  >  NC_002947/4896295‑4896589

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