Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,466,858 G→C 16.1% A45G (GCC→GGC)  bfd ← bacterioferritin‑associated ferredoxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,466,8580GC16.1% 521.7 / 76.1 190A45G (GCC→GGC) bfdbacterioferritin‑associated ferredoxin
Reads supporting (aligned to +/- strand):  ref base G (74/85);  new base C (15/16);  total (89/101)
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 199 total aligned reads displayed.

TAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTT  >  NC_000913/3466759‑3466957
                                                                                                   |                                                                                                   
 ggAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGc                                                                                                    >  1:1680478/2‑100 (MQ=255)
 aGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATTCTGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGc                                                                                                    <  2:2618704/100‑1 (MQ=255)
 aGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGc                                                                                                    >  2:2711666/1‑100 (MQ=255)
 aGAACGGCTTACGAGTAAGTCAAAAAACACACTCAGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCgg                                                                                                     >  2:1442405/1‑99 (MQ=255)
    aCGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCg                                                                                                 <  1:503015/100‑1 (MQ=255)
     cGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCgg                                                                                                <  1:2843091/100‑1 (MQ=255)
     cGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCgg                                                                                                <  1:3260284/100‑1 (MQ=255)
     cGGCGTACGAGTAAGTCAAAACACACACTCCGTTTATGCGGACTCCTTGAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCgg                                                                                                <  2:1671874/100‑1 (MQ=255)
      ggCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGa                                                                                               >  2:3162331/1‑100 (MQ=255)
       gCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAGCTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGac                                                                                              <  2:1097576/100‑1 (MQ=255)
       gCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGac                                                                                              <  1:1110362/100‑1 (MQ=255)
          tACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGAcaca                                                                                           >  2:2133927/1‑100 (MQ=255)
            cGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACAtt                                                                                         >  1:1373726/1‑100 (MQ=255)
            cGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACAtt                                                                                         >  1:1911965/1‑100 (MQ=255)
            cGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCCCGGCCCGGACACAtt                                                                                         >  2:3509885/1‑100 (MQ=255)
             gAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACAttt                                                                                        <  1:1406120/100‑1 (MQ=255)
             gAGTAAGTCAAAAAACACACTCCGTTTATGCCGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACAtt                                                                                         <  2:488672/99‑1 (MQ=255)
                  aGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTCATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGc                                                                                   >  1:3768912/1‑100 (MQ=255)
                   gTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCa                                                                                  >  1:980730/1‑100 (MQ=255)
                    tCAAAAAACACACTCGGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCAc                                                                                 >  2:975734/1‑100 (MQ=255)
                    tCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCAc                                                                                 >  2:928059/1‑100 (MQ=255)
                    tCAAAAAACACACTCCGTTTATGCGGACTCCTTAAAATCCGGCAGCTGCATTAATACATCCTCCATCACTTCAAGCGCGGCCCGGACACATTTCCCGCa                                                                                  >  2:3924443/1‑99 (MQ=255)
                     gaaaaaaCACACACCGTGTCTGCGGACCCATTAAACTCCGGCAGCCGCCTTAATTCAGCCTCCATCACCTCACGCGCGCCCCGGACACCTTTACCGCAc                                                                                 >  2:747644/2‑99 (MQ=255)
                      aaaaaaCACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTg                                                                               <  1:2521982/100‑1 (MQ=255)
                      aaaaaaCACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTg                                                                               >  1:3837781/1‑100 (MQ=255)
                      aaaaaaCACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTg                                                                               <  2:4111016/100‑1 (MQ=255)
                           acacacTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGc                                                                          >  1:3679154/1‑100 (MQ=255)
                            cacacTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGcc                                                                         <  2:233345/100‑1 (MQ=255)
                            cacacTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGcc                                                                         <  1:3208873/100‑1 (MQ=255)
                             acacTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCg                                                                        <  1:264625/100‑1 (MQ=255)
                               acTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGAc                                                                      <  2:73803/100‑1 (MQ=255)
                                cTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGAc                                                                      >  2:1185312/1‑99 (MQ=255)
                                 tCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCg                                                                    >  1:1057732/1‑100 (MQ=255)
                                   cGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGa                                                                  >  1:1350373/1‑100 (MQ=255)
                                     tttATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAAt                                                                <  1:792470/100‑1 (MQ=255)
                                       tATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATaa                                                              >  2:2051614/1‑100 (MQ=255)
                                       tATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATaa                                                              >  1:275883/1‑100 (MQ=255)
                                          gCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAAtt                                                           <  2:2961231/100‑1 (MQ=255)
                                          gCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAAtt                                                           <  2:125702/100‑1 (MQ=255)
                                           cGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAAttt                                                          >  1:3568162/1‑100 (MQ=255)
                                            ggACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAAtttt                                                         <  1:2336169/100‑1 (MQ=255)
                                              aCTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAAtttttt                                                       <  1:1470315/100‑1 (MQ=255)
                                              aCTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAAtttttt                                                       <  2:2451049/100‑1 (MQ=255)
                                              aCTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAAtttttt                                                       <  1:3092563/100‑1 (MQ=255)
                                               cTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAAttttttt                                                      <  2:653246/100‑1 (MQ=255)
                                                 ccTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTaa                                                    >  2:2197668/1‑100 (MQ=255)
                                                 ccTTAAACTCCGGCAGCTGCATTAATTCAACCTCCATCACTACACGCGCGGCCCGGACACATTCACCGCACTGATTGCCGACCGGAATAAATTTTTTTaa                                                    >  1:1194892/1‑100 (MQ=255)
                                                  cTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAAt                                                   >  2:1595733/1‑100 (MQ=255)
                                                  cTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAAt                                                   <  2:103449/100‑1 (MQ=255)
                                                    tAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTg                                                 >  1:1544064/1‑100 (MQ=255)
                                                     aaaCTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGc                                                <  1:562050/100‑1 (MQ=255)
                                                     aaaCTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGc                                                <  1:2167718/100‑1 (MQ=255)
                                                     aaaCTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGc                                                >  1:2311644/1‑100 (MQ=255)
                                                        cTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTgg                                             >  1:2320443/1‑100 (MQ=255)
                                                         tCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGa                                            >  1:1592253/1‑100 (MQ=255)
                                                         tCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGa                                            >  1:704417/1‑100 (MQ=255)
                                                               aGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGt                                      <  2:1723292/100‑1 (MQ=255)
                                                                 cTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTGTTTAATTGCTGGAACGAGTg                                     <  2:2537978/99‑1 (MQ=255)
                                                                  tGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAg                                   <  1:4106247/100‑1 (MQ=255)
                                                                  tGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAg                                   >  1:1222869/1‑100 (MQ=255)
                                                                  tGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAg                                   <  1:3871030/100‑1 (MQ=255)
                                                                   gCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAgg                                  >  1:3510139/1‑100 (MQ=255)
                                                                   gCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAg                                   >  2:2481227/1‑99 (MQ=255)
                                                                    cATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAggg                                 >  1:1972192/1‑100 (MQ=255)
                                                                    cATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAggg                                 <  2:3331459/100‑1 (MQ=255)
                                                                     aTTAATTCATCCGCCATCACTTCACGCGCGCCCCGGACACAGTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAgggg                                <  1:278155/100‑1 (MQ=255)
                                                                      ttAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGa                               <  2:3739976/100‑1 (MQ=255)
                                                                      ttAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGa                               >  1:2241690/1‑100 (MQ=255)
                                                                        aaTTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGaaa                             <  1:1178871/100‑1 (MQ=255)
                                                                         aTTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGaaaa                            >  1:653943/1‑100 (MQ=255)
                                                                           tCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACt                          >  2:1882718/1‑100 (MQ=255)
                                                                           tCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACt                          >  1:2103384/1‑100 (MQ=255)
                                                                            cATCCTCCATCACGTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTg                         >  2:2571777/1‑100 (MQ=255)
                                                                             aTCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACt                          >  2:3371707/1‑98 (MQ=255)
                                                                              tcctccATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGgc                       <  2:2807486/100‑1 (MQ=255)
                                                                              tcctccATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGgc                       <  1:1836745/100‑1 (MQ=255)
                                                                              tcctccATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGgc                       <  2:2109470/100‑1 (MQ=255)
                                                                               cctccATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGgcg                      <  1:4169711/100‑1 (MQ=255)
                                                                               cctccATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGgcg                      <  2:813411/100‑1 (MQ=255)
                                                                               cctccATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGgcg                      <  2:250977/100‑1 (MQ=255)
                                                                                 tccATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGca                    <  2:511197/100‑1 (MQ=255)
                                                                                   cATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGcaca                  >  1:3995646/1‑100 (MQ=255)
                                                                                     tCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAGTTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGc                <  1:2942586/100‑1 (MQ=255)
                                                                                      cACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGcc               <  1:3182673/100‑1 (MQ=255)
                                                                                       aCTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCt              <  1:1928141/100‑1 (MQ=255)
                                                                                         ttCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTgn            <  2:3053671/100‑2 (MQ=255)
                                                                                          tCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGAc           <  1:2988187/100‑1 (MQ=255)
                                                                                          tCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGAc           <  2:768177/100‑1 (MQ=255)
                                                                                           cACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACg          <  2:1902915/100‑1 (MQ=255)
                                                                                           cACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACg          >  2:3184545/1‑100 (MQ=255)
                                                                                           cACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACg          <  2:182763/100‑1 (MQ=255)
                                                                                               cgcgGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAAtt      >  1:646618/1‑100 (MQ=255)
                                                                                                gcgGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAAttt     <  1:1062528/100‑1 (MQ=255)
                                                                                                 cgGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAAtttt    <  2:2775606/100‑1 (MQ=255)
                                                                                                   gCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAAtttttt  <  1:45088/100‑1 (MQ=255)
                                                                                                   gCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAAtttttt  <  1:3410022/100‑1 (MQ=255)
                                                                                                   gCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAAtttttt  <  1:2593201/100‑1 (MQ=255)
                                                                                                   gCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAAtttttt  <  1:2565039/100‑1 (MQ=255)
                                                                                                   gCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAAtttttt  >  1:2398238/1‑100 (MQ=255)
                                                                                                   gCCCCGACACATTTACCGCACTGATTATCGTCCGGAATAAATTTTTTTTATTGCTGGAACTAGTGTGGCGAAAAATGGCGCACTGCCTGACGAAtgttt   >  2:3850284/1‑99 (MQ=255)
                                                                                                   |                                                                                                   
TAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTT  >  NC_000913/3466759‑3466957

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