Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,537,405 +GC intergenic (+61/‑394) b3407 → / → feoA hypothetical protein/ferrous iron transport protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,537,4041.C99.5% 708.4 / ‑3.8 207intergenic (+60/‑395)b3407/feoAhypothetical protein/ferrous iron transport protein A
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (94/112);  total (94/113)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01
*NC_0009133,537,4042.G99.5% 701.8 / ‑3.8 206intergenic (+60/‑395)b3407/feoAhypothetical protein/ferrous iron transport protein A
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (93/112);  total (93/113)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 228 total aligned reads displayed.

GCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CA‑‑CCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATA  >  NC_000913/3537312‑3537487
                                                                                               ||                                                                                   
 cgcTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CA‑‑c                                                                                    >  1:3109640/1‑93 (MQ=255)
 cgcTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CA‑‑c                                                                                    <  2:4763745/93‑1 (MQ=255)
  gcTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CA‑‑cg                                                                                   <  1:3046300/93‑2 (MQ=255)
  gcTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CA‑‑cg                                                                                   >  1:1458736/1‑92 (MQ=255)
    tGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CA‑‑cgcc                                                                                 >  1:124998/1‑90 (MQ=255)
     ggcggcggcAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CA‑‑cgccg                                                                                >  2:1929073/1‑89 (MQ=255)
     ggcggcggcAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CA‑‑cgccg                                                                                <  1:2095403/93‑5 (MQ=255)
           ggcAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTg                                                                            >  2:3746183/1‑93 (MQ=38)
             cAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGc                                                                          <  2:2891767/93‑1 (MQ=38)
              aaTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCa                                                                         >  2:1718745/1‑93 (MQ=37)
               aTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCaa                                                                        >  2:1029538/1‑93 (MQ=37)
               aTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCaa                                                                        <  2:5704027/93‑1 (MQ=37)
                tGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCaac                                                                       >  1:5470616/1‑93 (MQ=37)
                tGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCaac                                                                       >  2:2093979/1‑93 (MQ=37)
                   gCAAAAAACGTTAAACGCCCGTCCTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCaacgaa                                                                    <  2:1411044/93‑1 (MQ=37)
                   gCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAACTCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCaacgaa                                                                    <  2:3430750/93‑1 (MQ=37)
                   gCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCaacgaa                                                                    <  1:6204638/93‑1 (MQ=37)
                    cAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCaacgaac                                                                   <  2:4244680/93‑1 (MQ=38)
                    aaaaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCaacgaac                                                                   <  2:5961112/92‑1 (MQ=37)
                     aaaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCaacgaacg                                                                  <  1:2640340/93‑1 (MQ=37)
                      aaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGc                                                                 >  1:5175749/1‑93 (MQ=37)
                      aaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGc                                                                 >  1:5907412/1‑93 (MQ=37)
                       aaaaCGTTAAACGCCCGTCCTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGcc                                                                <  1:5114533/93‑1 (MQ=37)
                       aaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGcc                                                                <  2:6150385/93‑1 (MQ=37)
                        aaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCt                                                               <  1:4652413/93‑1 (MQ=37)
                         aaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCtt                                                              >  1:6027643/1‑93 (MQ=37)
                         aaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCtt                                                              <  1:4760232/93‑1 (MQ=37)
                             ttAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCt                                                          <  2:44494/93‑1 (MQ=37)
                             ttAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCt                                                          <  1:1710516/93‑1 (MQ=37)
                                 aCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACa                                                      >  2:2460032/1‑93 (MQ=38)
                                 aCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACa                                                      <  2:1566221/93‑1 (MQ=37)
                                  cGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACAt                                                     <  1:1104389/93‑1 (MQ=37)
                                    cccGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATaa                                                   >  1:4371759/1‑93 (MQ=38)
                                    cccGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATaa                                                   <  2:5928501/93‑1 (MQ=37)
                                     ccGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATaaa                                                  >  1:5193074/1‑93 (MQ=38)
                                      cGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATaaaa                                                 >  1:1016234/1‑93 (MQ=38)
                                       gTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAg                                                <  2:4794862/93‑1 (MQ=37)
                                       gTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAg                                                <  2:3036817/93‑1 (MQ=37)
                                        tACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGt                                               <  1:4673244/93‑1 (MQ=37)
                                          cTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGc                                             >  1:3891962/1‑93 (MQ=38)
                                           tGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGcc                                            <  2:657145/93‑1 (MQ=37)
                                           tGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGcc                                            <  1:1004822/93‑1 (MQ=37)
                                           tGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGcc                                            <  2:4290371/93‑1 (MQ=37)
                                            ggCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCg                                           <  2:5290184/93‑1 (MQ=37)
                                            ggCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCg                                           <  2:1764763/93‑1 (MQ=37)
                                             gCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCgg                                          <  1:5034301/93‑1 (MQ=37)
                                              ccTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGga                                         >  1:1938630/1‑93 (MQ=38)
                                              ccTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGga                                         <  1:2848171/93‑1 (MQ=37)
                                               cTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCACATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGgag                                        <  2:4532772/93‑1 (MQ=37)
                                                 aCGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGaa                                      >  1:934712/1‑93 (MQ=38)
                                                 aCGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGaa                                      <  1:2487332/93‑1 (MQ=37)
                                                  cGGTTCGAATTTGCACGAAATCGTAGGGCAAATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAat                                     <  1:4782624/93‑1 (MQ=37)
                                                   ggTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAata                                    >  2:3876827/1‑93 (MQ=38)
                                                   ggTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTC‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAata                                    <  1:5737597/93‑1 (MQ=38)
                                                    gTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAatat                                   >  1:3716714/1‑93 (MQ=38)
                                                     ttCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATAtc                                  >  1:3619791/1‑93 (MQ=37)
                                                     ttCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATAtc                                  >  2:5849905/1‑93 (MQ=38)
                                                     ttCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATAtc                                  <  1:1486962/93‑1 (MQ=37)
                                                      tCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATAtct                                 >  1:1376656/1‑93 (MQ=38)
                                                       cGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATAtctc                                <  1:6374571/93‑1 (MQ=37)
                                                       cGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATAtctc                                <  2:3247239/93‑1 (MQ=37)
                                                        gAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCGCG‑‑CCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTcc                               <  1:4019082/93‑1 (MQ=38)
                                                        gAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTcc                               <  2:5822606/93‑1 (MQ=37)
                                                        gAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTcc                               >  2:2776963/1‑93 (MQ=38)
                                                         aaTTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCg                              >  2:3823460/1‑93 (MQ=38)
                                                         aaTTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCg                              >  2:4622797/1‑93 (MQ=38)
                                                         aaTTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCg                              >  2:5516652/1‑93 (MQ=38)
                                                         aaTTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCg                              <  1:5415579/93‑1 (MQ=37)
                                                          aTTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCgg                             <  2:1797710/93‑1 (MQ=38)
                                                          aTTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCgg                             >  1:2242395/1‑93 (MQ=38)
                                                          aTTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCgg                             >  1:4582995/1‑93 (MQ=38)
                                                              gCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAtt                         <  1:3455593/93‑1 (MQ=38)
                                                              gCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAtt                         <  2:2156495/93‑1 (MQ=38)
                                                              gCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAtt                         <  2:5201516/93‑1 (MQ=38)
                                                              gCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAtt                         >  1:1734846/1‑93 (MQ=38)
                                                              gCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAtt                         <  1:885284/93‑1 (MQ=38)
                                                              gCACGAAATCGTAGGCCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAtt                         <  1:2959127/93‑1 (MQ=37)
                                                               cACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAttt                        <  1:1019429/93‑1 (MQ=38)
                                                                aCGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAtttt                       <  2:5632316/93‑1 (MQ=38)
                                                                aCGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAtttt                       <  2:3986844/93‑1 (MQ=38)
                                                                 cGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAttttt                      <  2:4104115/93‑1 (MQ=38)
                                                                   aaaTCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTAt                    <  2:3951970/93‑1 (MQ=38)
                                                                   aaaTCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTAt                    >  1:5594604/1‑93 (MQ=38)
                                                                    aaTCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCACTtttttt                   >  1:38492/1‑93 (MQ=37)
                                                                      tCGTAGGGCAGATAAGGCGTC‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTcc                 <  2:10853/93‑1 (MQ=38)
                                                                        gTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCcac               <  1:4270920/93‑1 (MQ=38)
                                                                         tAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCcaca              <  2:2856225/93‑1 (MQ=38)
                                                                          aGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGTCGGAGTATATCTCCGGCATTTTTATTCCACAg             >  2:327575/1‑93 (MQ=37)
                                                                          aGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAg             >  2:505397/1‑93 (MQ=38)
                                                                           gggCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGc            >  2:3623265/1‑93 (MQ=38)
                                                                            ggCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGcc           <  2:5597206/93‑1 (MQ=38)
                                                                             gCAGATAAGGCTTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCa          <  2:5398467/93‑1 (MQ=38)
                                                                             gCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGGATATCTCCGGCATTTTTATTCCACAGCCa          >  2:1472803/1‑93 (MQ=38)
                                                                             gCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCa          <  1:2165695/93‑1 (MQ=38)
                                                                               aGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCaaa        >  2:5125104/1‑93 (MQ=38)
                                                                                gATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAAc       <  1:3266140/93‑1 (MQ=38)
                                                                                 aTAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACt      <  1:563473/93‑1 (MQ=38)
                                                                                    aGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCat   <  1:1660735/93‑1 (MQ=39)
                                                                                    aGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCat   >  1:5786614/1‑93 (MQ=255)
                                                                                     ggCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCata  <  1:778903/93‑1 (MQ=39)
                                                                                               ||                                                                                   
GCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CA‑‑CCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATA  >  NC_000913/3537312‑3537487

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