Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 429,838 A→G 5.5% S112G (AGC→GGC)  secF → Sec translocon accessory complex subunit SecF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913429,8380AG 5.5% 225.1 / 2.2 73S112G (AGC→GGC) secFSec translocon accessory complex subunit SecF
Reads supporting (aligned to +/- strand):  ref base A (33/36);  new base G (2/2);  total (35/38)
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.86e-01

CGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTAC  >  NC_000913/429690‑429984
                                                                                                                                                    |                                                                                                                                                  
cGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAg                                                                                                                                                   >  1:335840/1‑150 (MQ=255)
cGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAg                                                                                                                                                   >  2:392625/1‑150 (MQ=255)
     aaCCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAgg                                                                                                                                              <  1:722792/150‑1 (MQ=255)
      aCCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGt                                                                                                                                             <  2:135462/150‑1 (MQ=255)
             taaaTTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAg                                                                                                                                      <  1:135870/149‑1 (MQ=255)
              aaaTTGACTTAATGCGTGATTCATTTCAAAAAGCCTGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAgg                                                                                                                                     <  1:1002763/150‑1 (MQ=255)
               aaTTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGt                                                                                                                                    >  2:1234778/1‑150 (MQ=255)
                aTTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTg                                                                                                                                   <  2:789239/150‑1 (MQ=255)
                   gACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGAtt                                                                                                                                <  1:1015370/150‑1 (MQ=255)
                   gACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGAtt                                                                                                                                <  2:187694/150‑1 (MQ=255)
                    aCGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTa                                                                                                                               <  1:1385281/150‑1 (MQ=255)
                     cGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTaa                                                                                                                              >  1:615794/1‑150 (MQ=255)
                      gTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAAc                                                                                                                             <  1:250193/150‑1 (MQ=255)
                         aTGCGTTATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGGTGGGCAGCCAGGTTCTGAAGGTGATTAACGaa                                                                                                                          <  2:1297736/150‑1 (MQ=255)
                          tGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTc                                                                                                                                           >  2:1414563/1‑132 (MQ=255)
                          tGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTc                                                                                                                                           <  1:1414563/132‑1 (MQ=255)
                          tGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGACGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAAt                                                                                                                         >  2:1369066/1‑150 (MQ=255)
                            cGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATcc                                                                                                                       <  2:75364/150‑1 (MQ=255)
                             gTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATcca                                                                                                                      >  1:911890/1‑150 (MQ=255)
                             gTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATcca                                                                                                                      <  1:949087/150‑1 (MQ=255)
                              tGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATccac                                                                                                                     >  1:1283575/1‑150 (MQ=255)
                              tGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGACGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATccac                                                                                                                     <  1:1369066/150‑1 (MQ=255)
                               gATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATccacc                                                                                                                    <  1:976721/150‑1 (MQ=255)
                                 tGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCaa                                                                                                                  >  1:429808/1‑150 (MQ=255)
                                        cAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAAt                                                                                                           >  1:360055/1‑150 (MQ=255)
                                          aaaaGCTGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATgc                                                                                                         <  2:628647/150‑1 (MQ=255)
                                            aaGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATgcag                                                                                                       >  1:1318820/1‑150 (MQ=255)
                                             aGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATgcagc                                                                                                      >  2:963310/1‑150 (MQ=255)
                                               ccGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGcagcag                                                                                                    <  2:1459662/150‑1 (MQ=255)
                                                  gTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGa                                                                                                 <  1:816943/150‑1 (MQ=255)
                                                      tGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCg                                                                                             <  1:1211338/150‑1 (MQ=255)
                                                        aaGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTa                                                                                           >  2:1215435/1‑150 (MQ=255)
                                                           agCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTg                                                                                        <  2:314148/150‑1 (MQ=255)
                                                           agCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTg                                                                                        <  2:911890/150‑1 (MQ=255)
                                                           agCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTg                                                                                        <  2:360055/150‑1 (MQ=255)
                                                                 tgctgcAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTtcg                                                                                  <  1:1091750/150‑1 (MQ=255)
                                                                   ctgcAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTtcgtc                                                                                <  2:567152/150‑1 (MQ=255)
                                                                    tgcAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTtcgtcg                                                                               >  1:144965/1‑150 (MQ=255)
                                                                     gcAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGTTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCgg                                                                              <  2:157428/150‑1 (MQ=255)
                                                                     gcAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCgg                                                                              <  2:1283575/150‑1 (MQ=255)
                                                                     gcAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCgg                                                                              >  1:364743/1‑150 (MQ=255)
                                                                        aaaCTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTcc                                                                           <  2:335840/150‑1 (MQ=255)
                                                                                gTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTgg                                                                   >  1:223222/1‑150 (MQ=255)
                                                                                 tAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTggg                                                                  >  1:117365/1‑150 (MQ=255)
                                                                                  agcagcCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTgggg                                                                 <  2:719933/150‑1 (MQ=255)
                                                                                   gcagcCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTggggg                                                                >  1:216667/1‑150 (MQ=255)
                                                                                      gcCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAg                                                             <  1:963310/150‑1 (MQ=255)
                                                                                             catcatGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTgc                                                      >  2:1171260/1‑150 (MQ=255)
                                                                                              atcatGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTgcg                                                     >  1:47720/1‑150 (MQ=255)
                                                                                                            gccgccTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGAt                                       >  1:1278485/1‑150 (MQ=255)
                                                                                                               gccTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGc                                    >  1:174242/1‑150 (MQ=255)
                                                                                                                 ctgctgAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGt                                  <  2:223222/150‑1 (MQ=255)
                                                                                                                  tgctgAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGtt                                 >  1:1352759/1‑150 (MQ=255)
                                                                                                                     tgAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGAt                              >  1:95742/1‑150 (MQ=255)
                                                                                                                      gAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATg                             >  2:156280/1‑150 (MQ=255)
                                                                                                                        aGGCGAAACCGGCGGTCAGGTGTTGGGCGGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCt                                                         >  1:1294087/1‑120 (MQ=255)
                                                                                                                        aGGCGAAACCGGCGGTCAGGTGTTGGGCGGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCt                                                         <  2:1299388/120‑1 (MQ=255)
                                                                                                                        aGGCGAAACCGGCGGTCAGGTGTTGGGCGGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCt                                                         >  1:1299388/1‑120 (MQ=255)
                                                                                                                        aGGCGAAACCGGCGGTCAGGTGTTGGGCGGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCt                                                         <  2:1294087/120‑1 (MQ=255)
                                                                                                                        aGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGTAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTg                                           <  1:991642/134‑1 (MQ=255)
                                                                                                                        aGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGTAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTg                                           >  2:991642/1‑134 (MQ=255)
                                                                                                                        aGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGc                           <  1:1462630/150‑1 (MQ=255)
                                                                                                                         ggCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCa                          >  2:219962/1‑150 (MQ=255)
                                                                                                                           cGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAgc                        >  1:123019/1‑150 (MQ=255)
                                                                                                                               aCCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCgctg                    <  1:156280/150‑1 (MQ=255)
                                                                                                                                    cggTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTc               <  2:615794/150‑1 (MQ=255)
                                                                                                                                    cggTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTc               <  1:417585/150‑1 (MQ=255)
                                                                                                                                            tgtTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCg       >  1:765223/1‑150 (MQ=255)
                                                                                                                                             gtTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCgt      >  1:1312103/1‑150 (MQ=255)
                                                                                                                                              ttGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCgtg     >  2:640685/1‑150 (MQ=255)
                                                                                                                                               tGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCgtgt    >  2:372283/1‑150 (MQ=255)
                                                                                                                                               tGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCgtgt    <  2:765223/150‑1 (MQ=255)
                                                                                                                                                 ggCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTAc  >  1:1037041/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGGCAGCCAGGTTCTGAAGGTGATTAACGAATCCACCAATCAGAATGCAGCAGTGAAGCGTATTGAGTTCGTCGGTCCGAGCGTGGGGGCAGACCTTGCGCAAACCGGTGCGATGGCGTTGATGGCAGCGCTGCTGTCTATCCTCGTGTAC  >  NC_000913/429690‑429984

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