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.8% 1501.5 / ‑4.6 468intergenic (+60/‑395)b3407/feoAhypothetical protein/ferrous iron transport protein A
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (308/159);  total (309/159)
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
*NC_0009133,537,4042.G99.8% 1510.1 / ‑4.5 470intergenic (+60/‑395)b3407/feoAhypothetical protein/ferrous iron transport protein A
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (308/161);  total (309/161)
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.87e-01

Only 100 of 505 total aligned reads displayed.

GCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑CCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATA  >  NC_000913/3537312‑3537487
                                                                                             ||                                                                                     
 cgcTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑c                                                                                      <  2:65813/93‑1 (MQ=255)
 cgcTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑c                                                                                      <  2:2249545/93‑1 (MQ=255)
    tGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑cgcc                                                                                   >  1:3500905/1‑90 (MQ=255)
    tGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑cgcc                                                                                   >  1:8794720/1‑90 (MQ=255)
     ggcggcggcAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑cgccg                                                                                  >  1:8625154/1‑89 (MQ=255)
           cgcAATGGGAAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTg                                                                            <  1:3096325/92‑1 (MQ=37)
             cAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGc                                                                          <  1:6689524/93‑1 (MQ=38)
             cAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGc                                                                          <  1:7179072/93‑1 (MQ=37)
              aaTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCa                                                                         <  1:1964352/93‑1 (MQ=37)
              aaTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCa                                                                         >  1:2781773/1‑93 (MQ=37)
               aTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaa                                                                        >  1:9288324/1‑93 (MQ=37)
               aTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaa                                                                        >  2:6729284/1‑93 (MQ=37)
               aTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaa                                                                        >  2:5636941/1‑93 (MQ=37)
               aTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaa                                                                        <  2:1027029/93‑1 (MQ=37)
                tGGGCAATAAGCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGGAGATAAGGCGTTCACGCCGCATCTGGCaac                                                                       <  1:181598/93‑1 (MQ=25)
                tGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaac                                                                       >  2:5941547/1‑93 (MQ=37)
                 gggCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaacg                                                                      >  1:1614904/1‑93 (MQ=37)
                 gggCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaacg                                                                      >  1:4972044/1‑93 (MQ=37)
                  ggCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaacga                                                                     <  2:2028103/93‑1 (MQ=37)
                   gaaaaaacCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaacgaa                                                                    <  2:9479574/91‑1 (MQ=37)
                   gCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaacgaa                                                                    <  2:4751791/93‑1 (MQ=37)
                   gCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaacgaa                                                                    <  2:7840954/93‑1 (MQ=37)
                     aaaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaacgaacg                                                                  >  2:6000164/1‑93 (MQ=37)
                     aaaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaacgaacg                                                                  >  2:9178930/1‑93 (MQ=37)
                      aaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGc                                                                 >  2:8280290/1‑93 (MQ=37)
                       aaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGcc                                                                <  1:6117943/93‑1 (MQ=37)
                        aaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCt                                                               <  1:4590182/93‑1 (MQ=37)
                         aaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCtt                                                              >  1:2388111/1‑93 (MQ=37)
                            gTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGcc                                                           <  2:433591/93‑1 (MQ=37)
                            gTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGcc                                                           >  2:7115630/1‑93 (MQ=37)
                            gTTAAACGCCCGTACTGCCCTACGGTTCGAATTGGCACGAAATCGGAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGcc                                                           <  1:3890505/93‑1 (MQ=37)
                             ttAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCt                                                          <  1:474320/93‑1 (MQ=37)
                             ttAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCt                                                          <  2:9576094/93‑1 (MQ=37)
                               aaaCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGa                                                        <  2:2887397/93‑1 (MQ=37)
                                aaCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGat                                                       >  1:398300/1‑92 (MQ=37)
                                aaCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGAc                                                       >  1:3427969/1‑93 (MQ=38)
                                aaCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGAc                                                       <  2:1400137/93‑1 (MQ=37)
                                    cccGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATaa                                                   >  2:7426289/1‑93 (MQ=37)
                                     ccGTAATGGCCTACGGTTCGAATTTGCAGGAAATCGTATGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATaaa                                                  <  2:3452725/93‑1 (MQ=37)
                                      cGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATaaaa                                                 >  1:9448312/1‑93 (MQ=37)
                                       gTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAg                                                >  2:8085541/1‑93 (MQ=38)
                                          cTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGGGc                                             >  1:7719891/1‑93 (MQ=37)
                                             gCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCgg                                          <  2:8754129/93‑1 (MQ=37)
                                             gCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCgg                                          >  1:4287387/1‑93 (MQ=37)
                                              cttACGGTTCGAATTGCCACGAAATGGTAGGTCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGga                                         <  2:4435146/91‑1 (MQ=21)
                                               cTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGgag                                        >  1:8355014/1‑93 (MQ=38)
                                                tACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGgaga                                       >  1:6894551/1‑93 (MQ=38)
                                                tACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGgaga                                       >  2:7055178/1‑93 (MQ=37)
                                                tACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGgaga                                       >  2:1367792/1‑93 (MQ=38)
                                                tACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGgaga                                       >  1:2331251/1‑93 (MQ=38)
                                                  cGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAat                                     >  2:3004825/1‑93 (MQ=38)
                                                   ggTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAata                                    >  1:595277/1‑93 (MQ=38)
                                                    gTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAatat                                   >  1:117436/1‑93 (MQ=38)
                                                     ttCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATAtc                                  <  1:5823650/93‑1 (MQ=37)
                                                     ttCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATAtc                                  >  1:5602540/1‑93 (MQ=38)
                                                      tCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATAtct                                 >  1:2108196/1‑93 (MQ=38)
                                                       cGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATAtctc                                >  1:5157466/1‑93 (MQ=38)
                                                         aaTTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCg                              >  1:7982183/1‑93 (MQ=38)
                                                         aaTTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAAAATATCCg                              >  2:3771158/1‑93 (MQ=37)
                                                         aaTTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑CGCAGCATCTGGCAACGAACGACTTGCCTGACATAAAAGTGCCGGAGAATATCTCCg                              >  1:5345994/1‑93 (MQ=37)
                                                          aTTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCgg                             <  2:3242719/93‑1 (MQ=38)
                                                           tttGCCCGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGc                            <  1:1522315/93‑1 (MQ=37)
                                                           tttGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGc                            >  1:6261454/1‑93 (MQ=38)
                                                           tttGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGc                            >  1:3750557/1‑93 (MQ=38)
                                                           tttGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGc                            >  2:2628313/1‑93 (MQ=38)
                                                           tttGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGc                            >  1:7536913/1‑93 (MQ=38)
                                                           tttGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGc                            >  2:3926664/1‑93 (MQ=38)
                                                            ttGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCa                           >  2:7337010/1‑93 (MQ=38)
                                                              gCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAtt                         <  2:9287045/93‑1 (MQ=38)
                                                              gCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAtt                         <  1:8433675/93‑1 (MQ=38)
                                                              gCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAtt                         >  2:5173766/1‑93 (MQ=38)
                                                                aCGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAtttt                       >  1:2893838/1‑93 (MQ=38)
                                                                  gAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATCTCTCCAGCATTTTTa                     >  2:1936923/1‑93 (MQ=37)
                                                                  gAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTa                     >  1:1120476/1‑93 (MQ=38)
                                                                  gAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTa                     >  1:3282805/1‑93 (MQ=38)
                                                                   aaaTGGTAGGGAAGATAGGCCTTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTAt                    <  2:9026063/93‑1 (MQ=255)
                                                                   aaaTCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTAt                    >  1:9205479/1‑93 (MQ=38)
                                                                   aaaTCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTAt                    <  2:1642169/93‑1 (MQ=38)
                                                                    aaTCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTAtt                   <  1:1310062/93‑1 (MQ=38)
                                                                     aTCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCCTTTTTATTc                  >  1:6584557/1‑93 (MQ=37)
                                                                        gTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAGCGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCcac               <  1:10004389/93‑1 (MQ=37)
                                                                         tAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCcaca              >  1:5511055/1‑93 (MQ=38)
                                                                          aGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAg             <  2:610348/93‑1 (MQ=38)
                                                                          aGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAg             >  1:7446586/1‑93 (MQ=38)
                                                                           gggCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTTCACAGc            >  2:111863/1‑93 (MQ=38)
                                                                            ggCAGATAAGGCGTGCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGcc           <  1:6036668/93‑1 (MQ=38)
                                                                             gCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCa          >  2:6466043/1‑93 (MQ=38)
                                                                              cAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCaa         >  1:9949577/1‑93 (MQ=38)
                                                                                gATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTTCCTGACATAAAAGTGCCGGCGAATATCTCCGGCATTTTTTTTCCACAGCCAAAc       >  1:8255174/1‑93 (MQ=37)
                                                                                gATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGGATATCTACGGCCTTTTTTTTTTACAGCCCAAc       >  2:5423822/1‑93 (MQ=25)
                                                                                gATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAAc       >  2:759996/1‑93 (MQ=39)
                                                                                gATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAAc       >  1:2553865/1‑93 (MQ=39)
                                                                                gATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTGTTATTCCACAGCCAAAc       >  1:8979334/1‑93 (MQ=38)
                                                                                 aTAAGGCGTTCACGCCGCATCTGGCAACGAACGGCGTGCCTGACATAAAAGTGCCGGAGACTATCTCCGGCATTTTTATTCCACAGCCAAACt      <  2:8341258/93‑1 (MQ=37)
                                                                                 aTAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACt      >  2:9723404/1‑93 (MQ=39)
                                                                                   aaGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCa    >  2:6273807/1‑93 (MQ=39)
                                                                                   aaGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCa    >  2:5765830/1‑93 (MQ=39)
                                                                                    aGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCat   <  1:9637279/93‑1 (MQ=39)
                                                                                    aGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCat   >  2:4927942/1‑93 (MQ=255)
                                                                                     ggCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCata  >  1:4489527/1‑93 (MQ=39)
                                                                                             ||                                                                                     
GCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑CCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATA  >  NC_000913/3537312‑3537487

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