Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,082,696 G→C 15.9% Q22E (CAG→GAG)  fdoH ← formate dehydrogenase‑O, Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,082,6960GC15.9% 651.2 / 101.8 233Q22E (CAG→GAG) fdoHformate dehydrogenase‑O, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base G (105/90);  new base C (7/30);  total (113/120)
Fisher's exact test for biased strand distribution p-value = 1.15e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

Only 100 of 239 total aligned reads displayed.

TGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAG  >  NC_000913/4082597‑4082795
                                                                                                   |                                                                                                   
tGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTg                                                                                                     >  1:3963932/1‑100 (MQ=255)
 gTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTg                                                                                                     <  2:614225/99‑1 (MQ=255)
  tCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAg                                                                                                   <  2:814343/100‑1 (MQ=255)
  tCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAg                                                                                                   >  1:525381/1‑100 (MQ=255)
   cGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAgg                                                                                                  <  1:2770901/100‑1 (MQ=255)
   cGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAgg                                                                                                  <  1:2288026/100‑1 (MQ=255)
    gTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGgc                                                                                                 >  1:3871997/1‑100 (MQ=255)
    gTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGgc                                                                                                 >  1:216764/1‑100 (MQ=255)
      tCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGgcgc                                                                                               <  2:991560/100‑1 (MQ=255)
        cACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCACCACTTCTTCCTGGAAGTCCCGCGCCTGAGgcg                                                                                                >  2:3740906/1‑97 (MQ=255)
         aCTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCggg                                                                                            <  1:1317982/100‑1 (MQ=255)
         aCTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCggg                                                                                            <  1:3526460/100‑1 (MQ=255)
          ctctGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCgggg                                                                                           <  1:1156558/100‑1 (MQ=255)
          ctctGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGGCGCggga                                                                                           <  2:1729173/100‑2 (MQ=255)
          ctctGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGGCGCggga                                                                                           <  2:119828/100‑2 (MQ=255)
           tctGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCggggg                                                                                          <  2:3687483/100‑1 (MQ=255)
            ctGAACACGCCACCTGACAGGCTTTACAGCCGATACATGTGGTAACGTCGATGAGTTTCGCCACGTCTTCCTGGAAGGCCCGCGCCTCAGGCGCggga                                                                                           <  1:3891202/98‑2 (MQ=25)
            ctGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGt                                                                                         <  1:2877151/100‑1 (MQ=255)
            ctGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGGCGCgggact                                                                                         <  1:1567909/100‑4 (MQ=255)
             tcaacgggCCACCTGACAGTCTTTACTGCCGATACAGGTGGTAACGTCGCTGAGTTTCGCCACTTCGCCCTGGAAGTGCCGCGCCTCAGGCGCggga                                                                                           <  2:3203679/90‑2 (MQ=255)
             tGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTc                                                                                        >  2:229589/1‑100 (MQ=255)
             tGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGGCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTc                                                                                        <  2:897168/100‑1 (MQ=255)
              gAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCa                                                                                       >  1:2403101/1‑100 (MQ=255)
               aaCACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAg                                                                                      >  2:3841018/1‑100 (MQ=255)
               aaCACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGCCCCGCGCCTCAGGCGCgggacttcc                                                                                      <  1:1680210/100‑7 (MQ=21)
                ccacGCCACCTGACAGGCTTTACAGCTGATACAGGTGGTAACGTCGATGAATTTAGCCACTTCTTCCAGCAAGTCCCGCGCCTGAGGCGAGGGGGTCAGa                                                                                     >  2:159484/2‑100 (MQ=21)
                acacGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGa                                                                                     >  1:2338803/1‑100 (MQ=255)
                acacGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGa                                                                                     <  1:2664681/100‑1 (MQ=255)
                 cacGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGAc                                                                                    >  1:543690/1‑100 (MQ=255)
                 cacGCCACCTGACAGGCTTTACAGCCAATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGAc                                                                                    >  1:2553784/1‑100 (MQ=255)
                  acGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGGCGCgggacttccagg                                                                                   <  1:4102200/100‑10 (MQ=25)
                    gCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGAcn                                                                                   <  2:2051415/98‑2 (MQ=255)
                    gCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGt                                                                                 <  1:1776841/100‑1 (MQ=255)
                     ccACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGtt                                                                                <  2:1242041/100‑1 (MQ=255)
                     ccACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGtt                                                                                >  2:3138330/1‑100 (MQ=255)
                     ccACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGtt                                                                                <  2:667234/100‑1 (MQ=255)
                     ccACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGtt                                                                                <  1:1482908/100‑1 (MQ=255)
                      cACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTa                                                                               <  1:2249119/100‑1 (MQ=255)
                      cACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGACGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTa                                                                               >  2:2243631/1‑100 (MQ=255)
                            aCAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTcgcn                                                                         <  1:1239261/100‑2 (MQ=255)
                            aCAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTcgcg                                                                         <  1:1898843/100‑1 (MQ=255)
                               ggCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGaa                                                                      >  1:3425712/1‑100 (MQ=255)
                                gCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAgatcggaaga                                                                     >  1:60897/1‑91 (MQ=255)
                                gCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAac                                                                     >  1:2961077/1‑100 (MQ=255)
                                gCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAac                                                                     >  2:3475158/1‑100 (MQ=255)
                                     aCAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAAcgacga                                                                <  1:1926175/100‑1 (MQ=255)
                                      cAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGaa                                                               <  2:233855/100‑1 (MQ=255)
                                        gCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATg                                                             >  1:2626739/1‑100 (MQ=255)
                                             gaCAGGTGGTAACGTCGATGAGGTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTAGCGGAACGACGAATGATATc                                                        <  2:2564759/99‑1 (MQ=255)
                                               cAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCtt                                                      >  2:1356271/1‑100 (MQ=255)
                                               cAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCtt                                                      >  1:717797/1‑100 (MQ=255)
                                                 ggtggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTga                                                    >  1:2698410/1‑100 (MQ=255)
                                                 ggtggtAACGTCGATGAGTTGCGCCACTGCTGCCTGGAAGTCCCGCGCCTGAGGGGCGGGGGTCAGAGCGTTAGTCGCGGAACGACGAATGATATCTTg                                                     <  2:1781098/99‑1 (MQ=255)
                                                      tAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTg                                               <  2:380269/100‑1 (MQ=255)
                                                      tAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTg                                               >  1:2213532/1‑100 (MQ=255)
                                                       aaCGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGa                                              >  2:2780122/1‑100 (MQ=255)
                                                         cGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATa                                            <  1:3026372/100‑1 (MQ=255)
                                                           tCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAg                                          <  1:2813716/100‑1 (MQ=255)
                                                           tCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAg                                          <  1:3716684/100‑1 (MQ=255)
                                                            cGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc                                         >  2:2844640/1‑100 (MQ=255)
                                                             gATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGACCGACGAATGATATCTTGAGATTGATAAGcc                                        <  2:3990658/100‑1 (MQ=255)
                                                              aTGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGTAACGACGAATGATATCTTGAGATTGATAAGCCa                                       <  2:2729275/100‑1 (MQ=255)
                                                               tGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCat                                      >  2:111490/1‑100 (MQ=255)
                                                                  gTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATg                                   >  1:55639/1‑100 (MQ=255)
                                                                  gTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATg                                   >  2:3020264/1‑100 (MQ=255)
                                                                   tttCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGATCAGACCGTTAGTCGCGGAAGGACGAATGATATCTTGAGATTGATAAGCCata                                     >  1:1821760/1‑97 (MQ=255)
                                                                     tCGCCACTTCTTCCTGGAATTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTCGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATgtcg                                >  1:1984025/1‑100 (MQ=255)
                                                                      cGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATgtcgt                               >  2:4152288/1‑100 (MQ=255)
                                                                       gCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGAGTGATAAGCCATATgtagtc                              >  1:3120796/1‑100 (MQ=255)
                                                                         cACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGtctc                            >  2:2356147/1‑100 (MQ=255)
                                                                          aCTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTcc                           <  1:3780862/100‑1 (MQ=255)
                                                                          aCTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTcc                           <  1:2237014/100‑1 (MQ=255)
                                                                             tcttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGtt                        <  2:2630018/100‑1 (MQ=255)
                                                                              cttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAACCATATGTCGTCTCCGtt                        >  1:413536/1‑99 (MQ=255)
                                                                                 ccTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTacac                    >  1:2571453/1‑100 (MQ=255)
                                                                                    ggAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCtt                 >  1:3322562/1‑100 (MQ=255)
                                                                                    ggAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCtt                 >  2:3913243/1‑100 (MQ=255)
                                                                                     gAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATGTCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCttt                >  2:2936244/1‑100 (MQ=255)
                                                                                      gacccccGCGCCTCAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCtttt               >  1:985172/5‑100 (MQ=255)
                                                                                      aaGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCtttt               >  1:2082421/1‑100 (MQ=255)
                                                                                        gTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTcc             >  2:1886207/1‑100 (MQ=255)
                                                                                        gTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTcc             <  2:444302/100‑1 (MQ=255)
                                                                                        cccccGCGCCTCAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTcc             >  2:3421856/3‑100 (MQ=255)
                                                                                         tCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCca            <  2:4086102/100‑1 (MQ=255)
                                                                                         ccccGCGCCTCAGGCGCGGGGGTCAGACCGGTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTCACACCTTTTCca            >  1:1615341/2‑100 (MQ=255)
                                                                                          cccGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCcac           >  2:3546983/1‑100 (MQ=255)
                                                                                          cccGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCca            <  2:543689/99‑1 (MQ=255)
                                                                                           ccGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCcaca          >  1:2735993/1‑100 (MQ=255)
                                                                                              cgcCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTc       <  1:916913/100‑1 (MQ=255)
                                                                                              cgcCTGAGGCGCGGGGGTCAGACCGTCAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTc       <  1:1278895/100‑1 (MQ=255)
                                                                                                ccTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTcac     <  1:2977966/100‑1 (MQ=255)
                                                                                                 cTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATTATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTcaca    <  1:1015359/100‑1 (MQ=255)
                                                                                                 cTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCCTATGTCGGCTCCGTTACACCTTTTCCACATTcaca    >  2:3860302/1‑100 (MQ=255)
                                                                                                 cTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTcaca    >  2:2170699/1‑100 (MQ=255)
                                                                                                 cTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTcaca    >  1:363679/1‑100 (MQ=255)
                                                                                                  tGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTcaca    >  2:1577391/1‑99 (MQ=255)
                                                                                                  tGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaa   >  1:1643489/1‑100 (MQ=255)
                                                                                                   gAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaag  <  2:3908656/100‑1 (MQ=255)
                                                                                                   gAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaag  >  1:3641207/1‑100 (MQ=255)
                                                                                                   gAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaag  <  1:1052090/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAG  >  NC_000913/4082597‑4082795

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