Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,083,693 A→C A106A (GCA→GCC)  efeB → deferrrochelatase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,083,6930AC98.9% 664.1 / ‑3.1 186A106A (GCA→GCC) efeBdeferrrochelatase, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (97/86);  total (99/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 192 total aligned reads displayed.

TGGTGGCGTTTGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGT  >  NC_000913/1083596‑1083788
                                                                                                 |                                                                                               
tggtggCGTTTGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGagccgc                                                                                               <  1:4864448/100‑1 (MQ=255)
 ggtggCGTTTGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAgccgca                                                                                              >  1:1412880/1‑100 (MQ=255)
 ggtggCGTTTGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAgccgca                                                                                              <  1:97582/100‑1 (MQ=255)
   tggCGTTTGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCAcc                                                                                            <  1:3447170/100‑1 (MQ=255)
    ggCGTTTGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCa                                                                                           >  2:2585694/1‑100 (MQ=255)
       gTTTGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGaa                                                                                        >  1:4310911/1‑100 (MQ=255)
        tctGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGaaa                                                                                       <  1:2608865/98‑1 (MQ=255)
         ttGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAAc                                                                                      <  1:5420772/100‑1 (MQ=255)
          tGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACg                                                                                     <  2:2558762/100‑1 (MQ=255)
            aTGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGcc                                                                                   <  2:2176183/100‑1 (MQ=255)
            aTGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGcc                                                                                   <  2:491043/100‑1 (MQ=255)
             tgtgCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCa                                                                                  <  2:2527833/100‑1 (MQ=255)
             tgtgCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCa                                                                                  <  2:1693719/100‑1 (MQ=255)
             tgtgCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCa                                                                                  >  1:2169442/1‑100 (MQ=255)
              gtgCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCaa                                                                                 <  1:516395/100‑1 (MQ=255)
               tgCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCaaa                                                                                <  1:2345543/100‑1 (MQ=255)
                gCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCTGCACCAGAAACGCCAAAt                                                                               >  2:2277641/1‑100 (MQ=255)
                 cTTGCCAGCGATAAAGCCGATCGTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCCGAAACGCCAAATc                                                                              <  2:5226061/100‑1 (MQ=255)
                  ttGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATcc                                                                             <  1:3270008/100‑1 (MQ=255)
                   tGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCcg                                                                            <  2:4830378/100‑1 (MQ=255)
                    gCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCcgc                                                                           >  2:3746930/1‑100 (MQ=255)
                     ccAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCcgcg                                                                          >  2:3672935/1‑100 (MQ=255)
                      cAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCcgcgn                                                                         <  1:4913995/100‑2 (MQ=255)
                      cAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCgcgc                                                                         >  2:4039146/1‑100 (MQ=255)
                      cAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCgcgc                                                                         <  2:2016638/100‑1 (MQ=255)
                       aGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGcc                                                                        >  1:3125065/1‑100 (MQ=255)
                       aGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGcc                                                                        <  2:2364950/100‑1 (MQ=255)
                       aGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGcc                                                                        >  2:3241237/1‑100 (MQ=255)
                        gCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCt                                                                       >  1:1752276/1‑100 (MQ=255)
                        gCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCt                                                                       >  1:1863708/1‑100 (MQ=255)
                         cGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTg                                                                      <  1:590360/100‑1 (MQ=255)
                          gATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGc                                                                     >  2:2901645/1‑100 (MQ=255)
                           aTAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGcc                                                                    >  1:4096425/1‑100 (MQ=255)
                           aTAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGcc                                                                    <  2:5365017/100‑1 (MQ=255)
                           aTAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGcc                                                                    <  2:5005347/100‑1 (MQ=255)
                            tAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGcca                                                                   <  1:4073498/100‑1 (MQ=255)
                             aaaGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGccac                                                                  <  2:3951558/100‑1 (MQ=255)
                              aaGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGccacc                                                                 <  1:3701916/100‑1 (MQ=255)
                              aaGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGccacc                                                                 <  1:2236804/100‑1 (MQ=255)
                              aaGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGccacc                                                                 <  1:3415922/100‑1 (MQ=255)
                              aaGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGccacc                                                                 <  1:2934773/100‑1 (MQ=255)
                               aGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGccacca                                                                >  1:9084/1‑100 (MQ=255)
                               aGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGccacca                                                                >  1:2459326/1‑100 (MQ=255)
                               aGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGccacca                                                                >  2:372559/1‑100 (MQ=255)
                               aGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGccacca                                                                >  1:3883093/1‑100 (MQ=255)
                                gCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCcaccac                                                               >  2:4612130/1‑100 (MQ=255)
                                gCCGATCTTGAGCGGTTGTTTCGCTTGTTGACGCAGCGATTTGCTTGTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCGAATCCGCGCCTGCcaccac                                                               >  1:4948749/1‑100 (MQ=255)
                                 ccGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACt                                                              >  1:2846118/1‑100 (MQ=255)
                                  cGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTc                                                             <  1:3213706/100‑1 (MQ=255)
                                  cGATCTTGAGCGGTTGTTTCGCTTGTGGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTc                                                             <  1:340105/100‑1 (MQ=255)
                                     tCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGAt                                                          <  2:206341/100‑1 (MQ=255)
                                      cTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGAtt                                                         >  2:3342010/1‑100 (MQ=255)
                                       ttGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTc                                                        <  1:4587367/100‑1 (MQ=255)
                                        tGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTcc                                                       <  2:4719157/100‑1 (MQ=255)
                                         gAGCGGTTGTTTCGCTTGTTGCCTCAGCGTTTTGCTTTTCTGCCTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCg                                                      <  2:3740861/100‑1 (MQ=255)
                                           gCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGc                                                    <  1:1179082/100‑1 (MQ=255)
                                             ggTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCAt                                                  >  2:1396550/1‑100 (MQ=255)
                                                tgttTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAGTCCGCGCCTGCCACCACTCGATTCCGGCATTCt                                               <  2:1541012/100‑1 (MQ=255)
                                                tgttTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCt                                               >  2:185336/1‑100 (MQ=255)
                                                tgttTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCt                                               >  1:4520614/1‑100 (MQ=255)
                                                tgttTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCt                                               <  2:4816155/100‑1 (MQ=255)
                                                   ttCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTgg                                            <  2:9033/100‑1 (MQ=255)
                                                   ttCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTgg                                            <  1:961955/100‑1 (MQ=255)
                                                   ttCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTgg                                            >  2:1967727/1‑100 (MQ=255)
                                                    tCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTggc                                           <  2:158830/100‑1 (MQ=255)
                                                       cTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTggcggc                                        >  2:854525/1‑100 (MQ=255)
                                                       cTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTggcggc                                        <  1:2742100/100‑1 (MQ=255)
                                                        ttgttgACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCt                                       <  1:2928688/100‑1 (MQ=255)
                                                           ttgACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACa                                    <  2:3510654/100‑1 (MQ=255)
                                                            tgACTCAGCGTTTTGCTTTTCTTACTCTTTGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACAt                                   >  2:4293820/1‑100 (MQ=255)
                                                             gACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACAtt                                  >  2:4797182/1‑100 (MQ=255)
                                                             gACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACAtt                                  >  1:3977534/1‑100 (MQ=255)
                                                             gACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACAtt                                  >  2:3618540/1‑100 (MQ=255)
                                                              aCTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATCCCGGCATTCTTGGCGGCTACATTg                                 >  2:446419/1‑100 (MQ=255)
                                                               ntCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTgc                                >  1:4443436/2‑100 (MQ=255)
                                                               cTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTgc                                <  2:5219619/100‑1 (MQ=255)
                                                               cTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTCCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTgc                                >  2:1673684/1‑100 (MQ=255)
                                                                 cAGCGTTTTGCTTTTCTGACTCAGGGCGGTGCCGCACCAGAAACGCCAAATCCGCGTCTGACACCACTCGATTAAGGCATTCTTGGCGGATACATTgcgc                              >  2:676877/1‑100 (MQ=255)
                                                                 cAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTgcgc                              >  1:288268/1‑100 (MQ=255)
                                                                   gCGTTTTGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGccc                            <  1:1992644/100‑1 (MQ=255)
                                                                        ttGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATaa                       <  1:4755712/100‑1 (MQ=255)
                                                                        ttGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATaa                       >  2:3046491/1‑100 (MQ=255)
                                                                        ttGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATaa                       >  1:2595701/1‑100 (MQ=255)
                                                                         tGCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAAt                      >  2:234263/1‑100 (MQ=255)
                                                                          gCTTTTCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAAtc                     >  2:306987/1‑100 (MQ=255)
                                                                              ttCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCAc                 <  1:1612311/100‑1 (MQ=255)
                                                                               tCTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCAcc                >  1:609116/1‑100 (MQ=255)
                                                                                cTGACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCa               <  1:375504/100‑1 (MQ=255)
                                                                                  gACTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATc             <  1:3773015/100‑1 (MQ=255)
                                                                                  gACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATc             <  1:2988213/100‑1 (MQ=255)
                                                                                    cTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCAc           >  2:4976058/1‑100 (MQ=255)
                                                                                    cTCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCAc           <  2:2926812/100‑1 (MQ=255)
                                                                                     tCAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACg          <  1:1502927/100‑1 (MQ=255)
                                                                                      cAGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCGTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGt         <  1:5241780/100‑1 (MQ=255)
                                                                                      cAGGGCGGAGACGCACCAGAAACGCAAAATCCGCGCCTGACACCAATCGATTCAGGCATTCTTGGCGGCTACATTGCGCCCGATAATATCACCATCAAGt         >  2:2755259/1‑100 (MQ=255)
                                                                                       aGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGtt        >  1:4030535/1‑100 (MQ=255)
                                                                                       aGGGCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGtt        <  2:5303006/100‑1 (MQ=255)
                                                                                         ggCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTAt      >  1:755007/1‑100 (MQ=255)
                                                                                          gCGGAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATc     <  1:1037023/100‑1 (MQ=255)
                                                                                            ggAGCCGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCgg   <  1:2205643/100‑1 (MQ=255)
                                                                                                 |                                                                                               
TGGTGGCGTTTGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGT  >  NC_000913/1083596‑1083788

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