Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,656,391 G→A 13.4% intergenic (‑33/+145) pepB ← / ← iscX aminopeptidase B/Fe(2+) donor and activity modulator for cysteine desulfurase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,656,3910GA13.4% 315.1 / 30.8 119intergenic (‑33/+145)pepB/iscXaminopeptidase B/Fe(2+) donor and activity modulator for cysteine desulfurase
Reads supporting (aligned to +/- strand):  ref base G (60/43);  new base A (10/6);  total (70/49)
Fisher's exact test for biased strand distribution p-value = 7.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.29e-01

Only 100 of 120 total aligned reads displayed.

TCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGACGCGGAAATTATT  >  NC_000913/2656243‑2656540
                                                                                                                                                    |                                                                                                                                                     
tCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGt                                                                                                                                                      >  2:1170541/1‑150 (MQ=255)
tCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGt                                                                                                                                                      >  2:844731/1‑150 (MQ=255)
      ccGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAg                                                                                                                                                >  1:472940/1‑150 (MQ=255)
       cGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGc                                                                                                                                               >  2:1095216/1‑150 (MQ=255)
       cGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATACTTCTTTGTAAAAGGGCCGCCAGTTGGCAGc                                                                                                                                               >  2:464729/1‑150 (MQ=255)
         ttCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGccc                                                                                                                                             >  1:1246576/1‑150 (MQ=255)
          tCAACTGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGcccc                                                                                                                                            <  2:234741/150‑1 (MQ=255)
            aaaTGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAg                                                                                                                                          >  1:164542/1‑150 (MQ=255)
             aaTGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAgg                                                                                                                                         <  1:596872/150‑1 (MQ=255)
               tGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCTGCCAACTGGCGGcccttttac                                                                                                                                       <  1:543777/150‑7 (MQ=255)
               tGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGGTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCTGCCAACTGGCGGcccttttac                                                                                                                                       <  1:1341018/150‑7 (MQ=255)
               tGCAGGGTAATGCCGTAATTATAAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGGCGGCAGGTTGGGTAGAGAGGGTAAGCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCTGCCAACTGGCGGcccttttac                                                                                                                                       <  1:662509/150‑7 (MQ=255)
               tGCACGGTAATGCCGGCATTATTAATGCTGTAAGTGGCGTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCTGCCAACTGGCGGcccttttac                                                                                                                                       <  1:1106489/150‑7 (MQ=255)
                gCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGTGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGAtt                                                                                                                                      <  1:567289/150‑1 (MQ=255)
                gCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCTGCCAACTGGCGGcccttttaca                                                                                                                                      <  1:406125/150‑8 (MQ=255)
                 cAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCTTCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGAttt                                                                                                                                     >  2:724227/1‑150 (MQ=255)
                  aGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGAtttt                                                                                                                                    >  2:622376/1‑150 (MQ=255)
                   gggTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGAttttt                                                                                                                                   >  2:1459302/1‑150 (MQ=255)
                   gggTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGAttttt                                                                                                                                   <  1:844731/150‑1 (MQ=255)
                             gTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATg                                                                                                                         <  2:1573461/150‑1 (MQ=255)
                                attattAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGAtt                                                                                                                      >  2:589015/1‑150 (MQ=255)
                                  tattaATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTcg                                                                                                                    >  2:788005/1‑150 (MQ=255)
                                    ttaATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTcgcg                                                                                                                  >  1:1535289/1‑150 (MQ=255)
                                    ttaATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTcgcg                                                                                                                  <  1:622376/150‑1 (MQ=255)
                                      aaTGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGtt                                                                                                                >  2:359696/1‑150 (MQ=255)
                                       aTGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTCTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTAGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTg                                                                                                               >  1:1673935/1‑150 (MQ=255)
                                         tcTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGGTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCa                                                                                                             <  1:149392/149‑1 (MQ=255)
                                         gCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCa                                                                                                             >  1:1019575/1‑150 (MQ=255)
                                         gCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCa                                                                                                             >  2:1627373/1‑150 (MQ=255)
                                         gCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCa                                                                                                             >  1:1054180/1‑150 (MQ=255)
                                         gCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCa                                                                                                             <  2:435152/150‑1 (MQ=255)
                                           tGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAACCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAgg                                                                                                           <  1:1480887/150‑1 (MQ=255)
                                                  tGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCTTCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAggcgg                                                                                                    <  1:724227/150‑1 (MQ=255)
                                                    cTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCa                                                                                                  >  2:864252/1‑150 (MQ=255)
                                                    cTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCa                                                                                                  >  1:707421/1‑150 (MQ=255)
                                                     tttttCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCaa                                                                                                 <  1:1244848/150‑1 (MQ=255)
                                                     tttttCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCaa                                                                                                 <  1:1255025/150‑1 (MQ=255)
                                                            ccccAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAagtgagtg                                                                                          >  1:1699492/1‑150 (MQ=255)
                                                            ccccAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAagtgagtg                                                                                          >  1:1634600/1‑150 (MQ=255)
                                                            ccccAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAagtgagtg                                                                                          <  2:809417/150‑1 (MQ=255)
                                                              ccAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGaa                                                                                        >  2:1707225/1‑150 (MQ=255)
                                                                  cgcgcgTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGccc                                                                                    <  1:728071/150‑1 (MQ=255)
                                                                  cgcgcgTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGccc                                                                                    <  2:1246576/150‑1 (MQ=255)
                                                                   gcgcgTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGcccc                                                                                   >  1:974581/1‑150 (MQ=255)
                                                                    cgcgTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCa                                                                                  >  2:1316750/1‑150 (MQ=255)
                                                                    cgcgTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCa                                                                                  <  2:1019575/150‑1 (MQ=255)
                                                                        gcGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAg                                                                              <  2:164542/149‑1 (MQ=255)
                                                                                ttGGGTAGAGAGGGTGATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATa                                                                      >  2:678804/1‑150 (MQ=255)
                                                                                  gggTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCTCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAATCCCCAGGAGCATAGATaaa                                                                    >  2:185681/1‑149 (MQ=255)
                                                                                   ggTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAATTGGCAGCCCCAGGATTTTTACGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACt                                                                   >  2:1698411/1‑150 (MQ=255)
                                                                                       gagagGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTAtgt                                                               >  1:995653/1‑150 (MQ=255)
                                                                                             gTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTgg                                                         >  2:1147366/1‑150 (MQ=255)
                                                                                               aaTCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTgggg                                                       <  2:1008497/150‑1 (MQ=255)
                                                                                               aaTCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTgggg                                                       <  1:1406116/150‑1 (MQ=255)
                                                                                               aaTCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTgggg                                                       <  1:1316750/150‑1 (MQ=255)
                                                                                                ttCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGTCAGCCCCAGGATTTATCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGt                                                      <  1:302250/149‑1 (MQ=255)
                                                                                                 tCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTg                                                     <  1:1707225/150‑1 (MQ=255)
                                                                                                 tCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTg                                                     <  1:1095216/150‑1 (MQ=255)
                                                                                                  cTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGa                                                    >  2:1091586/1‑150 (MQ=255)
                                                                                                    tCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGGTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAAc                                                  <  1:999443/150‑1 (MQ=255)
                                                                                                       tCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAg                                               >  1:888491/1‑150 (MQ=255)
                                                                                                         gCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGGTGCAGGACGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGCGACTGGGGTGAACGAgcg                                             <  1:589015/150‑1 (MQ=255)
                                                                                                           ttCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCa                                           <  1:1539770/150‑1 (MQ=255)
                                                                                                            tCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAg                                          >  1:927562/1‑150 (MQ=255)
                                                                                                            gcTGTCATTTAACTTATCCTTCGTTGTATAAGGGCCGCCAGGTGGCAGCCCCAGGCTATGTCCGCTAAATGCTGCGCGTGGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATCACTATGTGACTGGGGTGAACGAGCGCAg                                          <  1:1259835/149‑1 (MQ=255)
                                                                                                             cTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGc                                         >  2:1711597/1‑150 (MQ=255)
                                                                                                                 cATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAAc                                     <  1:1479737/150‑1 (MQ=255)
                                                                                                                 cATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAAc                                     >  1:177056/1‑150 (MQ=255)
                                                                                                                  aTTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACg                                    >  2:208074/1‑150 (MQ=255)
                                                                                                                   ttttAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGc                                   >  2:1059835/1‑150 (MQ=255)
                                                                                                                     ttAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGTCAGCCCCAGGATATTTCCGCTAAATGATTCGCATTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCAt                                 <  2:463424/150‑1 (MQ=255)
                                                                                                                        gTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTg                              >  2:88961/1‑150 (MQ=255)
                                                                                                                         ttATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGCGCGCAGCCAACGCATCTgc                             >  2:668038/1‑150 (MQ=255)
                                                                                                                             ccTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTgcggcg                         >  2:1092331/1‑150 (MQ=255)
                                                                                                                              cttcttTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGt                        >  2:7926/1‑150 (MQ=255)
                                                                                                                              cttcttTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGt                        >  1:294101/1‑150 (MQ=255)
                                                                                                                                tcttTGTATAAGGGCCTCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGa                      <  2:670059/150‑1 (MQ=255)
                                                                                                                                  tttGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAg                    >  1:1586916/1‑150 (MQ=255)
                                                                                                                                   ttGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGc                   >  2:45668/1‑150 (MQ=255)
                                                                                                                                     gTAAACGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGAATCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATTACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCAt                 >  2:736799/1‑150 (MQ=255)
                                                                                                                                     gTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCTCAGCCAACGCATCTGCGGCGTGatgaat                 >  2:1492311/1‑145 (MQ=255)
                                                                                                                                      tAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATg                <  2:1579802/150‑1 (MQ=255)
                                                                                                                                      tAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATg                <  2:472940/150‑1 (MQ=255)
                                                                                                                                      tAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATg                <  1:1192642/150‑1 (MQ=255)
                                                                                                                                       cctggggCTGCCAACTGGCGGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGCGCGCAGCCAACGCATCTGCGGCGTGAAGCATGa               >  2:1407809/5‑150 (MQ=38)
                                                                                                                                       cctggggCTGCCAACTGGCGGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCCGCCAACGCATCTGCGGCGTGAAGCATGa               >  1:300730/5‑150 (MQ=38)
                                                                                                                                       cctggggCTGCCAACTGGCGGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGACGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGa               >  2:931270/5‑150 (MQ=38)
                                                                                                                                       aaaaGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGAGAGTGAAGCCCCAGGAGCATAGATAACTATTTGAATGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCAtta               >  2:260933/1‑148 (MQ=255)
                                                                                                                                        ctggggCTGCCAACTGGCGGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGAc              >  2:1322297/4‑150 (MQ=38)
                                                                                                                                        ctggggCTGCCAACTGGCGGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGAc              >  2:131341/4‑150 (MQ=38)
                                                                                                                                        ctggggCTGCCAACTGGCGGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGAc              >  1:921668/4‑150 (MQ=38)
                                                                                                                                        aaaGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGa                      >  2:1331180/1‑142 (MQ=255)
                                                                                                                                        aaaGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGAc              >  2:849260/1‑150 (MQ=255)
                                                                                                                                          ggggCTGCCAACTGGCGGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGAcgc            >  1:1375563/2‑150 (MQ=255)
                                                                                                                                            ggCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAAAGATTCGCGTGCCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGAAAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGACGCgg          <  1:1427807/150‑1 (MQ=255)
                                                                                                                                              ccgccAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATTACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGACGCGGaa        >  1:1629087/1‑150 (MQ=255)
                                                                                                                                              ccgccAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATTACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTAACGCATGACGCGGaa        >  1:1622432/1‑150 (MQ=255)
                                                                                                                                              ccgccAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGACGCGGaa        >  1:191906/1‑150 (MQ=255)
                                                                                                                                               cgccAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGACGCGGaaa       >  1:1361471/1‑150 (MQ=255)
                                                                                                                                                    gTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGACGCGGAAattatt  <  1:678804/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTCTTTGTAAAAGGGCCGCCAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCATGACGCGGAAATTATT  >  NC_000913/2656243‑2656540

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