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.6% 825.8 / ‑5.9 258intergenic (+60/‑395)b3407/feoAhypothetical protein/ferrous iron transport protein A
Reads supporting (aligned to +/- strand):  ref base . (0/1);  major base C (133/122);  minor base T (0/1);  total (133/125)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009133,537,4042.G99.6% 830.7 / ‑4.8 258intergenic (+60/‑395)b3407/feoAhypothetical protein/ferrous iron transport protein A
Reads supporting (aligned to +/- strand):  ref base . (0/1);  major base G (132/124);  minor base T (1/0);  total (133/125)
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.90e-01

Only 100 of 283 total aligned reads displayed.

GCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑CCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCAT  >  NC_000913/3537312‑3537486
                                                                                             ||                                                                                    
 cgcTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑c                                                                                     >  2:1781796/1‑93 (MQ=255)
 cgcTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑c                                                                                     >  2:3890975/1‑93 (MQ=255)
  gcTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑cg                                                                                    >  1:941493/1‑92 (MQ=255)
   cTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑cgc                                                                                   <  1:4195117/93‑3 (MQ=255)
     ggcggcggcAGTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑cgccg                                                                                 <  1:3900514/93‑5 (MQ=255)
     ggcggcggcAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑cgccg                                                                                 <  1:5840236/93‑5 (MQ=255)
     ggcggcggcAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑cgccg                                                                                 >  2:4620130/1‑89 (MQ=255)
     ggcggcggcAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑cgccg                                                                                 >  1:1900724/1‑89 (MQ=255)
           ggcAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTg                                                                           >  1:3990944/1‑93 (MQ=38)
              aaTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCa                                                                        >  1:2777413/1‑93 (MQ=37)
              aaTGGGAAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCa                                                                        <  1:5762673/93‑1 (MQ=37)
               aTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaa                                                                       >  2:4003725/1‑93 (MQ=37)
                tGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaac                                                                      >  1:2941682/1‑93 (MQ=37)
                  ggCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaacga                                                                    <  2:2880458/93‑1 (MQ=38)
                   gCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaacgaa                                                                   <  2:6594218/93‑1 (MQ=38)
                    cAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaacgaac                                                                  <  2:1107330/93‑1 (MQ=37)
                     aaaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaacgaacg                                                                 >  2:4791266/1‑93 (MQ=37)
                      aaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGc                                                                >  1:53248/1‑93 (MQ=37)
                      aaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGc                                                                >  1:5157604/1‑93 (MQ=37)
                      aaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGc                                                                >  1:3206863/1‑93 (MQ=37)
                      aaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGc                                                                >  1:338996/1‑93 (MQ=37)
                       aaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCCACGAACGcc                                                               >  1:6249906/1‑93 (MQ=37)
                       aaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGcc                                                               >  1:5705026/1‑93 (MQ=37)
                        aaaCGTTAACCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCt                                                              <  1:3170630/93‑1 (MQ=37)
                        aaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCATGCCGCATCTGGAAACGAACGCCt                                                              <  2:193431/93‑1 (MQ=38)
                        aaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCt                                                              <  2:785351/93‑1 (MQ=37)
                        aaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCt                                                              <  2:559653/93‑1 (MQ=37)
                        aaaCGTTAAACGCCCGTACTGGCCTACGGCTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCt                                                              <  2:4290483/93‑1 (MQ=37)
                         aaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCtt                                                             >  2:3544080/1‑93 (MQ=37)
                         aaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCtt                                                             <  2:6356491/93‑1 (MQ=37)
                          aCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTg                                                            >  1:2310889/1‑93 (MQ=37)
                           cGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGc                                                           >  2:3161679/1‑93 (MQ=37)
                           cGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGc                                                           >  2:5413114/1‑93 (MQ=37)
                            gTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGcc                                                          <  2:6262046/93‑1 (MQ=37)
                             ttAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCt                                                         >  2:6053175/1‑93 (MQ=38)
                               aagcGCCCGTACTGCCGTACGGTTCAATTTAGCACGAATTCGGAGGGCAGATAGGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGa                                                       <  2:2304466/90‑1 (MQ=16)
                                aaCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGAc                                                      >  2:4228074/1‑93 (MQ=37)
                                aaCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGAc                                                      >  2:503491/1‑93 (MQ=38)
                                aaCGCCCAGACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGAc                                                      >  1:6609044/1‑93 (MQ=37)
                                 aCGCCGTTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACa                                                     <  2:2799470/93‑1 (MQ=25)
                                 aCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACa                                                     <  1:4627918/93‑1 (MQ=37)
                                    cccGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTCACAtta                                                  >  1:489853/1‑91 (MQ=37)
                                     ccGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATaaa                                                 >  2:6520883/1‑93 (MQ=38)
                                       gTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAg                                               >  1:5993073/1‑93 (MQ=37)
                                       gTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAg                                               <  2:4476671/93‑1 (MQ=37)
                                        tACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGt                                              <  1:537905/93‑1 (MQ=37)
                                         aCTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTg                                             <  2:3655830/93‑1 (MQ=37)
                                         aCTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTg                                             >  2:2215938/1‑93 (MQ=37)
                                          cTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGc                                            >  1:6343286/1‑93 (MQ=38)
                                          cTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGc                                            >  2:554081/1‑93 (MQ=38)
                                          cTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGc                                            <  1:2434961/93‑1 (MQ=37)
                                           tGGCCTACGGTTCGAATTTGCACGAAATCGTAGGTCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGcc                                           <  1:2531837/93‑1 (MQ=37)
                                           tGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGcc                                           <  1:2080163/93‑1 (MQ=37)
                                           tGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGcc                                           >  2:1691501/1‑93 (MQ=38)
                                            ggCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCg                                          >  2:5322414/1‑93 (MQ=37)
                                                tACGGTTCTAATTTGCACGAAATCGTGGGGCAAATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGgaga                                      <  1:1552265/93‑1 (MQ=37)
                                                tACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGggga                                      >  2:1297606/1‑93 (MQ=37)
                                                tACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGgaga                                      >  2:2232581/1‑93 (MQ=38)
                                                tACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGgaga                                      >  2:2633005/1‑93 (MQ=37)
                                                tACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGgaga                                      >  2:4855727/1‑93 (MQ=38)
                                                tACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGgaga                                      <  1:1075096/93‑1 (MQ=37)
                                                tACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGgaga                                      >  1:1726411/1‑93 (MQ=38)
                                                 aCGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGaa                                     <  2:2992182/93‑1 (MQ=37)
                                                 aCGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGaa                                     <  2:310240/93‑1 (MQ=37)
                                                    gTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGGGAatat                                  >  2:32916/1‑93 (MQ=37)
                                                    gTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAatat                                  >  1:1674507/1‑93 (MQ=38)
                                                    gTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAatat                                  >  2:5213723/1‑93 (MQ=38)
                                                       cGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATAtctc                               >  1:1311455/1‑93 (MQ=38)
                                                        gAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTcc                              <  2:6681520/93‑1 (MQ=37)
                                                        gAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTcc                              >  1:5897069/1‑93 (MQ=38)
                                                        gAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑CGCGGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTcc                              <  1:6662093/93‑1 (MQ=37)
                                                            ttGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCa                          <  2:3776455/93‑1 (MQ=38)
                                                             tGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAt                         <  1:3014121/93‑1 (MQ=38)
                                                               cACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAttt                       <  1:4120386/93‑1 (MQ=38)
                                                               cACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCAGAGAATATCTCCGGCAttt                       <  1:6537442/93‑1 (MQ=37)
                                                                aCGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAtttt                      <  1:1770455/93‑1 (MQ=38)
                                                                 cGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCCGCCttttt                     >  1:215811/1‑93 (MQ=37)
                                                                   aaaTCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTAt                   >  1:1493865/1‑93 (MQ=38)
                                                                    aaTCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTAtt                  <  2:1867626/93‑1 (MQ=38)
                                                                    aaTCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTAtt                  >  1:860460/1‑93 (MQ=38)
                                                                    aaTCGTAGGGAAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTAtt                  <  1:82027/93‑1 (MQ=38)
                                                                     aTCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTc                 <  2:5105284/93‑1 (MQ=38)
                                                                      tCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTcc                >  1:1108771/1‑93 (MQ=38)
                                                                      tCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTcc                <  2:5877512/93‑1 (MQ=38)
                                                                      aggTAGGCCAGATACGGCGTTCACGCCGCATCTGGCAACGAACGCATTGCCTGACATAAAACTGCCGGAGAATATCTCCGGCATTTTTATTcc                <  1:3713527/91‑1 (MQ=25)
                                                                         tAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCCGAGAATATCTCCGGCATTTTTATTCcaca             >  2:5757712/1‑93 (MQ=38)
                                                                          aGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAg            >  2:1645011/1‑93 (MQ=38)
                                                                           gggCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTACGGCATTTTTATTCCACAGc           >  1:2681439/1‑93 (MQ=38)
                                                                            ggCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGcc          >  2:5342805/1‑93 (MQ=38)
                                                                             gCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCa         <  2:3379896/93‑1 (MQ=38)
                                                                               aGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCaaa       <  1:1838252/93‑1 (MQ=38)
                                                                               aGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCaaa       >  2:4427472/1‑93 (MQ=38)
                                                                                gATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGGATATCTCCCGCATTTTTATTCCACAGCCAAAc      >  2:4756386/1‑93 (MQ=37)
                                                                                gATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAAc      <  1:6464487/93‑1 (MQ=38)
                                                                                  tAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTc    <  1:6116415/93‑1 (MQ=38)
                                                                                  tAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTc    >  2:2028966/1‑93 (MQ=39)
                                                                                  tAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTc    >  1:3806579/1‑93 (MQ=39)
                                                                                   aaGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCa   >  2:5712255/1‑93 (MQ=39)
                                                                                   aaGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCa   <  2:1243892/93‑1 (MQ=39)
                                                                                    aGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCat  >  2:1573508/1‑93 (MQ=255)
                                                                                             ||                                                                                    
GCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑CCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCAT  >  NC_000913/3537312‑3537486

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: