Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,539,817 A→T 5.3% intergenic (+111/‑346) yhgF → / → feoA putative transcriptional accessory protein/ferrous iron transporter, protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,539,8170AT 5.3% 988.8 / 31.2 272intergenic (+111/‑346)yhgF/feoAputative transcriptional accessory protein/ferrous iron transporter, protein A
Reads supporting (aligned to +/- strand):  ref base A (131/126);  new base T (7/8);  total (138/134)
Fisher's exact test for biased strand distribution p-value = 7.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 276 total aligned reads displayed.

ATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAACCTTGAGC  >  NC_000913/3539668‑3539965
                                                                                                                                                     |                                                                                                                                                    
ttggatgCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAata                                                                                                                                                      <  2:1702618/149‑1 (MQ=255)
  ggatgCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGTATAtc                                                                                                                                                    >  2:2546956/1‑150 (MQ=255)
  ggatgCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATAtc                                                                                                                                                    >  2:922857/1‑150 (MQ=255)
     tgCGCTGGCGGCGGCATTGGGCAAAATACTTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGATATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTcc                                                                                                                                                 <  2:3141250/150‑1 (MQ=255)
      gcgcTGGCGGCGGCAATGGGCAAAACACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCg                                                                                                                                                <  2:3360411/150‑1 (MQ=255)
        gcTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGc                                                                                                                                              >  1:2376737/1‑150 (MQ=255)
             cggcggcAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAtttt                                                                                                                                         <  1:3076575/150‑1 (MQ=255)
              ggcggcAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAttttt                                                                                                                                        >  2:2239585/1‑150 (MQ=255)
                cggcAATGGGAAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTAt                                                                                                                                      <  1:2808055/150‑1 (MQ=255)
                  gcAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTc                                                                                                                                    <  1:1419334/150‑1 (MQ=255)
                   cAATGGGCAAACAACTTAAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAACTCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAATGCCGGAGATATTCTCCGGCACTTTTAtgtc                                                                                                                                   <  1:1867590/150‑4 (MQ=38)
                   cAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAATGCCGGAGATATTCTCCGGCACTTTTAtgtc                                                                                                                                   <  2:581260/150‑4 (MQ=39)
                    attGGCCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGGAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGATATTCTCCGGCACTTTTAtgtca                                                                                                                                  <  2:3309696/148‑5 (MQ=38)
                    aaTGGGCAAACACCGTTAAACGCCCGTACTGGCCTACGGTTCGAATGTGCACGAACTCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGATATTCTCCGGCACTTTTAtgtca                                                                                                                                  <  2:3556916/150‑5 (MQ=38)
                    aaTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCca                                                                                                                                  <  1:4295963/150‑1 (MQ=255)
                      tGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCcaca                                                                                                                                <  2:3984596/150‑1 (MQ=255)
                       gggCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAg                                                                                                                               >  1:778735/1‑150 (MQ=255)
                       gggCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAg                                                                                                                               >  1:2061216/1‑150 (MQ=255)
                        ccaaaaaaaCCTAAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACCAAATCGTAGGGCAGATAAGGCGTACACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGc                                                                                                                              <  1:1128852/147‑1 (MQ=255)
                            aaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCaaa                                                                                                                          >  2:2757464/1‑150 (MQ=255)
                            aaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCaaa                                                                                                                          <  2:1813412/150‑1 (MQ=255)
                                  gTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCata                                                                                                                    >  1:3122445/1‑150 (MQ=255)
                                      caCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAatat                                                                                                                <  2:2245710/149‑1 (MQ=255)
                                       aCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAatata                                                                                                               >  1:3471282/1‑150 (MQ=255)
                                           ccGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTcc                                                                                                           >  2:1464513/1‑150 (MQ=255)
                                            cGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCg                                                                                                          <  2:2491403/150‑1 (MQ=255)
                                                     cTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATAttt                                                                                                 <  1:2308920/150‑1 (MQ=255)
                                                       tcGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATAAATTCCGGCAATATTTAt                                                                                               <  2:4056955/149‑1 (MQ=255)
                                                           ttCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcatt                                                                                           >  2:1365002/1‑150 (MQ=255)
                                                              gAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcatttca                                                                                        <  2:114197/150‑1 (MQ=255)
                                                               aaTTTGCACGACATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAACAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcatttcat                                                                                       <  2:1857660/150‑1 (MQ=255)
                                                                   tGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTaac                                                                                   <  1:2734925/150‑1 (MQ=255)
                                                                   tGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTaac                                                                                   <  1:3650806/150‑1 (MQ=255)
                                                                     cACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATACTATATTCCGGCAATATTTATCATTTCATTaacaa                                                                                 >  2:2309701/1‑150 (MQ=255)
                                                                      aCGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTaacaac                                                                                <  1:1938837/150‑1 (MQ=255)
                                                                        gAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTg                                                                              >  1:3286715/1‑150 (MQ=255)
                                                                        gAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTg                                                                              >  2:3429480/1‑150 (MQ=255)
                                                                          aaTCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGaa                                                                            <  2:2637690/150‑1 (MQ=255)
                                                                           aTCGTAGGGCAGTTAAGGCGTTCACGCCGCATCTGGCATCGAACGCCTTGCCTGACATAAAAGTGCCGTAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGaaa                                                                           <  2:867934/150‑1 (MQ=255)
                                                                           aTCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGaaa                                                                           <  2:2385396/150‑1 (MQ=255)
                                                                            tCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAAc                                                                          <  1:4193450/150‑1 (MQ=255)
                                                                              gTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATCTATTCCGGCACTCTTTATCATTTCATTAACAACTGAAACCt                                                                        >  2:3031069/1‑150 (MQ=255)
                                                                               tAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCtt                                                                       >  2:2815654/1‑150 (MQ=255)
                                                                                 gggCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCttaa                                                                     <  1:330463/150‑1 (MQ=255)
                                                                                  ggCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCttaat                                                                    >  1:2914440/1‑150 (MQ=255)
                                                                                   gCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCttaatt                                                                   >  2:1537585/1‑150 (MQ=255)
                                                                                    cAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCttaatta                                                                  <  1:4139829/150‑1 (MQ=255)
                                                                                     aGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGACTATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCttaattaa                                                                 >  2:3218230/1‑150 (MQ=255)
                                                                                     aGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCttaattaa                                                                 >  1:117122/1‑150 (MQ=255)
                                                                                      gATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATAATTATCATTTCATTAACAACTGAAACCTTAATTaaa                                                                <  1:1794181/150‑1 (MQ=255)
                                                                                        tAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACa                                                              >  2:3693271/1‑150 (MQ=255)
                                                                                        tAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACa                                                              >  1:2449692/1‑150 (MQ=255)
                                                                                          aGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACAtt                                                            >  1:1208739/1‑150 (MQ=255)
                                                                                           ggCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACCTTa                                                           >  1:392197/1‑150 (MQ=255)
                                                                                           ggCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTa                                                           <  1:985898/150‑1 (MQ=255)
                                                                                           ggCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTa                                                           <  1:1258146/150‑1 (MQ=255)
                                                                                           ggCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTa                                                           >  1:1701464/1‑150 (MQ=255)
                                                                                              gTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCATAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGcc                                                        >  2:28874/1‑150 (MQ=255)
                                                                                              gTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCATAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGcc                                                        >  2:2239830/1‑150 (MQ=255)
                                                                                              gTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCCAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACAGTAGcc                                                        >  1:1555689/1‑150 (MQ=255)
                                                                                               ttCACTCCGCATCTTGCAAAGAACGCCTTGCCTGACATAAAATTGCCTGAGAATATCTCCGGCATTTTTATTCCACATCCAAACTCATAATTTATTCCGGAAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCa                                                       <  1:560890/150‑1 (MQ=255)
                                                                                               ttCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCa                                                       <  1:3211305/150‑1 (MQ=255)
                                                                                               ttCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCa                                                       <  1:1755935/150‑1 (MQ=255)
                                                                                                  aCGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAAGTCATAATATATTCCGGCAATATTTATCATTTCATTATAAAGTGAAACCTTAATTAAACATTAGCCAGTc                                                    >  1:2534739/1‑150 (MQ=255)
                                                                                                  aCGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTc                                                    >  1:3997817/1‑150 (MQ=255)
                                                                                                  aCGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTc                                                    >  1:2975878/1‑150 (MQ=255)
                                                                                                  aCGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTc                                                    >  1:2282115/1‑150 (MQ=255)
                                                                                                  aCGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTc                                                    >  2:538993/1‑150 (MQ=255)
                                                                                                    gccgcATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATAATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTcag                                                  >  2:3320210/1‑148 (MQ=255)
                                                                                                      cgcATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCggg                                                >  2:3616349/1‑150 (MQ=255)
                                                                                                      cgcATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAATTTTATTCCACAGCCAAACTCATAATATATTCCCGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCgtg                                                <  2:2121717/150‑3 (MQ=255)
                                                                                                            tGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTc                                          >  1:1292779/1‑150 (MQ=255)
                                                                                                            tGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTc                                          >  1:1586157/1‑150 (MQ=255)
                                                                                                            tGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTc                                          >  1:2103765/1‑150 (MQ=255)
                                                                                                             ggCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCa                                         <  2:446629/150‑1 (MQ=255)
                                                                                                                 acgaacgCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTAt                                     <  1:480581/150‑1 (MQ=255)
                                                                                                                   gaccgCCTGGCCTGACATAACAGGGCCGGAGAATACCTCCGGCATGTTTATTCCACAGCCAAACTCATAATATCGTCCGGCAATATTTATCATTTCATGAACAACTGCAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTc                                   <  1:2403145/147‑1 (MQ=255)
                                                                                                                   gaacgCCTTGCCTGACATAAATGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTc                                   >  1:1462387/1‑150 (MQ=255)
                                                                                                                     acgCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGa                                 >  1:922579/1‑150 (MQ=255)
                                                                                                                       gCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAAt                               <  1:2261752/150‑1 (MQ=255)
                                                                                                                       gCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAAt                               <  2:1228068/150‑1 (MQ=255)
                                                                                                                         cTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATta                             >  1:1089503/1‑150 (MQ=255)
                                                                                                                         cTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATta                             >  2:759105/1‑150 (MQ=255)
                                                                                                                             ccTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATAtt                         >  1:257825/1‑150 (MQ=255)
                                                                                                                                 gaaTAAAAATGCCGGAGATATTCTCCGGCACTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCg                     >  1:2755554/3‑150 (MQ=39)
                                                                                                                                   aTAAAAATGCCGGAGATATTCTCCGGCACTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCt                   >  1:3696895/1‑150 (MQ=255)
                                                                                                                                    tAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTACACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTg                  <  1:2577457/150‑1 (MQ=255)
                                                                                                                                     aaaaGTGCCGGAGAATATCTCCGGCATATTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACACCTGAAACCTTAATTAAACATTAGCCAGGCCGGGTCCTTCACGCTTCGAAGTATATTGTCGCTGc                 >  2:2673376/1‑150 (MQ=255)
                                                                                                                                        aatGCCGGAGATATTCTCCGGCACTTTTATTCCACAGCCAATCTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAATCCTTAATTAATCATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGat              >  1:2004542/3‑150 (MQ=255)
                                                                                                                                          tGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGatat            <  1:3429480/150‑1 (MQ=255)
                                                                                                                                          tGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGatat            <  1:334484/150‑1 (MQ=255)
                                                                                                                                          tGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGatat            >  2:1917710/1‑150 (MQ=255)
                                                                                                                                          tGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAAGTATATTTTCGCTGCGatat            <  1:606807/150‑1 (MQ=255)
                                                                                                                                               gagaATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAACCt       <  1:3007572/150‑1 (MQ=255)
                                                                                                                                                  aaTATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAATCATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAACCTTGa    <  2:3815009/150‑1 (MQ=255)
                                                                                                                                                  aaTATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAACCTTGa    <  2:1603430/150‑1 (MQ=255)
                                                                                                                                                  aaTATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAACCTTGa    <  2:114518/150‑1 (MQ=255)
                                                                                                                                                  aaTATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAACCTTGa    <  2:1135086/150‑1 (MQ=255)
                                                                                                                                                  aaTATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAACCTTGa    <  2:2306184/150‑1 (MQ=255)
                                                                                                                                                  aaTATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAAACTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAACCTTGa    <  2:4342791/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATGGATGCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAACCTTGAGC  >  NC_000913/3539668‑3539965

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: