Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 901,421 G→A 5.4% intergenic (‑137/‑194) 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,4210GA 5.4% 586.8 / 14.6 186intergenic (‑137/‑194)cysZ/trxBsulfate transport protein CysZ/thioredoxin reductase
Reads supporting (aligned to +/- strand):  ref base G (85/91);  new base A (6/4);  total (91/95)
Fisher's exact test for biased strand distribution p-value = 5.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 189 total aligned reads displayed.

CAGGGGCTTGCATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAGCC  >  NC_002947/901272‑901569
                                                                                                                                                     |                                                                                                                                                    
cAGGGGCTTGCATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCg                                                                                                                                                      >  2:3716986/1‑150 (MQ=255)
  ggggCTTGCATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGAt                                                                                                                                                    <  2:3599096/150‑1 (MQ=255)
  ggggCTTGCATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGAt                                                                                                                                                    <  2:2710880/150‑1 (MQ=255)
   gggCTTGCATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATg                                                                                                                                                   <  2:4087903/150‑1 (MQ=255)
     gCTTGCATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGcc                                                                                                                                                 <  1:5284978/150‑1 (MQ=255)
       ttGCATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTa                                                                                                                                               >  2:1881860/1‑150 (MQ=255)
         gCATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACCAAAGGCACGCCGTGGGCGGGAAACGGTATTTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATg                                                                                                                                             <  2:4822370/150‑1 (MQ=255)
          cATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGTATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGa                                                                                                                                            <  1:5082051/150‑1 (MQ=255)
          cATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGa                                                                                                                                            <  2:2191570/150‑1 (MQ=255)
            tCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAgg                                                                                                                                          >  1:868751/1‑150 (MQ=255)
             cTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGt                                                                                                                                         >  1:3316974/1‑150 (MQ=255)
             cTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGt                                                                                                                                         <  2:1982291/150‑1 (MQ=255)
             cTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCGGGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGt                                                                                                                                         >  1:50580/1‑150 (MQ=255)
             cTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTGACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGt                                                                                                                                         <  2:2503160/150‑1 (MQ=255)
                gagaCTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGa                                                                                                                                      >  1:4380390/1‑150 (MQ=255)
                gagaCTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTGGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGa                                                                                                                                      >  1:1071895/1‑150 (MQ=255)
                  gaCTCCGCAGATAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGAtt                                                                                                                                    <  1:995106/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:4188056/150‑1 (MQ=255)
                        gCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTaa                                                                                                                              <  2:4498571/150‑1 (MQ=255)
                         cAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAg                                                                                                                             <  2:5245316/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)
                                 gAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTc                                                                                                                     >  2:1777502/1‑150 (MQ=255)
                                   aaaCGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCt                                                                                                                   >  2:2814555/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)
                                           tGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCttt                                                                                                           <  1:4835668/150‑1 (MQ=255)
                                                    cGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGAtt                                                                                                                                    <  2:3884426/116‑1 (MQ=255)
                                                           gCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGt                                                                                           >  2:2128722/1‑150 (MQ=255)
                                                              gcgcCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTgcg                                                                                        <  2:3668955/150‑1 (MQ=255)
                                                                   gacgaAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCaa                                                                                   >  1:3386222/1‑150 (MQ=255)
                                                                    acgaAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCttt                                                                                                           >  1:1784694/1‑125 (MQ=255)
                                                                    acgaAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCttt                                                                                                           <  2:1784694/125‑1 (MQ=255)
                                                                      gaAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACtt                                                                                >  1:4918812/1‑150 (MQ=255)
                                                                        aaGGCACGCCGTGGTCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGAGTGTCTAAGGATATTTCCTTAATATTTGCGACCTCGCTACAGTGCGCTCAACtttt                                                                              >  2:3826837/1‑150 (MQ=255)
                                                                              cGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCttt                                                                        >  1:1977418/1‑150 (MQ=255)
                                                                              cGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCttt                                                                        >  2:2083903/1‑150 (MQ=255)
                                                                                ccGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGc                                                                      <  1:272003/150‑1 (MQ=255)
                                                                                  gTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCgg                                                                    <  2:1231045/150‑1 (MQ=255)
                                                                                  atGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCgg                                                                    <  1:4480306/149‑1 (MQ=255)
                                                                                    gggCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGc                                                                  <  1:5000677/150‑1 (MQ=255)
                                                                                        gggAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTgc                                                              >  1:2886723/1‑150 (MQ=255)
                                                                                         ggAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTgcg                                                             >  1:4186002/1‑150 (MQ=255)
                                                                                          gAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTCCAGTGCGCTCAACTTTTCGCTTTGCGGGCCTgcgc                                                            >  2:2934656/1‑150 (MQ=255)
                                                                                          gAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTgcgc                                                            <  1:5441496/150‑1 (MQ=255)
                                                                                          gAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTgcgc                                                            <  1:51143/150‑1 (MQ=255)
                                                                                              cGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGttt                                                        >  1:4930778/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)
                                                                                                 gATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAg                                                     <  2:4918812/150‑1 (MQ=255)
                                                                                                   tGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGGGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGcc                                                   >  2:5331470/1‑150 (MQ=255)
                                                                                                    gTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCt                                                  >  1:533196/1‑150 (MQ=255)
                                                                                                      aaCAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTc                                                >  1:1520758/1‑150 (MQ=255)
                                                                                                      aaCAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTc                                                <  2:521598/150‑1 (MQ=255)
                                                                                                      aaCAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTc                                                >  1:2224695/1‑150 (MQ=255)
                                                                                                         aaaaGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTcctg                                             >  1:1366675/1‑148 (MQ=255)
                                                                                                         aaaaGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTcctg                                             >  1:2615588/1‑148 (MQ=255)
                                                                                                           aaGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCaa                                           >  1:3167911/1‑150 (MQ=255)
                                                                                                           aaGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCaa                                           >  1:2103266/1‑150 (MQ=255)
                                                                                                            aGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAg                                          >  1:3046083/1‑150 (MQ=255)
                                                                                                             gcgcCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGt                                         <  1:1887062/150‑1 (MQ=255)
                                                                                                                    gggCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGt                                  >  2:4072595/1‑150 (MQ=255)
                                                                                                                       cTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTggg                               <  1:3716986/150‑1 (MQ=255)
                                                                                                                         tGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTc                             <  2:1603051/150‑1 (MQ=255)
                                                                                                                           tttATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGttt                                                        <  2:3507193/121‑1 (MQ=255)
                                                                                                                           tttATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGttt                                                        >  1:3507193/1‑121 (MQ=255)
                                                                                                                            ttATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCtt                          <  1:408815/150‑1 (MQ=255)
                                                                                                                               tCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTctctaga                       <  2:868751/150‑7 (MQ=255)
                                                                                                                                 gATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTaa                     >  1:4619583/1‑150 (MQ=255)
                                                                                                                                  aTGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAAt                    >  1:2473502/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)
                                                                                                                                        tCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTtg              >  2:1153045/1‑150 (MQ=255)
                                                                                                                                        tCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTtg              <  1:2385369/150‑1 (MQ=255)
                                                                                                                                        tCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTtg              <  2:3925066/150‑1 (MQ=255)
                                                                                                                                        tCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTtg              >  2:3417674/1‑150 (MQ=255)
                                                                                                                                         cGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTtgg             <  1:626731/150‑1 (MQ=255)
                                                                                                                                          gcagcaTAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTtggc            >  2:1484/1‑150 (MQ=255)
                                                                                                                                           cagcaTAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTtggct           <  2:1071895/150‑1 (MQ=255)
                                                                                                                                           cagcaTAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTtggct           <  2:50580/150‑1 (MQ=255)
                                                                                                                                            ctcaTAGGCATCGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTa                                                      >  2:3172862/3‑106 (MQ=255)
                                                                                                                                            ctcaTAGGCATCGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTa                                                      <  1:3172862/104‑1 (MQ=255)
                                                                                                                                            ctcaTAGGCATCGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCtt                          >  1:2757089/3‑134 (MQ=255)
                                                                                                                                            ctcaTAGGCATCGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTtg              <  1:3836683/144‑1 (MQ=255)
                                                                                                                                                   gCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAGc   >  1:4294568/1‑150 (MQ=255)
                                                                                                                                                    cGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAGcc  <  2:3316979/150‑1 (MQ=255)
                                                                                                                                                    cGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAGcc  <  2:3316974/150‑1 (MQ=255)
                                                                                                                                                    cGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAGcc  <  2:3101886/150‑1 (MQ=255)
                                                                                                                                                    cGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAGcc  <  2:3865885/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAGGGGCTTGCATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAAGGCACGCCGTGGGCGGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAGCC  >  NC_002947/901272‑901569

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