Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,348,753 Δ1 bp coding (240/255 nt) yfaE → ferredoxin involved with ribonucleotide reductase diferric‑tyrosyl radical (Y*) cofactor maintenance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,348,7520T.98.5% 564.5 / 1.1 129coding (239/255 nt)yfaEferredoxin involved with ribonucleotide reductase diferric‑tyrosyl radical (Y*) cofactor maintenance
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base . (68/59);  total (69/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

Only 100 of 131 total aligned reads displayed.

CACGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTC  >  NC_000913/2348653‑2348846
                                                                                                   |                                                                                              
cacGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGatat                                                                                                >  1:3055834/1‑100 (MQ=255)
  cGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTTCCGGGCAAAAGGCGatatga                                                                                              >  2:282586/1‑98 (MQ=255)
        gTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAGAGGCGATATGAAATCga                                                                                       <  2:779136/100‑1 (MQ=255)
         ttGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCgag                                                                                      <  2:2334158/100‑1 (MQ=255)
          tGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCgaga                                                                                     >  2:1517949/1‑100 (MQ=255)
          tGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCgaga                                                                                     <  2:2644580/100‑1 (MQ=255)
          tGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCgaga                                                                                     <  2:2786465/100‑1 (MQ=255)
          ccgatctCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCgaga                                                                                     <  2:81093/94‑1 (MQ=255)
           gCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGAt                                                                                    <  1:2741718/100‑1 (MQ=255)
              ggTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGAtgtg                                                                                 >  1:1692959/1‑100 (MQ=255)
              ggTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGAtgtg                                                                                 >  1:3104499/1‑100 (MQ=255)
              ggTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGAtgtg                                                                                 >  1:544366/1‑100 (MQ=255)
                 cAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCgagaagatcgg                                                                              >  1:81095/1‑93 (MQ=255)
                    gTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTgg                                                                           >  2:794726/1‑100 (MQ=255)
                     ttGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGt                                                                          <  1:2546400/100‑1 (MQ=255)
                      ttccgatcttGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGtt                                                                         <  1:250094/92‑1 (MQ=255)
                      tGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGtt                                                                         >  2:1761363/1‑100 (MQ=255)
                       gACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTg                                                                        <  2:3026413/100‑1 (MQ=255)
                           ggATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGt                                                                    >  2:3138185/1‑100 (MQ=255)
                           ggATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGt                                                                    >  1:70561/1‑100 (MQ=255)
                           ggATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATCAAATCGAGATGTGAATTGGTTGTAGt                                                                    >  1:2973874/1‑100 (MQ=255)
                               tGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCa                                                                <  2:2071731/100‑1 (MQ=255)
                               tGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCa                                                                <  1:2278986/100‑1 (MQ=255)
                               tGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCa                                                                >  1:2314420/1‑100 (MQ=255)
                                gCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAg                                                               >  2:2060618/1‑100 (MQ=255)
                                gCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAg                                                               >  2:2706599/1‑100 (MQ=255)
                                gCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAg                                                               >  1:789024/1‑100 (MQ=255)
                                gCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGGCCAGa                                                              >  2:368479/1‑100 (MQ=255)
                                 ccGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGa                                                              >  1:3108144/1‑100 (MQ=255)
                                    aaCCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGCTGTGAATTGGTTGTAGTGCCAGATAc                                                           >  2:1125341/1‑100 (MQ=255)
                                       cGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAAc                                                        <  2:2532755/100‑1 (MQ=255)
                                       cGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAAc                                                        <  2:2159681/100‑1 (MQ=255)
                                        gTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACg                                                       >  1:1203702/1‑100 (MQ=255)
                                         ttAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCGAGATACAACGc                                                      <  2:2506641/100‑1 (MQ=255)
                                         ttAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGc                                                      <  1:1149848/100‑1 (MQ=255)
                                         ttAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGc                                                      >  1:1474708/1‑100 (MQ=255)
                                           aGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCtt                                                    >  1:2803947/1‑100 (MQ=255)
                                            gCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTa                                                   <  2:3185294/100‑1 (MQ=255)
                                            gCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTa                                                   >  1:142963/1‑100 (MQ=255)
                                            gCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTa                                                   <  2:670088/100‑1 (MQ=255)
                                            gCCTTTATTCAGCCGGGGGAAATTTTGCCCCGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTa                                                   <  1:26646/100‑1 (MQ=255)
                                             ccTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTAt                                                  >  1:895229/1‑100 (MQ=255)
                                              cTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTAtn                                                 <  1:60476/100‑2 (MQ=255)
                                              cTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTAtn                                                 <  1:239897/100‑2 (MQ=255)
                                              cTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATg                                                 >  1:142371/1‑100 (MQ=255)
                                              cTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATg                                                 >  2:1269522/1‑100 (MQ=255)
                                               tttATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGc                                                <  1:1944790/100‑1 (MQ=255)
                                                ttATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCg                                               <  2:600653/100‑1 (MQ=255)
                                                 tATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGt                                              >  2:2603316/1‑100 (MQ=255)
                                                 tATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGt                                              >  1:1587031/1‑100 (MQ=255)
                                                     cAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTa                                          >  2:265443/1‑100 (MQ=255)
                                                     cAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTa                                          >  1:1266106/1‑100 (MQ=255)
                                                     cAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTa                                          <  2:2875910/100‑1 (MQ=255)
                                                     cAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGAGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCCTa                                          >  2:1787067/1‑100 (MQ=255)
                                                      aGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTa                                          <  2:3157244/99‑1 (MQ=255)
                                                      aGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTAt                                         >  1:830640/1‑100 (MQ=255)
                                                      aGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTAt                                         >  1:3264391/1‑100 (MQ=255)
                                                      aGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTAt                                         <  1:265449/100‑1 (MQ=255)
                                                       gCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATc                                        >  2:2181407/1‑100 (MQ=255)
                                                           ggggAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATAAGACGCATAAGCGTTGTATCTGGc                                    <  2:752385/100‑1 (MQ=25)
                                                            gggAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGcc                                   >  1:2449581/1‑100 (MQ=255)
                                                            gggAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGcc                                   >  2:347450/1‑100 (MQ=255)
                                                            gggAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGcc                                   <  1:2333742/100‑1 (MQ=255)
                                                            gggAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGcc                                   <  1:2172725/100‑1 (MQ=255)
                                                             ggAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCt                                  <  2:238559/100‑1 (MQ=255)
                                                             ggAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCt                                  >  2:387203/1‑100 (MQ=255)
                                                              gAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTa                                 >  2:2358652/1‑100 (MQ=255)
                                                               aaaTTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTAc                                <  1:194655/100‑1 (MQ=255)
                                                                aaTTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACa                               >  1:2091294/1‑100 (MQ=255)
                                                                 aTTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACaa                              >  1:1833547/1‑100 (MQ=255)
                                                                  ttttGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAAc                             <  2:3108107/100‑1 (MQ=255)
                                                                   tttGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAAc                             <  1:927459/100‑1 (MQ=255)
                                                                     tGCCCTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGAtt                         <  2:1322534/100‑1 (MQ=255)
                                                                       cccTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTa                        >  1:1999310/1‑100 (MQ=255)
                                                                       cccTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTa                        <  2:1065657/100‑1 (MQ=255)
                                                                       cccTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTa                        >  1:2457759/1‑100 (MQ=255)
                                                                       cccTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTa                        <  2:1474693/100‑1 (MQ=255)
                                                                       cccTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTa                        <  1:347457/100‑1 (MQ=255)
                                                                        ccTGTTGTTGCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTAc                       >  2:1695832/1‑100 (MQ=255)
                                                                             tgtttCCGGGTAAAAGGCGATATGAAATCGAGATGTGAATTGGTTATAGTGACAGATATAACGCTTATGCGTTTTCTCTGGGCTAGATTGATTACGTgg                   >  2:2063579/1‑99 (MQ=255)
                                                                              gttgCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATgg                   >  1:1794776/1‑98 (MQ=255)
                                                                              gttgCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCg                 >  1:1001156/1‑100 (MQ=255)
                                                                              gttgCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCg                 >  2:2328248/1‑100 (MQ=255)
                                                                                tgctgggCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCgta               <  1:434629/96‑1 (MQ=255)
                                                                                tgCTGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCgta               <  1:434429/100‑1 (MQ=255)
                                                                                tgCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCgta               >  2:409133/1‑100 (MQ=255)
                                                                                 gCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCgtag              <  1:3023561/100‑1 (MQ=255)
                                                                                 gCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCgtag              >  2:1260872/1‑100 (MQ=255)
                                                                                 gCCGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCgtag              >  1:759634/1‑100 (MQ=255)
                                                                                  ccGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCgtagt             >  2:1575680/1‑100 (MQ=255)
                                                                                  ccGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCgtagt             >  2:1204239/1‑100 (MQ=255)
                                                                                  ccGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCgtagt             >  1:2517896/1‑100 (MQ=255)
                                                                                   cGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCgtagta            >  1:1767356/1‑100 (MQ=255)
                                                                                   cGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCgtagta            >  1:3105179/1‑100 (MQ=255)
                                                                                   cGGGCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCgtagta            <  1:160411/100‑1 (MQ=255)
                                                                                     ggCAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAAt          >  1:3045367/1‑100 (MQ=255)
                                                                                       cAAAAGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATAc        >  1:178921/1‑100 (MQ=255)
                                                                                         aaaGGCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGc      >  2:2472393/1‑100 (MQ=255)
                                                                                             gCGATATGAAATCGAGATGTGGATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTATTACTACGccatg  <  2:3055798/100‑5 (MQ=255)
                                                                                             gCGATATGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTc  <  1:3241503/100‑1 (MQ=255)
                                                                                                   |                                                                                              
CACGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTTGCCGGGCAAAAGGCGATATTGAAATCGAGATGTGAATTGGTTGTAGTGCCAGATACAACGCTTATGCGTCTTATCTGGCCTACAACGATTACATGGCGTAGTAATACGCTTTC  >  NC_000913/2348653‑2348846

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