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.C100.0% 615.0 / NA 192intergenic (+60/‑395)b3407/feoAhypothetical protein/ferrous iron transport protein A
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (118/74);  total (118/74)
*NC_0009133,537,4042.G99.5% 617.0 / ‑5.2 192intergenic (+60/‑395)b3407/feoAhypothetical protein/ferrous iron transport protein A
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base G (117/74);  minor base T (1/0);  total (118/74)
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.93e-01

Only 100 of 205 total aligned reads displayed.

GCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑C‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATA  >  NC_000913/3537312‑3537487
                                                                                             ||                                                                                     
gcgcTGGCGGCGGCAATGGGCAAAAACCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCa                                                                                         <  1:275055/93‑1 (MQ=255)
 cgcTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑c                                                                                      <  2:1620019/93‑1 (MQ=255)
   cTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑cgc                                                                                    <  2:1508499/93‑3 (MQ=255)
    tGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑cgcc                                                                                   >  1:2836255/1‑90 (MQ=255)
           ggcAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAAGCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTg                                                                            <  1:814032/93‑1 (MQ=38)
            gcAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTgg                                                                           >  2:4533501/1‑93 (MQ=38)
             cAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTATGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGc                                                                          >  1:3256006/1‑93 (MQ=37)
             cAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGc                                                                          >  2:485078/1‑93 (MQ=37)
             cAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGc                                                                          <  2:2102796/93‑1 (MQ=37)
              aaTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCa                                                                         >  1:3621792/1‑93 (MQ=37)
               aTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCaa                                                                        >  2:3823528/1‑93 (MQ=37)
                tGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCaac                                                                       <  2:2033673/93‑1 (MQ=38)
                   gCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCaacgaa                                                                    <  1:432949/93‑1 (MQ=38)
                    cAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCaacgaac                                                                   >  1:2816214/1‑93 (MQ=37)
                     aaaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCaacgaacg                                                                  >  1:728355/1‑93 (MQ=37)
                     aaaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCaacgaacg                                                                  >  2:3594029/1‑93 (MQ=37)
                     aaaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCaacgaacg                                                                  <  1:2641814/93‑1 (MQ=37)
                      aaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGc                                                                 <  1:1837367/93‑1 (MQ=37)
                      aaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGc                                                                 <  1:3181852/93‑1 (MQ=37)
                        aaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCt                                                               >  1:969731/1‑93 (MQ=37)
                        aaaCGTGAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCt                                                               <  1:2036669/93‑1 (MQ=37)
                         accGTTAAACGCCCGTATTGGCCTACGTTTGGAATTTGCACAAAACCGTAGGGCAGATAAGGCGTCCACGC‑‑CGCATCTGGCAACGAACGCCtt                                                              <  2:1745127/91‑1 (MQ=21)
                         aaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCtt                                                              <  1:176704/93‑1 (MQ=37)
                         aaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCtt                                                              >  2:2048869/1‑93 (MQ=37)
                          aCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTg                                                             >  2:66866/1‑93 (MQ=37)
                          aCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTg                                                             <  2:721798/93‑1 (MQ=37)
                          aCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTg                                                             >  1:4042069/1‑93 (MQ=37)
                          aCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTg                                                             >  1:3570934/1‑93 (MQ=37)
                           cGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGc                                                            <  1:3800511/93‑1 (MQ=37)
                              tAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTg                                                         >  2:1326670/1‑93 (MQ=37)
                              tAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTg                                                         >  2:1940313/1‑93 (MQ=38)
                                aaCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGAc                                                       <  1:894293/93‑1 (MQ=37)
                                 acccccGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTGCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACa                                                      <  2:3616818/90‑1 (MQ=37)
                                 aCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACa                                                      >  2:2211154/1‑93 (MQ=37)
                                 aCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAAAGCCTTGCCCGaaa                                                      >  2:14065/1‑91 (MQ=25)
                                  cGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACAt                                                     <  1:271224/93‑1 (MQ=37)
                                  cGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACAt                                                     >  1:1688435/1‑93 (MQ=38)
                                   gCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATa                                                    >  1:200351/1‑93 (MQ=38)
                                   ccccGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATa                                                    <  2:4259790/92‑1 (MQ=37)
                                     ccGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATaaa                                                  <  1:4405524/93‑1 (MQ=37)
                                      cGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATaaaa                                                 >  1:2105470/1‑93 (MQ=38)
                                      cGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATaaaa                                                 <  2:1175278/93‑1 (MQ=37)
                                      cGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATaaaa                                                 <  1:1949404/93‑1 (MQ=37)
                                      cGTACGGGCATCCGGTTCGAATTTGCACGAAATGGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATaaaa                                                 <  2:1351123/93‑1 (MQ=25)
                                       gTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAg                                                >  2:2962591/1‑93 (MQ=38)
                                       gTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAg                                                >  2:3220841/1‑93 (MQ=38)
                                        tACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGt                                               >  1:2239009/1‑93 (MQ=37)
                                        tACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGt                                               <  2:3133242/93‑1 (MQ=37)
                                         aCTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTg                                              >  2:4503574/1‑93 (MQ=38)
                                         aCTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTg                                              <  1:4348377/93‑1 (MQ=37)
                                         aCTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTg                                              >  2:3971548/1‑93 (MQ=37)
                                           tGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGcc                                            >  1:4527218/1‑93 (MQ=38)
                                             gCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCgg                                          >  2:2874228/1‑93 (MQ=38)
                                              ccTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCggg                                         >  2:1003014/1‑92 (MQ=37)
                                                tACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGgaga                                       >  2:2319984/1‑93 (MQ=38)
                                                tACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGgaga                                       >  2:588433/1‑93 (MQ=38)
                                                tACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGgaga                                       >  1:2350991/1‑93 (MQ=38)
                                                 aCGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGaa                                      >  2:648745/1‑93 (MQ=38)
                                                 aCGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGaa                                      <  1:1228752/93‑1 (MQ=37)
                                                  cGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAat                                     >  2:1614782/1‑93 (MQ=38)
                                                   ggTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAata                                    <  1:1467268/93‑1 (MQ=37)
                                                    gTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAatat                                   >  1:2619962/1‑93 (MQ=38)
                                                    gTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAatat                                   >  1:2534508/1‑93 (MQ=38)
                                                    gTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAatat                                   <  2:4309117/93‑1 (MQ=37)
                                                    gTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAatat                                   >  2:2240690/1‑93 (MQ=38)
                                                    gTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAatat                                   <  1:2093027/93‑1 (MQ=37)
                                                        gAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTTGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTcc                               <  1:68055/93‑1 (MQ=37)
                                                        gAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGGGAATATCTcc                               >  1:215802/1‑93 (MQ=37)
                                                        gAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTcc                               >  2:3085291/1‑93 (MQ=37)
                                                         aaTTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGACGGAGAATATCTCCg                              >  1:68080/1‑93 (MQ=37)
                                                          aTTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCgg                             >  2:4545163/1‑93 (MQ=38)
                                                           tttGCACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGc                            >  1:1178409/1‑93 (MQ=38)
                                                              gCACGAAATCGTAGGGCAGAGAAGGCGTTCA‑‑CGGCGCATCTGGCAACGACCTCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGTAtt                         <  2:1677887/93‑1 (MQ=255)
                                                               cACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAttt                        <  1:4647581/93‑1 (MQ=38)
                                                               cACGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAttt                        >  2:385177/1‑93 (MQ=38)
                                                                aCGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAtttt                       <  2:4617142/93‑1 (MQ=38)
                                                                aCGAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAtttt                       >  1:4438533/1‑93 (MQ=38)
                                                                  gAAATCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTa                     <  2:846267/93‑1 (MQ=38)
                                                                   aaaTCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTAt                    <  2:798027/93‑1 (MQ=38)
                                                                     aTCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGGGAGTATTTCCGGCAGTTTTAttt                  >  2:1799177/1‑92 (MQ=25)
                                                                     aTCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTc                  >  2:2417101/1‑93 (MQ=38)
                                                                     aTCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTc                  <  1:2307799/93‑1 (MQ=38)
                                                                      tCGTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCCTTTTTATTcc                 >  2:283638/1‑93 (MQ=37)
                                                                        gTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCcac               >  1:3072790/1‑93 (MQ=38)
                                                                        gTAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCcac               >  2:4722882/1‑93 (MQ=38)
                                                                         tAGGGCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCcaca              >  2:3401879/1‑93 (MQ=38)
                                                                           gggCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGc            <  1:1650669/93‑1 (MQ=38)
                                                                           gggCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGc            >  1:159032/1‑93 (MQ=38)
                                                                           gggCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGc            >  2:840891/1‑93 (MQ=38)
                                                                           gggCAGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGc            >  2:1554270/1‑93 (MQ=38)
                                                                               aGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGGATTTTTTTTCCACAGCCaaa        >  2:2568513/1‑93 (MQ=38)
                                                                               aGATAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCaaa        >  1:3104053/1‑93 (MQ=38)
                                                                                aaTAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAAc       <  1:2495954/92‑1 (MQ=38)
                                                                                  tAAGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGCGAATATCTCCGGCATTTTTATTCCACAGCCAAACTc     >  1:3716302/1‑93 (MQ=38)
                                                                                  ataGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTc     <  2:3259442/91‑1 (MQ=38)
                                                                                   aaGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCa    <  1:2431806/93‑1 (MQ=39)
                                                                                    aGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCat   <  1:333291/93‑1 (MQ=39)
                                                                                    aGGCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCat   >  1:3960820/1‑93 (MQ=255)
                                                                                     ggCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCata  >  2:2684515/1‑93 (MQ=255)
                                                                                     ggCGTTCACGC‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCata  >  2:2089310/1‑93 (MQ=255)
                                                                                             ||                                                                                     
GCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑C‑‑CGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATA  >  NC_000913/3537312‑3537487

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