Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,375,484 A→T 15.7% L160H (CTT→CAT) ‡ menC ← O‑succinylbenzoyl‑CoA synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,375,4840AT15.7% 552.2 / 84.4 217L160H (CTT→CAT) menCO‑succinylbenzoyl‑CoA synthase
Reads supporting (aligned to +/- strand):  ref base A (84/98);  new base T (11/23);  total (96/121)
Fisher's exact test for biased strand distribution p-value = 1.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

Only 100 of 229 total aligned reads displayed.

CACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCT  >  NC_000913/2375385‑2375583
                                                                                                   |                                                                                                   
cACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCCAACTGCTGACCTTTAAGCGGTGTCCCGGCCCGATTTGCGTCaa                                                                                                     >  2:3013405/1‑100 (MQ=255)
cACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCaa                                                                                                     >  1:1162702/1‑100 (MQ=255)
cACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCaa                                                                                                     >  1:2845620/1‑100 (MQ=255)
cACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCaa                                                                                                     >  1:2290139/1‑100 (MQ=255)
 aCGGCTCTTCGAGAAGCGCGATGCGGTCGCGAGAAGCCGGGTTAACTTATTAGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCaa                                                                                                     <  1:1306264/99‑1 (MQ=255)
  cGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGa                                                                                                   <  2:1287325/100‑1 (MQ=255)
   ggCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGAc                                                                                                  >  1:3400577/1‑100 (MQ=255)
      tctTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCa                                                                                               <  1:43470/100‑1 (MQ=255)
      tctTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCa                                                                                               <  1:4012993/100‑1 (MQ=255)
       ctTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCaa                                                                                              >  2:1650916/1‑100 (MQ=255)
         tCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAAt                                                                                            >  2:2641411/1‑100 (MQ=255)
         tCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAAt                                                                                            <  2:2734553/100‑1 (MQ=255)
          cGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCTTGACGCAAAtc                                                                                           <  2:918636/100‑2 (MQ=255)
          cGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAAt                                                                                            <  1:1406375/99‑1 (MQ=255)
           gagaAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCTTGACGCAAAtcg                                                                                          <  1:2012743/100‑3 (MQ=255)
           gagaAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCTTGACGCAAAtcg                                                                                          <  2:826671/100‑3 (MQ=255)
           gagaAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCTTGACGCAAAtcg                                                                                          <  2:1094206/100‑3 (MQ=255)
           gagaAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCTTGACGCAAAtcg                                                                                          <  1:1968511/100‑3 (MQ=255)
           gagaAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCTTGACGCAAAtcg                                                                                          <  1:545347/100‑3 (MQ=255)
           gagaAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCTTGACGCAAAtcg                                                                                          <  2:4186134/100‑3 (MQ=255)
           gagaAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGc                                                                                          <  1:1091489/100‑1 (MQ=255)
           gagaAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGc                                                                                          >  2:2067821/1‑100 (MQ=255)
            agaAACGCGATGCGGTCGCGATAATCCGGTTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCTTGACGCAAAtcgc                                                                                         <  1:1015456/100‑4 (MQ=255)
            agaAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCTTGACGCAAAtcgc                                                                                         <  2:4285244/100‑4 (MQ=255)
            agaAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCTTGACGCAAAtcgc                                                                                         <  1:2276943/100‑4 (MQ=255)
             gaAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCTTGACGCAAAtcgcg                                                                                        <  2:2300368/100‑5 (MQ=255)
             gaAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCTTGACGCAAAtcgcg                                                                                        <  2:1949565/100‑5 (MQ=255)
             gaAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAg                                                                                        <  1:2629932/100‑1 (MQ=255)
              aaaCGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCTTGACGCAAAtcgcgc                                                                                       <  1:1895264/100‑6 (MQ=255)
                aCGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATc                                                                                     <  2:4021434/100‑1 (MQ=255)
                aCGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACGGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCTTGACGCAAAtcgcgcc                                                                                      <  1:2877411/99‑7 (MQ=38)
                  gcgATGCGGTCGCGATAATCCGGGTTACCGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCg                                                                                   <  1:3602666/100‑1 (MQ=255)
                  gcgATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGTGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCg                                                                                   <  1:724153/100‑1 (MQ=255)
                  gcgATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCg                                                                                   >  1:265232/1‑100 (MQ=255)
                  gcgATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCg                                                                                   >  1:240001/1‑100 (MQ=255)
                  gcgATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCg                                                                                   <  1:652990/100‑1 (MQ=255)
                     aTGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGaa                                                                                <  1:2358566/100‑1 (MQ=255)
                     aTGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGaa                                                                                <  2:2228630/100‑1 (MQ=255)
                     aTGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGaa                                                                                <  1:1559588/100‑1 (MQ=255)
                       gCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAAtt                                                                              <  2:307302/100‑1 (MQ=255)
                        cGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTg                                                                             <  1:3860529/100‑1 (MQ=255)
                             gcgATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTcc                                                                        <  1:92001/100‑1 (MQ=255)
                               gATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAg                                                                      >  2:3254458/1‑100 (MQ=255)
                                aTAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGc                                                                     <  2:686929/100‑1 (MQ=255)
                                 tAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCa                                                                    <  2:2163850/100‑1 (MQ=255)
                                 tAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCa                                                                    >  1:888947/1‑100 (MQ=255)
                                  aaTCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATGTGCGTAAAGACGCAAATGCAGATCCGGAATTGCCTCGAGCaa                                                                   >  1:271075/1‑100 (MQ=255)
                                   aTCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAAc                                                                  <  2:2754053/100‑1 (MQ=255)
                                    tCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACa                                                                 <  2:841868/100‑1 (MQ=255)
                                     ccGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACa                                                                 <  1:2054862/99‑1 (MQ=255)
                                            aaCGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGTAACAGATTcacc                                                         >  1:3589774/1‑100 (MQ=255)
                                            aaCGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTcacc                                                         <  1:1644941/100‑1 (MQ=255)
                                              cGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTcaccac                                                       >  2:3631941/1‑100 (MQ=255)
                                              cGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTcaccac                                                       >  2:1624071/1‑100 (MQ=255)
                                              cGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTcaccac                                                       >  2:4111529/1‑100 (MQ=255)
                                                   ttGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGc                                                  >  1:1429233/1‑100 (MQ=255)
                                                   ttGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGc                                                  <  1:3254441/100‑1 (MQ=255)
                                                   ttGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGc                                                  >  2:384715/1‑100 (MQ=255)
                                                    tGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGcc                                                 <  1:3695928/100‑1 (MQ=255)
                                                      gCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGt                                               >  1:2246547/1‑100 (MQ=255)
                                                       cAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTc                                              >  2:2081587/1‑100 (MQ=255)
                                                       cAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTc                                              <  2:3569550/100‑1 (MQ=255)
                                                            tgctgACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCa                                         >  2:3277336/1‑100 (MQ=255)
                                                              ctgACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCAcc                                       >  2:2790275/1‑100 (MQ=255)
                                                               tgACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCg                                      <  2:1395579/100‑1 (MQ=255)
                                                                gACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGc                                     >  2:1540324/1‑100 (MQ=255)
                                                                  ccTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCCGCTCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCtt                                   >  1:760630/1‑100 (MQ=255)
                                                                   cTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTc                                  >  1:4139236/1‑100 (MQ=255)
                                                                   cTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTc                                  >  2:2477082/1‑100 (MQ=255)
                                                                    tttCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCg                                 <  1:2512634/100‑1 (MQ=255)
                                                                     ttCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGt                                <  2:3170119/100‑1 (MQ=255)
                                                                     ttCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGt                                <  2:1249759/100‑1 (MQ=255)
                                                                      tCAGCGGTGTCCAGGCGCGATTTGCGTCAAGCCGTATATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGt                                <  1:1761352/99‑1 (MQ=255)
                                                                       cAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTAc                              <  1:335003/100‑1 (MQ=255)
                                                                       cAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTAc                              <  2:3906179/100‑1 (MQ=255)
                                                                        aGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACa                             <  2:2465341/100‑1 (MQ=255)
                                                                         gCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACaa                            >  2:18290/1‑100 (MQ=255)
                                                                         gCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACaa                            >  2:593934/1‑100 (MQ=255)
                                                                             tgtCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATccc                        <  1:2199081/100‑1 (MQ=255)
                                                                               tCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCAc                      <  1:941234/100‑1 (MQ=255)
                                                                                ccAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACt                     >  2:3488156/1‑100 (MQ=255)
                                                                                 cAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACtt                    >  1:2580665/1‑100 (MQ=255)
                                                                                 cAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACtt                    >  2:358030/1‑100 (MQ=255)
                                                                                 cAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACtt                    >  1:419075/1‑100 (MQ=255)
                                                                                  aGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACttt                   >  2:2595264/1‑100 (MQ=255)
                                                                                   ggCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTg                  >  2:1731284/1‑100 (MQ=255)
                                                                                    gcgcGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGa                 >  2:17668/1‑100 (MQ=255)
                                                                                    gcgcGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGa                 <  2:3039955/100‑1 (MQ=255)
                                                                                     cgcgATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGAc                <  1:1075937/100‑1 (MQ=255)
                                                                                      tgcaTTTGCGTCTTGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGAcc               >  2:1680259/4‑100 (MQ=255)
                                                                                      gcgATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGAcc               <  2:1163073/100‑1 (MQ=255)
                                                                                       cgATTTGCGTCAAGACGCAAATGCAGATCCGGAATTACCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCt              >  2:3103125/1‑100 (MQ=255)
                                                                                         aTTTGCGTCTTGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTc            >  1:2964117/1‑100 (MQ=255)
                                                                                         aTTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTc            >  2:4245900/1‑100 (MQ=255)
                                                                                            tGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGTCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGcc         >  2:2055475/1‑100 (MQ=255)
                                                                                            tGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGcc         >  1:847924/1‑100 (MQ=255)
                                                                                               gTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACt      >  2:1390723/1‑100 (MQ=255)
                                                                                               gTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACt      <  2:753789/100‑1 (MQ=255)
                                                                                                tCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACtt     <  2:1558356/100‑1 (MQ=255)
                                                                                                   aGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTtct  >  1:3308258/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCT  >  NC_000913/2375385‑2375583

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: