Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,007,187 A→T 84.7% D309V (GAC→GTC)  ygeW → putative carbamoyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,007,1870AT84.7% 287.6 / 41.1 118D309V (GAC→GTC) ygeWputative carbamoyltransferase
Reads supporting (aligned to +/- strand):  ref base A (6/12);  new base T (42/57);  total (48/70)
Fisher's exact test for biased strand distribution p-value = 6.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01

Only 100 of 119 total aligned reads displayed.

AAGTGATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGT  >  NC_000913/3007088‑3007281
                                                                                                   |                                                                                              
caGTGACGGAAGAGCGTCCTGAATTGCTGCGTGCGAACGATCCCGCAGGCTTAACAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGa                                                                                                <  2:1586655/99‑1 (MQ=255)
  gTGATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCt                                                                                              <  1:54841/100‑1 (MQ=255)
     aaggAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGc                                                                                           >  2:1981868/3‑100 (MQ=255)
     aTGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGc                                                                                           <  1:1351602/100‑1 (MQ=255)
         aaGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTg                                                                                       <  2:2896217/100‑1 (MQ=255)
             gCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACt                                                                                   <  1:59182/100‑1 (MQ=255)
              cGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTg                                                                                  <  1:412162/100‑1 (MQ=255)
              cGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTg                                                                                  >  2:8956/1‑100 (MQ=255)
               gTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGa                                                                                 >  1:1407958/1‑100 (MQ=255)
                tACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGaa                                                                                >  2:2338904/1‑100 (MQ=255)
                tACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGaa                                                                                >  2:1792824/1‑100 (MQ=255)
                 aCTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGCACGCGCAACACAAAGTCTGGCATTGTACTGAag                                                                               <  1:1927165/100‑1 (MQ=255)
                 aCTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAag                                                                               <  2:2260522/100‑1 (MQ=255)
                 aCTGAACTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAag                                                                               <  1:856785/100‑1 (MQ=255)
                  cTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAaga                                                                              >  1:2505645/1‑100 (MQ=255)
                  cTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAaga                                                                              >  2:1779618/1‑100 (MQ=255)
                  cTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAaga                                                                              >  1:2848731/1‑100 (MQ=255)
                     aaTTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAgat                                                                           <  2:2572821/100‑1 (MQ=255)
                     aaTTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAgat                                                                           <  1:975660/100‑1 (MQ=255)
                       ttGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAgatga                                                                         <  2:907184/100‑1 (MQ=255)
                       ttGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAgatga                                                                         >  1:26838/1‑100 (MQ=255)
                         gctgcGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGatgatg                                                                       >  2:1655833/1‑100 (MQ=255)
                         gctgcGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGatgatg                                                                       <  1:1872367/100‑1 (MQ=255)
                         gctgcGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCAGTGTACTGAAGAGatgatg                                                                       <  1:1894114/100‑1 (MQ=255)
                          ctgcGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATgg                                                                      >  2:2807079/1‑100 (MQ=255)
                            gcgtgcgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGaa                                                                    >  1:2674807/1‑100 (MQ=255)
                                gcgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAAc                                                                   >  1:2035353/1‑97 (MQ=255)
                                 cgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGAc                                                               >  1:554045/1‑100 (MQ=255)
                                 cgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGAc                                                               >  1:2187997/1‑100 (MQ=255)
                                    aCGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCg                                                            >  1:2340932/1‑100 (MQ=255)
                                    aCGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCg                                                            <  1:211151/100‑1 (MQ=255)
                                     cGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGt                                                           >  2:408183/1‑100 (MQ=255)
                                     cGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACACGt                                                           >  2:241387/1‑100 (MQ=255)
                                     cGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGt                                                           >  2:34846/1‑100 (MQ=255)
                                        tCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGAt                                                        <  1:1068567/100‑1 (MQ=255)
                                         cACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATg                                                       <  2:852913/100‑1 (MQ=255)
                                         cACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATGGTACGGAAGAGATGATGGAACTGACCCg                                                            >  1:2147572/1‑95 (MQ=255)
                                           cGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGc                                                     <  2:2065949/100‑1 (MQ=255)
                                               ggCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAg                                                 <  2:811019/100‑1 (MQ=255)
                                                gCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGc                                                <  2:1124973/100‑1 (MQ=255)
                                                 cTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGcc                                               >  1:1544931/1‑100 (MQ=255)
                                                 cTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGcc                                               <  2:1312305/100‑1 (MQ=255)
                                                 cTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGcc                                               >  1:2392594/1‑100 (MQ=255)
                                                  ttAAAAGCACTTGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATGGTACTGAAGAGAAGATGGAACTGAACCGTGATGGCGAAGccc                                              <  1:2255285/100‑1 (MQ=255)
                                                  ttAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGccc                                              <  1:1905615/100‑1 (MQ=255)
                                                   tAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCt                                             <  2:706774/100‑1 (MQ=255)
                                                      aaGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTa                                          <  2:1659766/100‑1 (MQ=255)
                                                      aaGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTa                                          <  2:1466314/100‑1 (MQ=255)
                                                      aaGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTa                                          <  2:554048/100‑1 (MQ=255)
                                                      aaGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTa                                          <  2:2745481/100‑1 (MQ=255)
                                                        gCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACa                                        >  2:2293705/1‑100 (MQ=255)
                                                         cACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACAt                                       >  2:2780413/1‑100 (MQ=255)
                                                           cTGGAACAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATCGTACTGAAGCGATGATGGAACTGACCCGTGATGGCGCAGCCCTGTACATGc                                     >  1:277103/1‑100 (MQ=255)
                                                           cTGGAAACACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTTGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATg                                      <  1:503335/99‑1 (MQ=255)
                                                            tGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCa                                    <  1:146482/100‑1 (MQ=255)
                                                                 aaaCAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCActgcc                               <  2:2293243/100‑1 (MQ=255)
                                                                   aCAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCActgcctg                             <  2:307324/100‑1 (MQ=255)
                                                                    cAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCActgcctgc                            <  2:2909875/100‑1 (MQ=255)
                                                                    cAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCActgcctg                             <  2:194739/99‑1 (MQ=255)
                                                                    cAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCActgcctgc                            <  2:2238994/100‑1 (MQ=255)
                                                                    cAGGGTCTGGCACAGAACGCTCAACACAAAGGCTGGCATGGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTTTACATGCActgcctgc                            <  1:1837628/100‑1 (MQ=255)
                                                                     aGTGTCTGGCACAGAACGCGCAACACCAAGTCTGGCATCGTACTGAAGAGATGATAGAACTGACCCGTGATGGCGAAGCCCTGTACATGCAAtgcctgcc                           >  1:874111/1‑100 (MQ=255)
                                                                     aGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCActgcctgcn                           <  1:292454/100‑2 (MQ=255)
                                                                     aGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACtgcctgcc                           <  1:63730/100‑1 (MQ=255)
                                                                        gtCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGc                        <  1:975049/100‑1 (MQ=255)
                                                                        gtCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGc                        <  2:2801636/100‑1 (MQ=255)
                                                                          cTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGCGATGATGGAACCGACCCGTGATGGCGAAGCGGTGTACATGCACTGCCTGCCAGCTg                      >  1:2283008/1‑100 (MQ=255)
                                                                          cTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTg                      >  2:1484246/1‑100 (MQ=255)
                                                                          cTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTg                      >  1:1217684/1‑100 (MQ=255)
                                                                          cTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTg                      >  2:712080/1‑100 (MQ=255)
                                                                           tGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGa                     <  2:2812021/100‑1 (MQ=255)
                                                                           tGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGa                     <  2:2114328/100‑1 (MQ=255)
                                                                            ggCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGat                    >  2:2479283/1‑100 (MQ=255)
                                                                            ggCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGat                    >  1:934083/1‑100 (MQ=255)
                                                                              cacaGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGatat                  >  2:2488029/1‑100 (MQ=255)
                                                                              cacaGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATACACTGCCTGCCAGCTGatat                  <  2:614967/100‑1 (MQ=255)
                                                                               acaGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATc                 <  1:2752657/100‑1 (MQ=255)
                                                                                 aGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAg               <  1:1453189/100‑1 (MQ=255)
                                                                                 aGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAg               <  1:3336/100‑1 (MQ=255)
                                                                                  gAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAgc              >  2:1502498/1‑100 (MQ=255)
                                                                                   aaCGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAgcg             <  2:122157/100‑1 (MQ=255)
                                                                                    aCGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAgcgg            >  1:1569965/1‑100 (MQ=255)
                                                                                     cgcgCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAgcggc           <  1:2584916/100‑1 (MQ=255)
                                                                                      gcgcAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAgcggcg          <  2:27310/100‑1 (MQ=255)
                                                                                       cgcAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGt         >  2:932587/1‑100 (MQ=255)
                                                                                       cgcAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGt         >  2:1823962/1‑100 (MQ=255)
                                                                                        gcAACACAAAGTCTGGCATTGTACTGGAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTa        <  1:2230232/100‑1 (MQ=255)
                                                                                        gcAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCGCTGATATCAgctggcagg       <  2:1039265/100‑8 (MQ=255)
                                                                                         cAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTAt       >  1:1015066/1‑100 (MQ=255)
                                                                                         cAAAACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCGCTGATATCAgctggcaggc      <  1:2718257/100‑9 (MQ=255)
                                                                                           acacAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATcc     <  2:1544937/100‑1 (MQ=255)
                                                                                           acacAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATcc     <  1:1189825/100‑1 (MQ=255)
                                                                                           acacAAAGACTGGCATTGTACAGAATAGAGGATGGACCTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAgcggcg          <  1:2624212/95‑1 (MQ=255)
                                                                                             acaAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTg   <  1:932584/100‑1 (MQ=255)
                                                                                             acaAAGTCNGGCATTGTACTGAAGAGNTGATNGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGNATCCTg   >  2:2555550/1‑100 (MQ=255)
                                                                                             acaAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTg   >  1:416372/1‑100 (MQ=255)
                                                                                              caAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGt  >  1:2105551/1‑100 (MQ=255)
                                                                                              caAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCCTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGt  >  1:1589458/1‑100 (MQ=255)
                                                                                              caAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGt  <  2:460645/100‑1 (MQ=255)
                                                                                                aaGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGt  >  2:19352/1‑98 (MQ=255)
                                                                                                   |                                                                                              
AAGTGATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGT  >  NC_000913/3007088‑3007281

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: