Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,007,187 A→T 65.1% 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,1870AT65.1% 172.7 / 165.7 156D309V (GAC→GTC) ygeWputative carbamoyltransferase
Reads supporting (aligned to +/- strand):  ref base A (24/30);  new base T (42/59);  total (67/89)
Fisher's exact test for biased strand distribution p-value = 7.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

Only 100 of 156 total aligned reads displayed.

TGATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTA  >  NC_000913/3007091‑3007282
                                                                                                |                                                                                               
tGATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTg                                                                                              >  1:3280477/1‑100 (MQ=255)
tGATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTg                                                                                              >  1:3415661/1‑100 (MQ=255)
 gATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTgg                                                                                             <  1:2254243/100‑1 (MQ=255)
  aTGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGc                                                                                            <  1:1052947/100‑1 (MQ=255)
   tGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCa                                                                                           >  2:839558/1‑100 (MQ=255)
   tGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCa                                                                                           <  1:1694686/100‑1 (MQ=255)
   tGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCa                                                                                           <  2:2262985/100‑1 (MQ=255)
     gAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCAtt                                                                                         >  1:2728747/1‑100 (MQ=255)
      aaGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTg                                                                                        <  2:2496630/100‑1 (MQ=255)
       agagCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGt                                                                                       <  2:3058542/100‑1 (MQ=255)
          gCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACt                                                                                    >  1:1035238/1‑100 (MQ=255)
          gCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACt                                                                                    >  2:2577851/1‑100 (MQ=255)
          gCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTAc                                                                                     <  1:3271695/99‑1 (MQ=255)
              aCTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAag                                                                                <  2:3146587/100‑1 (MQ=255)
               cTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGTATTGTACTGAaga                                                                               <  1:1294905/100‑1 (MQ=255)
               cTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAaga                                                                               >  2:286541/1‑100 (MQ=255)
                  aaTTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAgat                                                                            <  2:1148486/100‑1 (MQ=255)
                   aTTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAgatg                                                                           >  2:1139488/1‑100 (MQ=255)
                      gctgcGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGatgatg                                                                        <  2:310109/100‑1 (MQ=255)
                      gctgcGTGCGAACGATCACGAAGGCTTAAAAAAAAAAAAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGatgatg                                                                        <  1:497204/100‑1 (MQ=255)
                        tgcgtgcgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGa                                                                      <  2:210780/100‑1 (MQ=255)
                          cgtgcgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAAc                                                                    <  1:1610821/100‑1 (MQ=255)
                          cgtgcgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAAc                                                                    <  1:1166466/100‑1 (MQ=255)
                          cgtgcgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAAc                                                                    <  2:1589350/100‑1 (MQ=255)
                           gtgcgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACt                                                                   >  2:2398734/1‑100 (MQ=255)
                           gtgcgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACt                                                                   <  1:3090905/100‑1 (MQ=255)
                            tgcgAACGATCACGAAGGCTTAAAGGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTg                                                                  <  1:1631526/100‑1 (MQ=255)
                            tgcgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACAAAAAGACTGGCATTGTACTGAAGAGATGATGGAACTg                                                                  >  1:1422987/1‑100 (MQ=255)
                             gcgAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGa                                                                 >  2:584776/1‑100 (MQ=255)
                                 aCGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCg                                                             <  2:3415687/100‑1 (MQ=255)
                                 aCGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCg                                                             <  1:1036541/100‑1 (MQ=255)
                                  cGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACAGACCCGt                                                            >  1:3159204/1‑100 (MQ=255)
                                  cGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGt                                                            >  1:491106/1‑100 (MQ=255)
                                   gATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTg                                                           >  2:1765079/1‑100 (MQ=255)
                                   gATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTg                                                           <  1:584772/100‑1 (MQ=255)
                                   gATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTg                                                           >  1:600141/1‑100 (MQ=255)
                                    aTCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGa                                                          >  2:1992602/1‑100 (MQ=255)
                                      cACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATg                                                        <  2:75405/100‑1 (MQ=255)
                                       aCGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATgg                                                       >  1:2856260/1‑100 (MQ=255)
                                        cGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGc                                                      >  1:197706/1‑100 (MQ=255)
                                        cGAAGGCTTAAAACCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGCCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGc                                                      >  1:1456092/1‑100 (MQ=255)
                                         gAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCg                                                     <  1:812135/100‑1 (MQ=255)
                                          aaGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGgaga                                                    >  1:2958569/1‑100 (MQ=255)
                                            ggCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAg                                                  <  1:99251/100‑1 (MQ=255)
                                            ggCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAg                                                  >  2:1315863/1‑100 (MQ=255)
                                              cTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGcc                                                >  2:69992/1‑100 (MQ=255)
                                              cTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGcc                                                >  2:1659233/1‑100 (MQ=255)
                                              cTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGcc                                                >  1:2161358/1‑100 (MQ=255)
                                               ttAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGccc                                               <  2:3449572/100‑1 (MQ=255)
                                                tAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCt                                              >  1:895196/1‑100 (MQ=255)
                                                tAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCt                                              >  2:1941206/1‑100 (MQ=255)
                                                 aaaaGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTg                                             <  2:661573/100‑1 (MQ=255)
                                                 aaaaGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTg                                             <  2:3280501/100‑1 (MQ=255)
                                                  aaaGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGt                                            >  2:1279857/1‑100 (MQ=255)
                                                   aaGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTa                                           <  2:627833/100‑1 (MQ=255)
                                                   aaGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTa                                           <  1:67286/100‑1 (MQ=255)
                                                    aGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTAc                                          >  2:2147027/1‑100 (MQ=255)
                                                     gCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATCGTACTGAAGAGGTGATGGAACTGACCCGTGATGGCGAAGCCCTGTACa                                         >  2:27473/1‑100 (MQ=255)
                                                       aCTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATg                                       <  2:2625484/100‑1 (MQ=255)
                                                       aCTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATg                                       <  2:2093404/100‑1 (MQ=255)
                                                        cTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGc                                      <  2:486614/100‑1 (MQ=255)
                                                         tGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCa                                     >  2:2269331/1‑100 (MQ=255)
                                                         tGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCa                                     >  1:2444698/1‑100 (MQ=255)
                                                          ggAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCAc                                    >  2:3163393/1‑100 (MQ=255)
                                                              aaaCAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCActgcc                                <  1:469327/100‑1 (MQ=255)
                                                                aCAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCActgcctg                              <  1:1523511/100‑1 (MQ=255)
                                                                aCAGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCActgcctg                              <  2:1037853/100‑1 (MQ=255)
                                                                aCAGGGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCActgcctg                              <  1:760604/100‑1 (MQ=255)
                                                                  aGTGTCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACtgcctgcc                            <  2:799800/100‑1 (MQ=255)
                                                                     gtCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGc                         <  1:3434693/100‑1 (MQ=255)
                                                                     gtCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGc                         >  2:1219423/1‑100 (MQ=255)
                                                                     gtCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACACGCACTGCCTGCCAGc                         >  2:893574/1‑100 (MQ=255)
                                                                     gtCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGc                         <  2:2529137/100‑1 (MQ=255)
                                                                      tCTGGCACAGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCt                        <  2:3017341/100‑1 (MQ=255)
                                                                       cTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTg                       >  2:1865449/1‑100 (MQ=255)
                                                                        tGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGa                      >  2:2427444/1‑100 (MQ=255)
                                                                          gaacaGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTg                       <  1:1872216/95‑1 (MQ=255)
                                                                           cacaGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGatat                   >  1:448063/1‑100 (MQ=255)
                                                                             caGAACGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCa                 >  2:1846011/1‑100 (MQ=255)
                                                                                 aCGCGCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAgcgg             <  2:2856283/100‑1 (MQ=255)
                                                                                  cgctCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACAGGCACTGCCTGCCAGCTGATATCAgcggc            <  1:1554747/100‑1 (MQ=255)
                                                                                  cgcgCAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAgcggc            >  2:3173349/1‑100 (MQ=255)
                                                                                   gcgcAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAgcggcg           <  2:545253/100‑1 (MQ=255)
                                                                                   gcgcAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAgcggcg           <  2:1792767/100‑1 (MQ=255)
                                                                                   gcgcAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAgcggcg           <  2:3099146/100‑1 (MQ=255)
                                                                                   gcgcAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAgcggcg           <  2:2743867/100‑1 (MQ=255)
                                                                                    cgcAACACAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGt          >  1:261622/1‑100 (MQ=255)
                                                                                     gcAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTa         <  1:2491317/100‑1 (MQ=255)
                                                                                     gcAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTa         >  2:2211684/1‑100 (MQ=255)
                                                                                        acacAAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATcc      >  1:654072/1‑100 (MQ=255)
                                                                                        acacAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATcc      <  1:1069509/100‑1 (MQ=255)
                                                                                         cacaAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCt     >  1:1708590/1‑100 (MQ=255)
                                                                                         cacaAAGACTGGCATTGTACTGAAGAGCTGATGGAACTGACCCGTGATGCCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCt     >  1:1570055/1‑100 (MQ=255)
                                                                                          acaAAGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTg    <  2:189796/100‑1 (MQ=255)
                                                                                          acaAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTg    <  2:2566225/100‑1 (MQ=255)
                                                                                           caAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACTCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGt   >  1:471736/1‑100 (MQ=255)
                                                                                           caAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGt   >  1:2642284/1‑100 (MQ=255)
                                                                                           caAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGt   >  2:964335/1‑100 (MQ=255)
                                                                                            aaaGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTa  <  1:1992585/100‑1 (MQ=255)
                                                                                            aaaGTCTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACCGCCTGCCAGCTGATATCAGCGGCGTATCCTGTa  <  2:233136/100‑1 (MQ=255)
                                                                                                |                                                                                               
TGATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTA  >  NC_000913/3007091‑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: