Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,007,187 A→T 82.4% 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,1870AT82.4% 255.0 / 41.0 106D309V (GAC→GTC) ygeWputative carbamoyltransferase
Reads supporting (aligned to +/- strand):  ref base A (9/10);  new base T (38/49);  total (47/59)
Fisher's exact test for biased strand distribution p-value = 8.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

GTGATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTA  >  NC_000913/3007090‑3007282
                                                                                                 |                                                                                               
gTGATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACt                                                                                               <  2:1570361/100‑1 (MQ=255)
 tGATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTg                                                                                              <  1:1247803/100‑1 (MQ=255)
   aTGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGc                                                                                            >  2:2094700/1‑100 (MQ=255)
   aTGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGc                                                                                            <  1:316544/100‑1 (MQ=255)
     ggAAGAGCTTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCAt                                                                                          <  1:2169806/100‑1 (MQ=255)
      gAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCAtt                                                                                         >  1:2065570/1‑100 (MQ=255)
      gAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCAtt                                                                                         >  2:2000303/1‑100 (MQ=255)
       aaGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTg                                                                                        <  2:935113/100‑1 (MQ=255)
       aaGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTg                                                                                        <  2:1658745/100‑1 (MQ=255)
          agCGTACTGAATTTCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTAc                                                                                     <  2:948711/100‑1 (MQ=255)
           gCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACt                                                                                    >  2:440164/1‑100 (MQ=255)
             gTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGa                                                                                  <  2:1767966/100‑1 (MQ=255)
             gTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGa                                                                                  <  1:2311453/100‑1 (MQ=255)
             gTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGa                                                                                  >  1:850104/1‑100 (MQ=255)
              tACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGaa                                                                                 >  2:61049/1‑100 (MQ=255)
                cTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAaga                                                                               >  2:1223238/1‑100 (MQ=255)
                cTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAaga                                                                               <  2:1032382/100‑1 (MQ=255)
                cTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACGGGCAAAACAGTGTCTGGCACAGCACGCGCAACAGAAAGTCTGGCCTTGTACTGacg                                                                                >  2:226981/1‑97 (MQ=255)
                 tGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAagag                                                                              <  1:1963123/100‑1 (MQ=255)
                  gAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAgaga                                                                             >  2:848762/1‑100 (MQ=255)
                  gAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAgaga                                                                             <  1:2269397/100‑1 (MQ=255)
                  gAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAgaga                                                                             >  1:819929/1‑100 (MQ=255)
                   aaTTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAgat                                                                            >  2:192155/1‑100 (MQ=255)
                    aTTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAgatg                                                                           >  1:67934/1‑100 (MQ=255)
                       gctgcGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGatgatg                                                                        <  2:850099/100‑1 (MQ=255)
                        ctgcGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATgg                                                                       <  1:846998/100‑1 (MQ=255)
                         tgcgtgcgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGatgatg                                                                        >  2:1123299/1‑98 (MQ=255)
                           cgtgcgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAAc                                                                    <  2:1965356/100‑1 (MQ=255)
                           cgtgcgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCGGGCATTGTACTGAAGAGATGATGGAAc                                                                    <  2:85134/100‑1 (MQ=255)
                            gtgcgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACt                                                                   >  1:1706946/1‑100 (MQ=255)
                               cgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGAc                                                                >  1:1629383/1‑100 (MQ=255)
                               cgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGAc                                                                >  2:1210539/1‑100 (MQ=255)
                               cgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGAc                                                                >  1:1349874/1‑100 (MQ=255)
                                gAACGATCACGAAGGCTTAAAAGCACTGGAAAAACCGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGAcc                                                               >  2:1156606/1‑100 (MQ=255)
                                gAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGAcc                                                               >  2:89557/1‑100 (MQ=255)
                                  aCGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCg                                                             >  2:1680358/1‑100 (MQ=255)
                                   cGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGt                                                            >  2:2180317/1‑100 (MQ=255)
                                   cGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGt                                                            >  1:1012951/1‑100 (MQ=255)
                                   cGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGt                                                            >  1:1307618/1‑100 (MQ=255)
                                   cGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGt                                                            >  1:1609500/1‑100 (MQ=255)
                                   cGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGt                                                            >  2:2542336/1‑100 (MQ=255)
                                    gATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTg                                                           <  1:1129170/100‑1 (MQ=255)
                                     aTCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGa                                                          >  1:1886809/1‑100 (MQ=255)
                                     aTCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGa                                                          >  2:194043/1‑100 (MQ=255)
                                      tCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGAt                                                         <  1:1998709/100‑1 (MQ=255)
                                      tCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGAt                                                         <  2:1007159/100‑1 (MQ=255)
                                        aCGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATgg                                                       >  1:649106/1‑100 (MQ=255)
                                         cGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGc                                                      >  1:2443424/1‑100 (MQ=255)
                                          gAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCg                                                     <  1:1398856/100‑1 (MQ=255)
                                           aaGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGa                                                    <  2:844221/100‑1 (MQ=255)
                                           aaGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGa                                                    <  2:844084/100‑1 (MQ=255)
                                           aaGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGa                                                    >  1:2201559/1‑100 (MQ=255)
                                             ggCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAg                                                  >  2:2557795/1‑100 (MQ=255)
                                              gCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGc                                                 <  2:888935/100‑1 (MQ=255)
                                              gCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGc                                                 >  2:2373730/1‑100 (MQ=255)
                                              gCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATCGCGAAGc                                                 <  1:2316223/100‑1 (MQ=255)
                                               cTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGcc                                                >  2:1543157/1‑100 (MQ=255)
                                                ttAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGccc                                               <  1:300207/100‑1 (MQ=255)
                                                 tAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGCTGGCGAAGccat                                              >  2:1130250/1‑98 (MQ=255)
                                                  aaaaGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTg                                             <  2:559284/100‑1 (MQ=255)
                                                    aaGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTa                                           <  1:2094711/100‑1 (MQ=255)
                                                    aaGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTa                                           <  1:760586/100‑1 (MQ=255)
                                                     aGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTAc                                          <  2:740644/100‑1 (MQ=255)
                                                      gCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACa                                         <  1:1225148/100‑1 (MQ=255)
                                                        aCTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATg                                       <  1:2257046/100‑1 (MQ=255)
                                                         cTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGc                                      <  1:1123304/100‑1 (MQ=255)
                                                          tGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCa                                     <  1:2504717/100‑1 (MQ=255)
                                                           ggAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCAc                                    >  2:1186064/1‑100 (MQ=255)
                                                           ggAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCAc                                    >  1:310326/1‑100 (MQ=255)
                                                            gAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCAct                                   <  2:1227202/100‑1 (MQ=255)
                                                               aaaCAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCActgcc                                <  2:67934/100‑1 (MQ=255)
                                                               aaaCAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCActgcc                                <  1:1676484/100‑1 (MQ=255)
                                                                 aCAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCActgcctg                              <  2:529060/100‑1 (MQ=255)
                                                                   aGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACtgcctgcc                            <  1:64605/100‑1 (MQ=255)
                                                                    gtgtCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCa                           <  2:1834953/100‑1 (MQ=255)
                                                                    gtgtCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCa                           <  1:2326933/100‑1 (MQ=255)
                                                                     tgtCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAg                          <  1:119300/100‑1 (MQ=255)
                                                                     tgtCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAg                          <  2:2490493/100‑1 (MQ=255)
                                                                      gtCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGc                         >  1:823262/1‑100 (MQ=255)
                                                                      gtCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGc                         >  1:68050/1‑100 (MQ=255)
                                                                      gtCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGc                         >  2:107609/1‑100 (MQ=255)
                                                                        cTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTg                       <  2:570758/100‑1 (MQ=255)
                                                                        cTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTg                       <  2:1453903/100‑1 (MQ=255)
                                                                         tGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGa                      >  2:2149885/1‑100 (MQ=255)
                                                                         tGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGa                      <  2:1041718/100‑1 (MQ=255)
                                                                         tGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGa                      >  2:2083914/1‑100 (MQ=255)
                                                                           gCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGata                    >  1:2544400/1‑100 (MQ=255)
                                                                           gCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGata                    <  1:1223243/100‑1 (MQ=255)
                                                                                   cgcgCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAgcggc            <  1:349660/100‑1 (MQ=255)
                                                                                   cgcgCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATACCAgcggc            >  1:2198836/1‑100 (MQ=255)
                                                                                    gcgcAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAgcggcg           <  2:1911869/100‑1 (MQ=255)
                                                                                    gcgcAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAgcggcg           <  1:649427/100‑1 (MQ=255)
                                                                                       cAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTAt        <  1:1907393/100‑1 (MQ=255)
                                                                                       cAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTAt        >  1:394704/1‑100 (MQ=255)
                                                                                         acacAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATcc      <  2:206254/100‑1 (MQ=255)
                                                                                         acacAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATcc      <  2:1620084/100‑1 (MQ=255)
                                                                                         acacAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATcc      <  1:2558153/100‑1 (MQ=255)
                                                                                         acacAAAGTCTGGCAGTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATcc      <  1:2557808/100‑1 (MQ=255)
                                                                                         acacAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATcc      <  1:547521/100‑1 (MQ=255)
                                                                                             aaaGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTa  >  1:1219028/1‑100 (MQ=255)
                                                                                                 |                                                                                               
GTGATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTA  >  NC_000913/3007090‑3007282

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