Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,007,187 A→T 32.5% 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,1870AT32.5% 212.5 / 149.9 166D309V (GAC→GTC) ygeWputative carbamoyltransferase
Reads supporting (aligned to +/- strand):  ref base A (59/53);  new base T (24/30);  total (83/83)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 167 total aligned reads displayed.

AAGTGATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAG  >  NC_000913/3007088‑3007285
                                                                                                    |                                                                                                  
   tGATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTg                                                                                                 <  1:4298268/100‑1 (MQ=255)
    gATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTgg                                                                                                <  1:1407825/100‑1 (MQ=255)
    gATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTgg                                                                                                <  2:238921/100‑1 (MQ=255)
         aaGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTg                                                                                           <  2:1587633/100‑1 (MQ=255)
         aaGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTg                                                                                           <  2:1910762/100‑1 (MQ=255)
         aaGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTg                                                                                           <  1:2610904/100‑1 (MQ=255)
         aaGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTg                                                                                           <  1:4011277/100‑1 (MQ=255)
          agagCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGt                                                                                          >  1:4096500/1‑100 (MQ=255)
          agagCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGt                                                                                          >  1:3005847/1‑100 (MQ=255)
            agCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTAc                                                                                        <  1:2149074/100‑1 (MQ=255)
            agCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTAc                                                                                        <  2:3635892/100‑1 (MQ=255)
            agCGTAATGAATTTCTGCGTCCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTAc                                                                                        <  1:926907/100‑1 (MQ=255)
               gTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGa                                                                                     <  1:3521086/100‑1 (MQ=255)
               gTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGa                                                                                     >  1:23848/1‑100 (MQ=255)
               gTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGa                                                                                     <  1:3221467/100‑1 (MQ=255)
                 aCTGAATTGCTGGGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAag                                                                                   <  1:991335/100‑1 (MQ=255)
                 aCTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGCACGCGCAACACAAAGTCTGGCATTGTACTGAag                                                                                   >  1:450317/1‑100 (MQ=255)
                 aCTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAag                                                                                   >  1:3088197/1‑100 (MQ=255)
                 aCTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAag                                                                                   >  2:608663/1‑100 (MQ=255)
                  cTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACCGAACGCGCAACACAAAGACTGACATTGTACTGAaga                                                                                  >  1:68488/1‑100 (MQ=255)
                     aaTTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAgat                                                                               <  2:2213487/100‑1 (MQ=255)
                         gctgcGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGatgatg                                                                           <  2:2140940/100‑1 (MQ=255)
                          ctgcGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATgg                                                                          <  1:1301003/100‑1 (MQ=255)
                           tgcgtgcgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGa                                                                         <  2:3685419/100‑1 (MQ=255)
                            gcgtgcgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGaa                                                                        <  2:4297704/100‑1 (MQ=255)
                             cgtgcgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAAc                                                                       <  2:1975063/100‑1 (MQ=255)
                              gtgcgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACt                                                                      <  1:1286441/100‑1 (MQ=255)
                              gtgcgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACt                                                                      >  2:1186450/1‑100 (MQ=255)
                               tgcgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTg                                                                     >  2:3370113/1‑100 (MQ=255)
                                gcgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGa                                                                    <  2:3973114/100‑1 (MQ=255)
                                gcgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGGACTGAAGAGATGATGGAACTGa                                                                    <  2:3563442/100‑1 (MQ=255)
                                 cgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGAc                                                                   >  1:2385399/1‑100 (MQ=255)
                                 cgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGa                                                                    >  1:753339/1‑99 (MQ=255)
                                  gAACGATCACGAAGGCTTAAAAGCACTGGAGAATCAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGAcc                                                                  >  1:2895105/1‑100 (MQ=255)
                                    aCGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCg                                                                <  2:735325/100‑1 (MQ=255)
                                    aCGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCg                                                                >  2:530161/1‑100 (MQ=255)
                                    aCGATCACGAAGGCTTAAAAGCACAGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCg                                                                <  1:2754367/100‑1 (MQ=255)
                                         cACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATg                                                           >  1:2290267/1‑100 (MQ=255)
                                           cGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGc                                                         <  2:1227470/100‑1 (MQ=255)
                                            gAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCg                                                        <  2:3721509/100‑1 (MQ=255)
                                             aaGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGa                                                       <  2:1797833/100‑1 (MQ=255)
                                             aaGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGa                                                       >  2:4283260/1‑100 (MQ=255)
                                               ggCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGaa                                                      >  1:4325643/1‑99 (MQ=255)
                                               ggCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAg                                                     <  1:388683/100‑1 (MQ=255)
                                                gCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGc                                                    >  2:890983/1‑100 (MQ=255)
                                                 cTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGcc                                                   >  1:2639677/1‑100 (MQ=255)
                                                 cTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGcc                                                   >  2:2675740/1‑100 (MQ=255)
                                                  ttAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGccc                                                  <  2:2834071/100‑1 (MQ=255)
                                                   tAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCt                                                 >  2:306386/1‑100 (MQ=255)
                                                   tAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCt                                                 >  1:2166844/1‑100 (MQ=255)
                                                     aaaGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGt                                               >  1:1691973/1‑100 (MQ=255)
                                                     aaaGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGt                                               <  1:1507628/100‑1 (MQ=255)
                                                     aaaGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGt                                               >  1:2837945/1‑100 (MQ=255)
                                                      aaGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTa                                              <  2:397961/100‑1 (MQ=255)
                                                      aaGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTa                                              >  2:1125169/1‑100 (MQ=255)
                                                       aGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTAc                                             <  2:697220/100‑1 (MQ=255)
                                                        gCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACa                                            <  2:2550011/100‑1 (MQ=255)
                                                        gCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACa                                            <  2:1025542/100‑1 (MQ=255)
                                                        gCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGAGAAGCCCTGTACa                                            >  2:3943277/1‑100 (MQ=255)
                                                         cACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACAt                                           >  1:2979076/1‑100 (MQ=255)
                                                         cACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACAt                                           >  2:2876984/1‑100 (MQ=255)
                                                          aCTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATg                                          >  1:929968/1‑100 (MQ=255)
                                                           cTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGc                                         >  1:3599355/1‑100 (MQ=255)
                                                           cTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGc                                         >  2:679372/1‑100 (MQ=255)
                                                           cTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGc                                         >  2:3226798/1‑100 (MQ=255)
                                                           cTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGc                                         <  2:3776765/100‑1 (MQ=255)
                                                            tGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCa                                        >  2:3410372/1‑100 (MQ=255)
                                                             ggAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCAc                                       >  2:2625291/1‑100 (MQ=255)
                                                             ggAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACCGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCAc                                       >  1:420044/1‑100 (MQ=255)
                                                              gAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACCGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCAct                                      >  1:1039222/1‑100 (MQ=255)
                                                              gAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCAct                                      <  2:3520036/100‑1 (MQ=255)
                                                                      gtgtCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCa                              >  2:1281836/1‑100 (MQ=255)
                                                                      gtgtCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGAGGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCa                              <  1:1959508/100‑1 (MQ=255)
                                                                        gtCTGGCACATAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGCTGCCGAAGCCCTGTACATGCACTGCCTGCCAGc                            <  1:2895486/100‑1 (MQ=255)
                                                                        gtCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGc                            <  2:3648366/100‑1 (MQ=255)
                                                                          cTGGCACAGAACGCGCAACACAAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCt                           >  2:3983561/1‑100 (MQ=255)
                                                                          cTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTg                          <  2:1324140/100‑1 (MQ=255)
                                                                          cTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTg                          <  1:1108070/100‑1 (MQ=255)
                                                                            ggCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGat                        <  2:4395955/100‑1 (MQ=255)
                                                                            ggCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGat                        >  2:1890518/1‑100 (MQ=255)
                                                                             gCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGata                       <  2:3122584/100‑1 (MQ=255)
                                                                               acaGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATc                     <  1:576610/100‑1 (MQ=255)
                                                                               acaGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATc                     <  2:2217650/100‑1 (MQ=255)
                                                                                caGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCa                    <  1:4107722/100‑1 (MQ=255)
                                                                                 aGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAg                   <  1:147682/100‑1 (MQ=255)
                                                                                  gAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAgc                  >  1:3129769/1‑100 (MQ=255)
                                                                                  gAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAgc                  <  2:4235895/100‑1 (MQ=255)
                                                                                   aaCGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAgcg                 <  1:3703197/100‑1 (MQ=255)
                                                                                   aaCGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAgc                  <  1:1605364/99‑1 (MQ=255)
                                                                                    atgcgcAACACAAAGTGTGGCATTGTACTGAAGAGATGATGAAACTTACCCATTCTGACGCATCCCGGTAGCTTCACTGCCTGCCAGCTGATATCAgcgg                >  1:1111993/3‑100 (MQ=255)
                                                                                     cgcgCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAgcggc               >  1:3860020/1‑100 (MQ=255)
                                                                                      gcgcAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAgcggcg              <  2:2907161/100‑1 (MQ=255)
                                                                                        gcAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTa            >  2:1499060/1‑100 (MQ=255)
                                                                                         cAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCGCTGATATCAgctggcaggc          <  2:54919/100‑9 (MQ=255)
                                                                                          aaCACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATc          <  2:3268291/100‑1 (MQ=255)
                                                                                           acacAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATcc         <  1:3272216/100‑1 (MQ=255)
                                                                                            cacaAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCt        >  1:969673/1‑100 (MQ=255)
                                                                                              caAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGt      >  1:2955822/1‑100 (MQ=255)
                                                                                                   gACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAaag  >  1:1849469/1‑100 (MQ=255)
                                                                                                   gACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAaag  <  2:1954073/100‑1 (MQ=255)
                                                                                                    |                                                                                                  
AAGTGATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAG  >  NC_000913/3007088‑3007285

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: