Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,241,973 Δ1 bp 7.7% intergenic (+32/+26) grxD → / ← bfr‑II monothiol glutaredoxin/bacterioferritin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,241,9730A. 7.7% 338.4 / 22.6 117intergenic (+32/+26)grxD/bfr‑IImonothiol glutaredoxin/bacterioferritin
Reads supporting (aligned to +/- strand):  ref base A (55/53);  new base . (4/5);  total (59/58)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

Only 100 of 122 total aligned reads displayed.

GTGGGTAGCCGGTGAGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCACGTCA  >  NC_002947/1241824‑1242122
                                                                                                                                                      |                                                                                                                                                     
gtgGGTAGCCGGTGAGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGAAAGCATAAAAAAACCCGC‑GATCa                                                                                                                                                       <  2:2926153/150‑1 (MQ=255)
gtgGGTAGCCGGTGAGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGCCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCa                                                                                                                                                       <  1:104314/150‑1 (MQ=255)
  gggTAGCCGGTGAGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGc                                                                                                                                                     <  1:1233767/150‑1 (MQ=255)
   ggTAGCCGGTGAGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCg                                                                                                                                                    >  1:1278857/1‑150 (MQ=255)
    gTAGCCGGTGAGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCgg                                                                                                                                                   >  1:2139698/1‑150 (MQ=255)
        ccGGTGAGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGttt                                                                                                                                               >  2:1057670/1‑150 (MQ=255)
          ggTGAGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGttttt                                                                                                                                             <  1:1719497/150‑1 (MQ=255)
            tGAGCTGGTTGGCGGCCGCGAAAACATGCTGGAAATGTTCGAGAATGGTGAGCTGCCAACCCTGATCGAGGACGCAGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGttttttt                                                                                                                                           <  2:1077759/150‑1 (MQ=255)
            tGAGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCTGATC‑GCGGGttttttt                                                                                                                                           <  2:250926/150‑1 (MQ=255)
             gAGCTGGTTGGCGGCAGCGACATGATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCAAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTa                                                                                                                                          <  2:1302220/150‑1 (MQ=255)
                  ggTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGG‑TTTTTTATGGGGc                                                                                                                                    <  2:2517474/150‑1 (MQ=255)
                        cggcAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATa                                                                                                                               <  1:2260637/150‑1 (MQ=255)
                        cggcAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATa                                                                                                                               <  2:1948780/150‑1 (MQ=255)
                         ggcAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAg                                                                                                                              <  2:1019759/150‑1 (MQ=255)
                            aGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTATTTTTATGGCGCGAATAGTTa                                                                                                                           >  2:1582828/1‑150 (MQ=255)
                               gACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTActc                                                                                                                        >  1:965587/1‑150 (MQ=255)
                                aCATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTActct                                                                                                                       >  2:1835813/1‑150 (MQ=255)
                                  atcatGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAACAGTTACTCTTc                                                                                                                     <  1:2058623/150‑1 (MQ=255)
                                      tGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCa                                                                                                                                                       <  2:2511537/112‑1 (MQ=255)
                                       gCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGGTTTTTTATGGGGCGAATAGTTACTCTTCGCCCa                                                                                                                >  2:976993/1‑150 (MQ=255)
                                            aaaTGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGc                                                                                                           <  2:2317564/150‑1 (MQ=255)
                                              aTGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGa                                                                                                         <  1:1057670/150‑1 (MQ=255)
                                                   cGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCa                                                                                                    <  2:833606/150‑1 (MQ=255)
                                                     agaAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCTGATC‑GCGGGttttttt                                                                                                                                           <  2:1280187/109‑1 (MQ=39)
                                                     agaAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCTGATC‑GCGGGttttttt                                                                                                                                           >  1:1280187/1‑109 (MQ=39)
                                                      gaAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGt                                                                                                 >  2:2369862/1‑150 (MQ=255)
                                                            tGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTtct                                                                                           <  1:1835813/150‑1 (MQ=255)
                                                            tGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTtct                                                                                           >  1:2793799/1‑150 (MQ=255)
                                                              ggCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCATGGCTTCTGGAGCCTGATCTGCCAGCAACCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGCTAGTTCTCg                                                                                         <  1:1582828/149‑1 (MQ=255)
                                                               gCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGACTTCAGAAGCATCATCATCAGCCCAACATAAACGAACAGGC‑AAAGAGCGTGTTTTTTGGTGTTGCTATTAGTTACTCTTCGCCCATCTGCTACTGCAGGTAGTTCTCGa                                                                                        >  2:2648493/1‑150 (MQ=255)
                                                               gCTGCAAACCCTCATCAAGGACGCGGCTCCCAAGGCTAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGAGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGa                                                                                        <  2:639730/150‑1 (MQ=255)
                                                                cTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGAt                                                                                       >  1:995898/1‑150 (MQ=255)
                                                                 tGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAACCCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATa                                                                                      <  1:1741905/150‑1 (MQ=255)
                                                                  gCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATAc                                                                                     <  2:1278857/150‑1 (MQ=255)
                                                                    aaaCCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCTGATC‑GCGGGtttttt                                                                                                                                            <  2:1709754/93‑1 (MQ=38)
                                                                    aaaCCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCTGATC‑GCGGGtttttt                                                                                                                                            >  1:1709754/1‑93 (MQ=38)
                                                                     aaCCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGa                                                                                  >  1:1156900/1‑150 (MQ=255)
                                                                     aaCCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGa                                                                                  >  1:602316/1‑150 (MQ=255)
                                                                      aCCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATTTGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGAt                                                                                 >  1:1573729/1‑150 (MQ=255)
                                                                        ccTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATc                                                                                <  2:1562670/149‑1 (MQ=255)
                                                                        ccTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCt                                                                               >  1:2564602/1‑150 (MQ=255)
                                                                        ccTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCt                                                                               >  2:2685048/1‑150 (MQ=255)
                                                                         cTGATCGAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCtt                                                                              >  2:2554784/1‑150 (MQ=255)
                                                                         cTGATCGAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCtt                                                                              >  2:2645729/1‑150 (MQ=255)
                                                                         cTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GAGCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCtt                                                                              >  1:843659/1‑150 (MQ=255)
                                                                           gATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGt                                                                            <  1:1113013/150‑1 (MQ=255)
                                                                           gATCAAGGACGCGACTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAGCCCGC‑GATCAGCGGGTTTGTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGt                                                                            <  2:2046746/150‑1 (MQ=255)
                                                                             tCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCTC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCg                                                                          >  2:1436427/1‑150 (MQ=255)
                                                                                      cgGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGc                                                                                                                   >  1:3029823/1‑100 (MQ=255)
                                                                                           gccaaggccaagGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGa                                                                                        >  2:2110790/1‑122 (MQ=255)
                                                                                           gccaaggccaagGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGa                                                                                        <  1:2110790/122‑1 (MQ=255)
                                                                                           gccaaggccaagGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTg                                                            <  2:602316/150‑1 (MQ=255)
                                                                                            ccaaggccaaggCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGccc                                                                                                                 <  2:281241/96‑1 (MQ=255)
                                                                                            ccaaggccaaggCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGccc                                                                                                                 >  1:281241/1‑96 (MQ=255)
                                                                                             caaggccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAgg                                                                                                  >  1:310796/1‑110 (MQ=255)
                                                                                             caaggccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAgg                                                                                                  <  2:310796/110‑1 (MQ=255)
                                                                                             caaggccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTggg                                                          >  1:2166304/1‑150 (MQ=255)
                                                                                             caaggccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGGTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTggg                                                          >  1:2794566/1‑150 (MQ=255)
                                                                                              aaggccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTCCTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTtggc                                                         <  2:2422981/150‑2 (MQ=255)
                                                                                              aaggccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGt                                                         <  1:2645729/150‑1 (MQ=255)
                                                                                              aaggccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGt                                                         <  1:2554784/150‑1 (MQ=255)
                                                                                                ggccaaggcTGCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGc                                                                                                           <  2:1949796/98‑1 (MQ=255)
                                                                                                 tccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGtt                                                             >  1:1319970/2‑143 (MQ=255)
                                                                                                 gccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCTC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTc                                                      <  1:1436427/150‑1 (MQ=255)
                                                                                                 gccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTc                                                      >  1:99782/1‑150 (MQ=255)
                                                                                                 gccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAACAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTc                                                      <  2:593158/150‑1 (MQ=255)
                                                                                                    aaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTccag                                                   >  1:1635499/1‑150 (MQ=255)
                                                                                                    aaggcTTCTGAAGCCTGATCTTCCGGCAAGAATAAAAAAACCCGC‑GATAAGCGGGGTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATAAGCCCCAGTTGGGTTTccag                                                   >  1:385044/1‑150 (MQ=255)
                                                                                                       gcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGcccc                                                                 >  2:74289/1‑133 (MQ=255)
                                                                                                       gcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGcccc                                                                 <  1:74289/133‑1 (MQ=255)
                                                                                                         ttCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGt                                              >  1:662766/1‑150 (MQ=255)
                                                                                                          tCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTc                                             <  2:965587/150‑1 (MQ=255)
                                                                                                            tGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGa                                           >  2:704849/1‑150 (MQ=255)
                                                                                                              aaGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCGTAATCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGAtg                                         >  1:2493407/1‑150 (MQ=255)
                                                                                                                gCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGAtgtt                                       >  2:3065302/1‑149 (MQ=255)
                                                                                                                 ccTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGc                                      >  1:73327/1‑150 (MQ=255)
                                                                                                                   tGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTggg                                                          <  1:2538627/128‑1 (MQ=255)
                                                                                                                      tCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGAt                                                                        >  2:342714/1‑111 (MQ=255)
                                                                                                                       cTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCGCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTcc                                >  1:1566939/1‑150 (MQ=255)
                                                                                                                        ttCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGc                                      >  2:1346102/1‑143 (MQ=255)
                                                                                                                        ttCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGc                                      <  1:1346102/143‑1 (MQ=255)
                                                                                                                        ttCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCt                               >  2:3071749/1‑150 (MQ=255)
                                                                                                                        ttCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCt                               >  1:2020617/1‑150 (MQ=255)
                                                                                                                         tCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTc                              <  1:2695398/150‑1 (MQ=255)
                                                                                                                                 gCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAgg                                                                                                  <  1:360164/74‑1 (MQ=255)
                                                                                                                                  cATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAg                     >  2:2911783/1‑150 (MQ=255)
                                                                                                                                  cATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAg                     <  2:2984167/150‑1 (MQ=255)
                                                                                                                                   aTAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAgg                    <  2:1243882/150‑1 (MQ=255)
                                                                                                                                       aaaaaCCCGCTGATC‑GCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAg                                                                                                   >  2:3017834/1‑67 (MQ=25)
                                                                                                                                       aaaaaCCCGCTGATC‑GCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAg                                                                                                   <  1:3017834/67‑1 (MQ=25)
                                                                                                                                       aaaaaCCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGAt                                          >  2:2131620/1‑124 (MQ=255)
                                                                                                                                        aaaaCCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGAt                                                                                 >  2:1956943/1‑84 (MQ=255)
                                                                                                                                        aaaaCCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGAt                                                                                 <  1:1956943/84‑1 (MQ=255)
                                                                                                                                         aaaCCCGCTGATC‑GCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTcc                                                     <  1:1671660/111‑1 (MQ=255)
                                                                                                                                           aCCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTc            <  1:2369862/150‑1 (MQ=255)
                                                                                                                                           aCCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTc            <  1:2911783/150‑1 (MQ=255)
                                                                                                                                             ccGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGCCTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGa          <  1:2807621/150‑1 (MQ=255)
                                                                                                                                               gc‑gATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGc        >  2:1165462/1‑150 (MQ=255)
                                                                                                                                                  gATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCAc      >  1:664753/1‑150 (MQ=255)
                                                                                                                                                      aGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCACGTCa  >  2:584564/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GTGGGTAGCCGGTGAGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCACGTCA  >  NC_002947/1241824‑1242122

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