Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,034,751 Δ1 bp 37.0% coding (1498/1794 nt) hyaB → hydrogenase 1, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,034,7500C.37.0% 381.1 / 282.7 184coding (1497/1794 nt)hyaBhydrogenase 1, large subunit
Reads supporting (aligned to +/- strand):  ref base C (50/66);  new base . (31/37);  total (81/103)
Fisher's exact test for biased strand distribution p-value = 7.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

Only 100 of 192 total aligned reads displayed.

GCAGTATTTCTTCGACAAGCTGATGACCAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGATCTCTACCAGTGCGTGGTGCCGACC  >  NC_000913/1034651‑1034849
                                                                                                   |                                                                                                   
gCAGTATTTCTTCGACAAGCTGATGACCAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGc                                                                                                     >  2:2711391/1‑100 (MQ=255)
 cAGTATTTCTTCGACAAGCTGATGACCAACCTGAAAAACGGCAATCTCACCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGcc                                                                                                    >  2:3193188/1‑100 (MQ=255)
    tATTTCTTCGACAAGCTGATGACCAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGcgtgg                                                                                                 <  2:3143657/100‑5 (MQ=255)
     aTTTCTTCGACAAGCTGATGACCAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCgtgg                                                                                                <  1:1762682/100‑1 (MQ=255)
      tttCTTCGACAAGCTGATGACCAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCgtggt                                                                                               <  1:1137811/100‑1 (MQ=255)
        tcttcGACAAGCTGATGACCAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGgtgt                                                                                             >  1:362334/1‑100 (MQ=255)
        tcttcGACAAGCTGATGACCAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGgtgt                                                                                             >  1:335313/1‑100 (MQ=255)
        tcttcGACAAGCTGATGACCAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCGTGGTGTc                                                                                            >  2:4223806/1‑100 (MQ=255)
         cttcGACAAGCTGATGACCAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCGTGGTGTCg                                                                                           >  2:202754/1‑100 (MQ=255)
         cttcGACAAGCTGATGACCAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCGTGGTGTCg                                                                                           <  2:2957558/100‑1 (MQ=255)
           tcGACAAGCTGATGACCAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCgg                                                                                          <  2:2197733/100‑1 (MQ=255)
             gACAAGCTGATGACCAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGtt                                                                                        <  2:4176348/100‑1 (MQ=255)
             gACAAGCTGATGACCAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCGTGGTGTCGGttt                                                                                       <  1:172804/100‑1 (MQ=255)
              aCAAGCTGATGACCAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGttt                                                                                       <  2:3624669/100‑1 (MQ=255)
                 aGCTGATGACCAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTa                                                                                     >  1:2719097/1‑99 (MQ=255)
                 aGCTGATGACCAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCGTGGTGTCGGTTTTAcc                                                                                   <  2:1496683/100‑1 (MQ=255)
                  gCTGATGACCAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTAcc                                                                                   <  2:2767490/100‑1 (MQ=255)
                   cTGATGACCAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCg                                                                                  >  1:447373/1‑100 (MQ=255)
                    tgatgaCCAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGa                                                                                 >  2:1748030/1‑100 (MQ=255)
                     gatgaCCAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGaa                                                                                <  2:3247536/100‑1 (MQ=255)
                     gatgaCCAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCGTGGTGTCGGTTTTACCGAAg                                                                               >  1:1821554/1‑100 (MQ=255)
                     gatgaCCAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAg                                                                               >  1:3281617/1‑100 (MQ=255)
                       tgaCCAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAgc                                                                              >  1:3648560/1‑100 (MQ=255)
                       tgaCCAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCGTGGTGTCGGTTTTACCGAAgcg                                                                             >  1:1749889/1‑100 (MQ=255)
                          ccAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGcgcc                                                                           >  2:3590588/1‑100 (MQ=255)
                          ccAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCGTGGTGTCGGTTTTACCGAAGcgccg                                                                          >  2:3935357/1‑100 (MQ=255)
                          ccAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCGTGGTGTCGGTTTTACCGAAGcgccg                                                                          <  2:3251799/100‑1 (MQ=255)
                           cAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAATCTGGCCGACAGAGTGCGTGGTGTCGGTTTTACCGAAgcgctg                                                                          <  2:2887691/99‑3 (MQ=255)
                           cAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGcgccg                                                                          >  2:2279979/1‑100 (MQ=255)
                           cAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGcgccg                                                                          >  2:3580977/1‑100 (MQ=255)
                            aaCCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCGTGGTGTCGGTTTTACCGAAGCGCcgcg                                                                        >  1:1332903/1‑100 (MQ=255)
                              ccTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCGTGGTGTCGGTTTTACCGAAGCGCcgcgcg                                                                      >  2:1959909/1‑100 (MQ=255)
                               cTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGGCGACAGAGTGCGTGGTGTCGGTTTTACCGAAGCGCcgcgcg                                                                      >  1:3000584/1‑99 (MQ=255)
                               cTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCcgcgcg                                                                      >  1:2087325/1‑100 (MQ=255)
                                tGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCGTGGTGTCGGTTTTACCGAAGCGCCGCGCggg                                                                    <  2:3181977/100‑1 (MQ=255)
                                tGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGAAACCTGGCCGACAGAGTGCCGTGGCGGCGGTTGTACCGAAgcgctggacgg                                                                     >  1:2680345/1‑93 (MQ=255)
                                   aaaaCGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCg                                                                 >  2:287874/1‑100 (MQ=255)
                                    aaaCGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGt                                                                <  2:3455271/100‑1 (MQ=255)
                                       cGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTa                                                              <  2:437220/100‑1 (MQ=255)
                                        ggCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAg                                                             <  2:593572/100‑1 (MQ=255)
                                        ggCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAg                                                             <  1:3431874/100‑1 (MQ=255)
                                         gCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAgg                                                            <  2:91482/100‑1 (MQ=255)
                                         gCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAgg                                                            <  1:2795181/100‑1 (MQ=255)
                                         gCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAgg                                                            <  1:2682471/100‑1 (MQ=255)
                                         gCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGc                                                           <  1:1256203/100‑1 (MQ=255)
                                         gCAATCGCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGGGCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAggn                                                           <  1:2011778/100‑2 (MQ=255)
                                          cAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGc                                                           <  1:2053454/100‑1 (MQ=255)
                                           aaTCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCACCGCGCGGGGCGTTAGGcc                                                          <  2:1576235/100‑1 (MQ=255)
                                              ctggCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCAc                                                        <  1:3237244/96‑1 (MQ=255)
                                               tcGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTgg                                                     <  2:3684618/100‑1 (MQ=255)
                                                cGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTgg                                                     >  2:3889717/1‑100 (MQ=255)
                                                 gCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTggg                                                    >  1:3380684/1‑100 (MQ=255)
                                                   cACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGgcc                                                  <  1:1630719/100‑1 (MQ=255)
                                                     cTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGgccgcc                                               <  1:2853988/100‑1 (MQ=255)
                                                       gCTTCCACGGAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCAt                                             <  2:3281629/100‑1 (MQ=255)
                                                           ccACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTcg                                          <  2:2620499/100‑1 (MQ=255)
                                                           ccACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTcgc                                         <  2:1106187/100‑1 (MQ=255)
                                                           ccACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTcgc                                         >  1:1950505/1‑100 (MQ=255)
                                                            cACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTcgcg                                        >  1:918573/1‑100 (MQ=255)
                                                            cACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTcgcg                                        <  2:2321529/100‑1 (MQ=255)
                                                             aCGGAAAAATGGGAACCTGCAAGCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTcgcg                                        >  1:2473992/1‑100 (MQ=255)
                                                             aCGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTcgcg                                        <  1:393012/100‑1 (MQ=255)
                                                               ggAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATg                                     <  2:1158013/100‑1 (MQ=255)
                                                                gAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATg                                     <  2:4253472/100‑1 (MQ=255)
                                                                  aaaaTGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGCTGGc                                   >  1:3661207/1‑100 (MQ=255)
                                                                  aaaaTGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGc                                   >  1:706469/1‑100 (MQ=255)
                                                                  aaaaTGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGc                                   >  1:1790624/1‑100 (MQ=255)
                                                                  aaaaTGGGAACCTGCAACCTGGCCGACAGAGTGCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGGTAGGCCACTGGGCCGCCATTCGCGATGGCa                                  >  1:3848209/1‑100 (MQ=255)
                                                                     aTGGGAGCCTGCAACCTGGCCGACAGAGTGCGTGGTGGCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGTCCGCCATCGCGaatggcggc                               <  1:2614803/99‑9 (MQ=255)
                                                                     aTGGGAACCTGCAACCTGGCCGACAGAGTGCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGa                               >  2:1795367/1‑100 (MQ=255)
                                                                     aTGGGAACCTGCAACCTGGCCGACAGAGTGCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGa                               >  2:3101652/1‑100 (MQ=255)
                                                                      tGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTCGGCCACTGGGCCGCCATTCGCGCTGGCAAGa                               >  1:1004528/1‑100 (MQ=255)
                                                                        ccgatcTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGAtt                             <  1:3523709/95‑1 (MQ=255)
                                                                           aCCTGCAACCTGGCCGACAGAGTGCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGAtc                         <  1:2520741/100‑1 (MQ=255)
                                                                            ccTGCAACCTGGCCGACAGAGTGCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGAtct                        >  2:2909738/1‑100 (MQ=255)
                                                                            ccTGCAACCTGGCCGACAGAGTGCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGAtct                        <  2:1781501/100‑1 (MQ=255)
                                                                               gCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCTGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGATctct                      <  2:2052162/100‑1 (MQ=255)
                                                                               gCAACCTGGCCGACAGAGTGCCGTGGTGTCGCTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGATctct                      <  1:3182388/100‑1 (MQ=255)
                                                                                cAACCTGGCCGACAGAGTGCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGATCTCTAc                    <  1:1488331/100‑1 (MQ=255)
                                                                                 aaCCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGATCTCTa                     >  1:1048107/1‑99 (MQ=255)
                                                                                  aCCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGATCTCTAcc                   <  2:2437009/100‑1 (MQ=255)
                                                                                  aCCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGATCTCTAcc                   <  2:1734009/100‑1 (MQ=255)
                                                                                  aCCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGATCTCTAcc                   >  1:2429277/1‑100 (MQ=255)
                                                                                  aCCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGATCTCTAcc                   >  1:3900425/1‑100 (MQ=255)
                                                                                  aCCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGATCTCTAcc                   >  1:1431692/1‑100 (MQ=255)
                                                                                   ccTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGATCTCTACCa                  <  2:1445968/100‑1 (MQ=255)
                                                                                      ggCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGATCTCTACCAGTg               <  2:1522781/100‑1 (MQ=255)
                                                                                      ggCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGATCTCTACCAGTg               <  1:2898084/100‑1 (MQ=255)
                                                                                       gCCGACAGAGTGCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGATCTCTACCAGTGCg             <  2:522806/100‑1 (MQ=255)
                                                                                         cGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGATCTCTACCAGTGCgt            <  1:2991823/100‑1 (MQ=255)
                                                                                         cGACAGAGTGCCGTGGTGTCGCATTTACCGGAGCGCCACGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGATCTCTACCAGTg               <  1:1674111/97‑1 (MQ=255)
                                                                                         cGACAGAGTGCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGATCTCTACCAGTGCgtg           >  2:1679035/1‑100 (MQ=255)
                                                                                          gACAGAGTGCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGATCTCTACCAGTGCgtgg          <  1:223231/100‑1 (MQ=255)
                                                                                            cAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGATCTCTACCAGTGCgtggt         <  2:2473998/100‑1 (MQ=255)
                                                                                             agagTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGATCTCTACCAGTGCgtggtg        <  1:1015622/100‑1 (MQ=255)
                                                                                              gagTGCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGATCTCTACCAGTGCGTGGTGcc      <  1:99879/100‑1 (MQ=255)
                                                                                                gTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGATCTCTACCAGTGCGTGGTGCCg     >  1:1369101/1‑100 (MQ=255)
                                                                                                  gCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGATCTCTACCAGTGCGTGGTGCCGac   >  1:3385760/1‑100 (MQ=255)
                                                                                                  gCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGATCTCTACCAGTGCGTGGTGCCGac   >  1:4179075/1‑100 (MQ=255)
                                                                                                   ccGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGATCTCTACCAGTGCGTGGTGCCGacc  <  1:3769937/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCAGTATTTCTTCGACAAGCTGATGACCAACCTGAAAAACGGCAATCTCGCCACTGCTTCCACGGAAAAATGGGAACCTGCAACCTGGCCGACAGAGTGCCGTGGTGTCGGTTTTACCGAAGCGCCGCGCGGGGCGTTAGGCCACTGGGCCGCCATTCGCGATGGCAAGATTGATCTCTACCAGTGCGTGGTGCCGACC  >  NC_000913/1034651‑1034849

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