Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,241,967 Δ1 bp 7.9% 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,9670A. 7.9% 241.3 / 13.8 76intergenic (+32/+26)grxD/bfr‑IImonothiol glutaredoxin/bacterioferritin
Reads supporting (aligned to +/- strand):  ref base A (33/37);  new base . (4/2);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01

GTAGCCGGTGAGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCAC  >  NC_002947/1241822‑1242112
                                                                                                                                                  |                                                                                                                                                 
gTAGCCGGTGAGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCg                                                                                                                                                >  7:35199/1‑149 (MQ=255)
  aGCCGGTGAGCTGGTTGGCGGCAGCTACTTAATGCTGTAAATGTTCGAGAAGGGTGAGCTTCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAATATACCCGCGAGCAGCgg                                                                                                                                               <  8:99957/148‑1 (MQ=255)
    ccGGTGAGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCTGATCGCGGGtt                                                                                                                                            <  8:131689/149‑1 (MQ=255)
          aGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGttttttt                                                                                                                                       <  8:268489/148‑1 (MQ=255)
          aGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTa                                                                                                                                      <  6:232612/149‑1 (MQ=255)
           gCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCTGATCGCGGGTTTTTTTAt                                                                                                                                     <  2:106933/149‑1 (MQ=255)
            cTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATg                                                                                                                                    <  8:207713/149‑1 (MQ=255)
               gctGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATgg                                                                                                                                   <  7:156586/145‑1 (MQ=255)
               gTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATgggg                                                                                                                                 <  3:198747/149‑1 (MQ=255)
                ttGGCTGCAGCGACATCATGCTGGAAATGTTCGAGAATGGTGAGCTTCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGc                                                                                                                                <  8:198738/149‑1 (MQ=255)
                ttGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGc                                                                                                                                >  5:173768/1‑149 (MQ=255)
                   gcggcAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGaa                                                                                                                             <  6:173768/149‑1 (MQ=255)
                          cGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTa                                                                                                                       >  7:221062/1‑148 (MQ=255)
                          cGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTAc                                                                                                                      <  7:183571/149‑1 (MQ=255)
                              atcatGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGACGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGGTTTTTTATGGGGCGAATAGTTACTCtt                                                                                                                  >  8:120678/1‑149 (MQ=255)
                                     tGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCAt                                                                                                           >  2:64502/1‑149 (MQ=255)
                                     tGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATc                                                                                                          >  4:206267/1‑149 (MQ=255)
                                         aaTGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGc                                                                                                       >  5:47074/1‑149 (MQ=255)
                                                gagaAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTg                                                                                                  >  4:284623/1‑147 (MQ=255)
                                                 agaAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAg                                                                                               <  1:64502/149‑1 (MQ=255)
                                                    aGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATATTCCGGCAAGCATAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAg                                                                                           <  3:206267/149‑1 (MQ=255)
                                                     gggTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAg                                                                                           >  4:258011/1‑149 (MQ=255)
                                                       gTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGtt                                                                                         >  5:128219/1‑149 (MQ=255)
                                                       gTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGt                                                                                          >  1:155995/1‑148 (MQ=255)
                                                         taGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTtct                                                                                       >  6:185919/2‑149 (MQ=255)
                                                         gAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTtct                                                                                       >  6:299248/1‑149 (MQ=255)
                                                            cTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGa                                                                                    <  5:299248/149‑1 (MQ=255)
                                                             tGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGa                                                                                    >  8:291921/1‑148 (MQ=255)
                                                               cAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATa                                                                                  <  6:257740/148‑1 (MQ=255)
                                                               cAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATa                                                                                  >  5:257740/1‑148 (MQ=255)
                                                               cAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCAATCTGCGACTGCAGGTAGTTCTCGATAc                                                                                 >  8:86717/1‑149 (MQ=255)
                                                                    ccTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATc                                                                            >  1:70616/1‑149 (MQ=255)
                                                                    ccTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATc                                                                            >  3:283394/1‑149 (MQ=255)
                                                                       gATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTg                                                                         >  4:269062/1‑149 (MQ=255)
                                                                         tCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTc                                                                       <  2:155995/149‑1 (MQ=255)
                                                                         tCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTc                                                                       <  6:128219/149‑1 (MQ=255)
                                                                           aaGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATc                                                                            <  4:93558/142‑1 (MQ=255)
                                                                           aaGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATc                                                                            >  3:93558/1‑142 (MQ=255)
                                                                           aaGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCg                                                                      <  5:18572/148‑1 (MQ=255)
                                                                           aaGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCg                                                                      <  6:47074/148‑1 (MQ=255)
                                                                             ggACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATc                                                                   <  8:221062/149‑1 (MQ=255)
                                                                                cgcgGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCATGTCGATCAGc                                                                >  6:149658/1‑149 (MQ=255)
                                                                                      tGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGt                                                          >  7:132834/1‑149 (MQ=255)
                                                                                        ccaaggccaaggCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTg                                                        <  8:35199/149‑1 (MQ=255)
                                                                                             gccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGAt                                                                             >  3:248782/1‑123 (MQ=255)
                                                                                             gccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGAt                                                                             <  4:248782/123‑1 (MQ=255)
                                                                                             gccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGttt                                                   >  4:301661/1‑149 (MQ=255)
                                                                                             gccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGttt                                                   >  7:80543/1‑149 (MQ=255)
                                                                                              ccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTc                                                  >  7:28242/1‑149 (MQ=255)
                                                                                                aaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTc                                                  >  3:112429/1‑147 (MQ=255)
                                                                                                   gcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTccagc                                              >  1:109980/1‑148 (MQ=255)
                                                                                                    cTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTccagcca                                            <  4:112429/149‑1 (MQ=255)
                                                                                                       cTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTATATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTGTCCAGCCAGTc                                         >  4:60529/1‑149 (MQ=255)
                                                                                                        tGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGt                                                                        <  1:162498/117‑1 (MQ=255)
                                                                                                        tGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGt                                                                        >  2:162498/1‑117 (MQ=255)
                                                                                                        tGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTc                                         >  8:222442/1‑148 (MQ=255)
                                                                                                          aaGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGAt                                      <  7:73659/149‑1 (MQ=255)
                                                                                                          aaGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCATGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCg                                        <  5:149658/147‑1 (MQ=255)
                                                                                                           aGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGAt                                      <  3:258011/148‑1 (MQ=255)
                                                                                                                      ccGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCt                           <  3:284623/148‑1 (MQ=255)
                                                                                                                          gccccATAAAAAAACCCGCTGATCGCGGGTGTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGAc                      >  7:273494/5‑149 (MQ=39)
                                                                                                                              cATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGa                  <  2:109980/149‑1 (MQ=255)
                                                                                                                                tAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAgg                <  4:283394/149‑1 (MQ=255)
                                                                                                                                tAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAgg                <  2:70616/149‑1 (MQ=255)
                                                                                                                                 aaaaaaaCCCGCTGATCCCGGGTTTTTTTATGGGGCGAATATTTACTCTTCGCCCATCAGAGACTGCAGGTAGTTCTCGATACCTATCTTGAAGATAATCCCCAGTTCGCTTACCAGCCACTCGATGTGCTCTTCCTCGGACTCGAGGa               >  4:212962/1‑149 (MQ=255)
                                                                                                                                  aaaaaaCCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGat              <  8:28242/149‑1 (MQ=255)
                                                                                                                                  aaaaaaCCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGa               <  3:269062/148‑1 (MQ=255)
                                                                                                                                    aaaaCCCGCTGATCGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGata             >  6:32753/1‑148 (MQ=255)
                                                                                                                                    aaaaCCCGCTGATCGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGata             >  4:78849/1‑148 (MQ=255)
                                                                                                                                         ccGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTc        >  3:129221/1‑148 (MQ=255)
                                                                                                                                         ccGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTc        >  3:49503/1‑148 (MQ=255)
                                                                                                                                         ccGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTc        <  4:129221/148‑1 (MQ=255)
                                                                                                                                           gcgATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTTTCGATCAGCCGCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAg     <  7:291921/149‑1 (MQ=255)
                                                                                                                                            cgATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGc    <  4:49503/149‑1 (MQ=255)
                                                                                                                                              gATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCa   <  5:185919/149‑1 (MQ=255)
                                                                                                                                               aTCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCAc  <  7:222442/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GTAGCCGGTGAGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCAC  >  NC_002947/1241822‑1242112

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