Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,466,858 G→C 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,8580GC97.0% 353.7 / 2.0 101A45G (GCC→GGC) bfdbacterioferritin‑associated ferredoxin
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base C (55/43);  total (56/45)
Fisher's exact test for biased strand distribution p-value = 5.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

Only 100 of 103 total aligned reads displayed.

AACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTT  >  NC_000913/3466762‑3466954
                                                                                                |                                                                                                
aaCGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGcccc                                                                                               >  1:1139265/1‑100 (MQ=255)
  cGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCgg                                                                                             >  1:1839260/1‑100 (MQ=255)
  cGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCgg                                                                                             >  1:714208/1‑100 (MQ=255)
     cTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGac                                                                                           >  2:867341/1‑99 (MQ=255)
       tACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGAcaca                                                                                        >  1:2236268/1‑100 (MQ=255)
        aCGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACAt                                                                                       <  2:1139263/100‑1 (MQ=255)
         cGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACAtt                                                                                      <  1:281972/100‑1 (MQ=255)
          gAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACAttt                                                                                     <  2:1990872/100‑1 (MQ=255)
          gAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACAttt                                                                                     >  1:1964961/1‑100 (MQ=255)
           aGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTa                                                                                    >  1:1746012/1‑100 (MQ=255)
           aGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTa                                                                                    >  2:563795/1‑100 (MQ=255)
           aGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTa                                                                                    >  1:716913/1‑100 (MQ=255)
            gTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTAc                                                                                   >  1:384303/1‑100 (MQ=255)
             tAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTAcc                                                                                  >  2:2237480/1‑100 (MQ=255)
                gTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCa                                                                               >  2:2301349/1‑100 (MQ=255)
                 tCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCAc                                                                              >  2:2266195/1‑100 (MQ=255)
                 tCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCAc                                                                              >  1:2022234/1‑100 (MQ=255)
                 tCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCAc                                                                              <  1:1349085/100‑1 (MQ=255)
                  cAAAAAACACACTCCGTTTATGCGGCCTCCTTAAACTCCAGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACt                                                                             >  2:2092192/1‑100 (MQ=255)
                  cAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACt                                                                             >  1:594282/1‑100 (MQ=255)
                  cAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACt                                                                             <  2:673525/100‑1 (MQ=255)
                   aaaaaaCACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTg                                                                            >  2:478681/1‑100 (MQ=255)
                   aaaaaaCACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTg                                                                            >  1:1968869/1‑100 (MQ=255)
                       aaCACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTg                                                                        <  2:1864628/100‑1 (MQ=255)
                          acacTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCg                                                                     <  2:1964954/100‑1 (MQ=255)
                           cacTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGa                                                                    >  1:2071157/1‑100 (MQ=255)
                             cTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGAcc                                                                  >  2:1919887/1‑100 (MQ=255)
                               ccGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCgg                                                                <  1:412132/100‑1 (MQ=255)
                               ccGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCgg                                                                <  2:2071150/100‑1 (MQ=255)
                               ccGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCgg                                                                >  1:1528472/1‑100 (MQ=255)
                                cGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGa                                                               >  2:13920/1‑100 (MQ=255)
                                    tATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATaa                                                           >  1:1500311/1‑100 (MQ=255)
                                      tGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAAt                                                         <  1:1587054/100‑1 (MQ=255)
                                       gCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAAtt                                                        >  1:1734581/1‑100 (MQ=255)
                                         ggACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAAtttt                                                      <  2:1963840/100‑1 (MQ=255)
                                         ggACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAAtttt                                                      <  2:438252/100‑1 (MQ=255)
                                         ggACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAAtttt                                                      <  1:2335651/100‑1 (MQ=255)
                                          gACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAAttttt                                                     <  2:384301/100‑1 (MQ=255)
                                          gACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAAttttt                                                     <  2:2236260/100‑1 (MQ=255)
                                           aCTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTCCCGCACTGATTGCCGACCGGAATAAAtttttt                                                    <  1:2218284/100‑1 (MQ=255)
                                           aCTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAAtttttt                                                    <  1:1378909/100‑1 (MQ=255)
                                            cTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAAttttttt                                                   >  2:301799/1‑100 (MQ=255)
                                              ccTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTaa                                                 >  2:2258204/1‑100 (MQ=255)
                                              ccTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTaa                                                 >  2:913134/1‑100 (MQ=255)
                                              ccTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTaa                                                 >  1:2052782/1‑100 (MQ=255)
                                              ccTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTCCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTaa                                                 >  2:2185401/1‑100 (MQ=255)
                                               cTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAAt                                                >  2:744927/1‑100 (MQ=255)
                                               cTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAAt                                                <  2:642354/100‑1 (MQ=255)
                                               cTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAAt                                                <  2:1839255/100‑1 (MQ=255)
                                                ttAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAAtt                                               >  1:896682/1‑100 (MQ=255)
                                                 tAAACTCCGGCAGCTGCATTAATTTATCCTCCATCACTTCACGCGAGCCCCGGACACCTTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAAtt                                               >  2:1860270/1‑99 (MQ=255)
                                                  aaaCTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGc                                             >  2:27441/1‑100 (MQ=255)
                                                  aaaCTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGc                                             >  2:1870885/1‑100 (MQ=255)
                                                   aaCTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCt                                            >  1:1200086/1‑100 (MQ=255)
                                                    aCTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCt                                            >  2:930240/1‑99 (MQ=255)
                                                    aCTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTg                                           >  2:1114723/1‑100 (MQ=255)
                                                     cTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCGCGCGCGCACCGGACACATTTACACCAGTGATTGCCGACCGGTATAAATTTTTTTAATTGCTg                                           <  2:440606/99‑1 (MQ=255)
                                                     cTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTgg                                          <  1:1710010/100‑1 (MQ=255)
                                                     cTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTgg                                          <  1:2266203/100‑1 (MQ=255)
                                                       ccGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACAGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGaa                                        >  2:1641187/1‑100 (MQ=255)
                                                       ccGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACAGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGaa                                        >  2:1641695/1‑100 (MQ=255)
                                                        cGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAAc                                       <  1:1954126/100‑1 (MQ=255)
                                                             gCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTg                                  <  1:478683/100‑1 (MQ=255)
                                                             gCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTg                                  >  1:402831/1‑100 (MQ=255)
                                                             gCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTg                                  >  1:1250670/1‑100 (MQ=255)
                                                              cTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGa                                 <  2:714207/100‑1 (MQ=255)
                                                              cTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCACACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGa                                 <  2:313602/100‑1 (MQ=255)
                                                               tGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAg                                >  1:208067/1‑100 (MQ=255)
                                                                ctaTTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAgg                               <  1:1375531/98‑1 (MQ=255)
                                                                 cATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAggg                              >  1:1912534/1‑100 (MQ=255)
                                                                  aTTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAgggg                             <  1:13921/100‑1 (MQ=255)
                                                                  aTTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAgggg                             <  2:1200083/100‑1 (MQ=255)
                                                                  aTTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAggag                             >  2:1375528/1‑98 (MQ=255)
                                                                    tAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGaa                           >  2:336219/1‑100 (MQ=255)
                                                                     aaTTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGaaa                          >  1:1779615/1‑100 (MQ=255)
                                                                       ttCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAAc                        <  1:1919894/100‑1 (MQ=255)
                                                                         cATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTg                      >  2:2122006/1‑100 (MQ=255)
                                                                          aTCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGAGGAAAACTgg                     <  2:1647373/100‑1 (MQ=255)
                                                                           tcctccATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGgc                    >  1:641780/1‑100 (MQ=255)
                                                                            cctccATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTg                      >  2:819230/1‑97 (MQ=255)
                                                                            cctccATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGgcg                   <  1:1870891/100‑1 (MQ=255)
                                                                            cctccATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGgcg                   >  2:595580/1‑100 (MQ=255)
                                                                             ctccATCACTTCACGGGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGgcgc                  <  1:2301358/100‑1 (MQ=255)
                                                                             ctccATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGgcgc                  <  1:563796/100‑1 (MQ=255)
                                                                             ctccATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGgcgc                  <  2:2052775/100‑1 (MQ=255)
                                                                             ctccATCACTTCACGCACGCCCCGTACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGgcgc                  <  2:2331180/100‑1 (MQ=255)
                                                                              tccATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGca                 <  2:133974/100‑1 (MQ=255)
                                                                                cATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGcaca               >  1:468398/1‑100 (MQ=255)
                                                                                 aTCGCTTCACGCGCGCCCCGGACACTTTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAg              <  2:225306/100‑1 (MQ=255)
                                                                                 aTCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAg              <  1:2122013/100‑1 (MQ=255)
                                                                                  tCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGc             >  1:1226387/1‑100 (MQ=255)
                                                                                  tCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGc             <  1:913136/100‑1 (MQ=255)
                                                                                  tCACTCCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGc             <  1:2185408/100‑1 (MQ=255)
                                                                                   cACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAGAACTGGCGCACAGcc            >  1:2305678/1‑100 (MQ=255)
                                                                                   cACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGcc            >  2:706322/1‑100 (MQ=255)
                                                                                   cACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGcc            <  1:27442/100‑1 (MQ=255)
                                                                                    aCTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGccn           <  1:1996531/100‑2 (MQ=255)
                                                                                        cACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAACAAATTGTTTTAATTGCTGGAACGAGTGAGAGGAAAACTGGCGCACAGCCAGACg       >  2:750933/1‑100 (MQ=255)
                                                                                         aCGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGa      <  2:208065/100‑1 (MQ=255)
                                                                                             gcgcCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAAttt  <  2:1968862/100‑1 (MQ=255)
                                                                                                |                                                                                                
AACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTT  >  NC_000913/3466762‑3466954

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: ‑