Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 901,422 A→T 5.6% intergenic (‑138/‑193) cysZ ← / → trxB sulfate transport protein CysZ/thioredoxin reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947901,4220AT 5.6% 629.8 / 17.1 188intergenic (‑138/‑193)cysZ/trxBsulfate transport protein CysZ/thioredoxin reductase
Reads supporting (aligned to +/- strand):  ref base A (86/91);  new base T (6/5);  total (92/96)
Fisher's exact test for biased strand distribution p-value = 7.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 188 total aligned reads displayed.

GGGGCTTGCATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAGCC  >  NC_002947/901274‑901569
                                                                                                                                                    |                                                                                                                                                   
ggggCTTGCATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTCCCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAACGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGAt                                                                                                                                                    <  2:3599149/150‑1 (MQ=255)
 gggCTTGCATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATg                                                                                                                                                   <  2:4087903/150‑1 (MQ=255)
   gCTTGCATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGcc                                                                                                                                                 <  1:5284978/150‑1 (MQ=255)
     ttGCATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTa                                                                                                                                               <  2:2024055/150‑1 (MQ=255)
       gCATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACCAAAGGCACGCCGTGGGCGGGAAACGGTATTTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATg                                                                                                                                             <  2:4822370/150‑1 (MQ=255)
        cATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGa                                                                                                                                            <  2:39030/150‑1 (MQ=255)
        cATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGa                                                                                                                                            <  2:2191570/150‑1 (MQ=255)
        aaTCTCGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTGGCCGCGCCGACTAACGGCACGCCGTGGGCGGGAAACGGGATGTAACAAACGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGa                                                                                                                                            <  1:5030034/149‑1 (MQ=255)
          tCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAgg                                                                                                                                          >  2:4042141/1‑150 (MQ=255)
           cTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGt                                                                                                                                         >  1:3316979/1‑150 (MQ=255)
           cTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGt                                                                                                                                         >  1:3316974/1‑150 (MQ=255)
           cTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTGACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGt                                                                                                                                         <  2:2503160/150‑1 (MQ=255)
            tAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTg                                                                                                                                        >  1:4447594/1‑150 (MQ=255)
              gagaCTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTGGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGa                                                                                                                                      >  1:1071895/1‑150 (MQ=255)
              gagaCACCGCAGAGAAGGAAAACGCGCGGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGa                                                                                                                                      <  1:3409204/150‑1 (MQ=255)
                gaCTCCGCAGAGAAGGAAAACGCGCGGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGAtt                                                                                                                                    <  1:800556/150‑1 (MQ=255)
                gaCGCCGCAGAGAAGGAAAACGCGCGGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGAtt                                                                                                                                    <  1:5069066/150‑1 (MQ=255)
                  cTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCGGCATAGGCGATGCCTATGAGGTGGATTGt                                                                                                                                  <  2:3435065/150‑1 (MQ=255)
                    ccGCAGAGAAGGAAAACGCGCGGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCt                                                                                                                                <  1:3229896/150‑1 (MQ=255)
                      gCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCTTGGTCGGGAAACGGGATGTAAAAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTaa                                                                                                                              <  2:4431599/150‑1 (MQ=255)
                      gCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTaa                                                                                                                              >  1:4503490/1‑150 (MQ=255)
                      gCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTaa                                                                                                                              <  2:4498571/150‑1 (MQ=255)
                      gCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTaa                                                                                                                              <  2:4188056/150‑1 (MQ=255)
                       cAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAg                                                                                                                             <  2:5245316/150‑1 (MQ=255)
                          agaAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGat                                                                                                                          <  1:4936341/150‑1 (MQ=255)
                           gaAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGata                                                                                                                         >  1:2322775/1‑150 (MQ=255)
                            aaGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGatat                                                                                                                        >  1:3057199/1‑150 (MQ=255)
                             aGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCGTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATAtt                                                                                                                       >  1:3925066/1‑150 (MQ=255)
                              ggAAATCGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATAttt                                                                                                                      >  2:4997515/1‑150 (MQ=255)
                              ggAAATCGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATAttt                                                                                                                      >  2:1795953/1‑150 (MQ=255)
                               gAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTc                                                                                                                     >  2:1777857/1‑150 (MQ=255)
                               gAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTc                                                                                                                     >  2:1777502/1‑150 (MQ=255)
                                 aaaCGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCt                                                                                                                   >  2:626731/1‑150 (MQ=255)
                                    cgcgCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTaa                                                                                                                <  1:2227399/150‑1 (MQ=255)
                                      cgcTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATc                                                                                                              <  1:3304513/150‑1 (MQ=255)
                                         tGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTACGGATATTTCCTTAATCttt                                                                                                           <  2:3480789/150‑1 (MQ=255)
                                         tGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCttt                                                                                                           <  1:4835668/150‑1 (MQ=255)
                                                  cGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGGTTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCg                                                                                                  >  2:2922160/1‑150 (MQ=255)
                                                  cGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCg                                                                                                  >  2:3085609/1‑150 (MQ=255)
                                                       ttGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAACGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACa                                                                                             <  2:3569145/150‑1 (MQ=255)
                                                         gCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGt                                                                                           >  1:285600/1‑150 (MQ=255)
                                                                  acgaAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCttt                                                                                                           >  1:1784694/1‑125 (MQ=255)
                                                                   cgaAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACt                                                                                 >  1:3030903/1‑150 (MQ=255)
                                                                    gaAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACtt                                                                                >  1:4918812/1‑150 (MQ=255)
                                                                      aaGGCACGCCGTGGTCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGAGTGTCTAAGGATATTTCCTTAATATTTGCGACCTCGCTACAGTGCGCTCAACtttt                                                                              >  2:3826837/1‑150 (MQ=255)
                                                                            cGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCttt                                                                        >  2:2084193/1‑150 (MQ=255)
                                                                            cGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCttt                                                                        >  1:1977418/1‑150 (MQ=255)
                                                                              ccGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGc                                                                      <  1:272003/150‑1 (MQ=255)
                                                                                gTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCgg                                                                    <  1:4198612/150‑1 (MQ=255)
                                                                                gTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCgg                                                                    <  2:1231045/150‑1 (MQ=255)
                                                                                      gggAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTgc                                                              <  1:90345/150‑1 (MQ=255)
                                                                                       ggAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTgcg                                                             >  1:4186002/1‑150 (MQ=255)
                                                                                        gAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTgcgc                                                            <  1:5441496/150‑1 (MQ=255)
                                                                                        gAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTgcgc                                                            <  1:51143/150‑1 (MQ=255)
                                                                                            cGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGttt                                                        >  1:4930775/1‑150 (MQ=255)
                                                                                              ggATGTAACAACAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTa                                                      <  2:57024/150‑1 (MQ=255)
                                                                                              ggATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTa                                                      <  2:4744939/150‑1 (MQ=255)
                                                                                              ggATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTa                                                      <  2:2322775/150‑1 (MQ=255)
                                                                                                 tGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGGGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGcc                                                   >  2:5331470/1‑150 (MQ=255)
                                                                                                  gTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCt                                                  >  1:533196/1‑150 (MQ=255)
                                                                                                    aaCAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTc                                                <  2:521598/150‑1 (MQ=255)
                                                                                                    aaCAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTc                                                >  1:1520758/1‑150 (MQ=255)
                                                                                                       aaaaGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTcctg                                             >  1:1366675/1‑148 (MQ=255)
                                                                                                       aaaaGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTcctg                                             >  1:2615588/1‑148 (MQ=255)
                                                                                                         aaGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCaa                                           >  1:2103266/1‑150 (MQ=255)
                                                                                                         aaGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCaa                                           >  1:2103285/1‑150 (MQ=255)
                                                                                                         aaGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCaa                                           >  1:3167911/1‑150 (MQ=255)
                                                                                                           gcgcCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGt                                         <  1:1887062/150‑1 (MQ=255)
                                                                                                                  gggCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGt                                  >  2:4072595/1‑150 (MQ=255)
                                                                                                                    gCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTgg                                <  2:2558913/150‑1 (MQ=255)
                                                                                                                     cTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTggg                               <  1:3716986/150‑1 (MQ=255)
                                                                                                                       tGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTc                             <  2:1603051/150‑1 (MQ=255)
                                                                                                                          ttATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCtt                          <  1:408815/150‑1 (MQ=255)
                                                                                                                            aTCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCgg                                                                    <  2:3315239/106‑1 (MQ=255)
                                                                                                                             tCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTctctaga                       <  2:868751/150‑7 (MQ=255)
                                                                                                                               gATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTaa                     >  1:4619583/1‑150 (MQ=255)
                                                                                                                               gATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTaa                     >  2:1022753/1‑150 (MQ=255)
                                                                                                                                aTGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAAt                    >  1:2473502/1‑150 (MQ=255)
                                                                                                                                aTGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAAt                    >  2:4972702/1‑150 (MQ=255)
                                                                                                                                aTGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAAt                    >  1:192344/1‑150 (MQ=255)
                                                                                                                                 tGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAAtt                   >  1:3605813/1‑150 (MQ=255)
                                                                                                                                  gAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTa                                                      <  1:4511437/114‑1 (MQ=255)
                                                                                                                                  gAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAAtat                  >  1:3229246/1‑148 (MQ=255)
                                                                                                                                     aTCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCtt               <  1:1096391/150‑1 (MQ=255)
                                                                                                                                      tCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTATTTTCTtg              >  2:1152997/1‑150 (MQ=255)
                                                                                                                                      tCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTtg              >  2:3417674/1‑150 (MQ=255)
                                                                                                                                      tCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTtg              <  1:2385369/150‑1 (MQ=255)
                                                                                                                                       cGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTtgg             <  1:626731/150‑1 (MQ=255)
                                                                                                                                        gcagcaTAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTtggc            >  2:1484/1‑150 (MQ=255)
                                                                                                                                         cctcaTAGGCATCGCCTATGCGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTtggct           >  1:4331252/4‑150 (MQ=255)
                                                                                                                                         cagcaTAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTtggct           <  2:50580/150‑1 (MQ=255)
                                                                                                                                          ctcaTAGGCATCGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTa                                                      >  2:3172801/3‑106 (MQ=255)
                                                                                                                                          ctcaTAGGCATCGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTa                                                      <  1:3172862/104‑1 (MQ=255)
                                                                                                                                          ctcaTAGGCATCGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCtt                          <  2:2757089/132‑1 (MQ=255)
                                                                                                                                          ctcaTAGGCATCGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCtt                          >  1:2757089/3‑134 (MQ=255)
                                                                                                                                          ctcaTAGGCATCGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTtg              <  1:3836683/144‑1 (MQ=255)
                                                                                                                                                  cGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAGcc  <  2:3865885/150‑1 (MQ=255)
                                                                                                                                                  cGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAGcc  <  2:3707218/150‑1 (MQ=255)
                                                                                                                                                  cGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAGcc  <  2:3316979/150‑1 (MQ=255)
                                                                                                                                                   gATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTaa                     >  2:1908576/1‑130 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGGGCTTGCATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAGCC  >  NC_002947/901274‑901569

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