Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 908,807 A→C L163R (CTG→CGG)  ltaE ← L‑allo‑threonine aldolase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913908,8070AC96.5% 479.4 / 2.4 141L163R (CTG→CGG) ltaEL‑allo‑threonine aldolase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base C (62/74);  total (65/76)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 143 total aligned reads displayed.

CAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATT  >  NC_000913/908710‑908903
                                                                                                 |                                                                                                
cAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACTGCATTAAAGATGCGCGCGCCGTCAACATGCCgc                                                                                                >  2:3586584/1‑100 (MQ=255)
 aGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCcgcg                                                                                               >  2:379557/1‑100 (MQ=255)
  ggCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCgcgc                                                                                              <  2:2174842/100‑1 (MQ=255)
  ggCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCgcgc                                                                                              >  1:2447783/1‑100 (MQ=255)
  ggCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCgcgc                                                                                              >  1:4014858/1‑100 (MQ=255)
   gCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGcc                                                                                             <  1:3691087/100‑1 (MQ=255)
        tGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGAtt                                                                                        >  2:1985504/1‑100 (MQ=255)
        tGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACAGTCAACATGCAGCGCCAGAtt                                                                                        >  2:1794463/1‑100 (MQ=255)
         ggTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTg                                                                                       <  1:175313/100‑1 (MQ=255)
         ggTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTg                                                                                       <  1:90772/100‑1 (MQ=255)
         ggTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTg                                                                                       <  2:2106017/100‑1 (MQ=255)
          gTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTgc                                                                                      >  2:3250732/1‑100 (MQ=255)
           tGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTgcg                                                                                     <  2:3579771/100‑1 (MQ=255)
           tGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTgcg                                                                                     <  1:2311431/100‑1 (MQ=255)
           tGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTgcg                                                                                     <  2:3097233/100‑1 (MQ=255)
             aaCGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGTCGCGCCAGATTGCGCt                                                                                   <  1:2669658/100‑1 (MQ=255)
                 aaTCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTcgcg                                                                               <  1:1659569/100‑1 (MQ=255)
                 aaTCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTcgcg                                                                               <  2:1908731/100‑1 (MQ=255)
                 aaTCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTcgcg                                                                               <  2:1852767/100‑1 (MQ=255)
                  aTCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCgg                                                                              <  1:3083520/100‑1 (MQ=255)
                   tCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCggg                                                                             <  2:2116146/100‑1 (MQ=255)
                   tCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCggg                                                                             <  1:2883947/100‑1 (MQ=255)
                   tCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCggg                                                                             >  1:2284323/1‑100 (MQ=255)
                   tCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCggg                                                                             >  2:1470406/1‑100 (MQ=255)
                   tCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCggg                                                                             >  2:1289515/1‑100 (MQ=255)
                    cacaATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGt                                                                            >  1:1647690/1‑100 (MQ=255)
                     acaATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTa                                                                           >  1:181562/1‑100 (MQ=255)
                     acaATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTa                                                                           <  2:736681/100‑1 (MQ=255)
                       aaTATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTaaa                                                                         >  1:57435/1‑100 (MQ=255)
                            tGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTccc                                                                    <  2:2971351/100‑1 (MQ=255)
                            tGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTccc                                                                    >  2:256036/1‑100 (MQ=255)
                             tcGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCa                                                                   <  1:1883604/99‑1 (MQ=255)
                             gCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCa                                                                   <  1:3191879/100‑1 (MQ=255)
                             gCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCa                                                                   >  2:2679075/1‑100 (MQ=255)
                             gCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCCCGCACCGTCAACATGCCGCGGCAGATTGCGCTCGCGGGTAAATTCCCa                                                                   <  2:3980635/100‑1 (MQ=255)
                             atctGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCa                                                                   <  2:2780979/97‑1 (MQ=255)
                              cGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCaa                                                                  >  2:1853403/1‑99 (MQ=255)
                              cGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCAt                                                                  <  2:3265180/100‑1 (MQ=255)
                              cGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCAt                                                                  <  2:3467082/100‑1 (MQ=255)
                              cGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCAt                                                                  <  1:1102148/100‑1 (MQ=255)
                              cGTGATCTCTTTCACTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCCCTCGCGGGTAAATTCCCAt                                                                  <  2:3662176/100‑1 (MQ=255)
                               gTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATg                                                                 <  2:954662/100‑1 (MQ=255)
                               gTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATg                                                                 >  1:3503167/1‑100 (MQ=255)
                                tGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGc                                                                <  1:4060549/100‑1 (MQ=255)
                                 gATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCCCCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATCct                                                               >  2:2805750/1‑100 (MQ=255)
                                 gATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGct                                                               <  1:1789003/100‑1 (MQ=255)
                                 gATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGct                                                               <  1:1133306/100‑1 (MQ=255)
                                 gATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCaagat                                                               >  1:2823609/1‑96 (MQ=255)
                                  aTCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGctt                                                              <  2:319274/100‑1 (MQ=255)
                                       tttCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTc                                                         <  2:3148866/100‑1 (MQ=255)
                                       tttCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTc                                                         <  1:2464974/100‑1 (MQ=255)
                                        ntCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCa                                                        >  1:3235055/2‑100 (MQ=255)
                                         tCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAg                                                       <  1:821541/100‑1 (MQ=255)
                                          cAGTTCGCAGCCGTACGCCACCACGGCATGAAAGATGCGCGCACCGTCAACATGCCGCGCCAGAGTGCGCTCGCGGGTACATTCCCATGCTTCTTTCAgg                                                      <  2:3777213/100‑1 (MQ=255)
                                          cAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAgg                                                      <  2:3862516/100‑1 (MQ=255)
                                          cAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAgg                                                      <  1:2995476/100‑1 (MQ=255)
                                            gTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTa                                                    >  2:1241850/1‑100 (MQ=255)
                                            gTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTa                                                    <  1:2389771/100‑1 (MQ=255)
                                                 cAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGGGCTCGCGGGTAACTTCCCGTGCTTCTTTCAGGTATTccc                                               >  1:3820871/1‑100 (MQ=255)
                                                    ccGTAAGCCACCACGGCATTAAATATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCcgcg                                            <  2:1559410/100‑1 (MQ=255)
                                                     cGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCgg                                           <  2:852992/100‑1 (MQ=255)
                                                      gTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGc                                          >  1:1093286/1‑100 (MQ=255)
                                                       tAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCa                                         <  1:366148/100‑1 (MQ=255)
                                                       tAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCa                                         >  2:2882785/1‑100 (MQ=255)
                                                           ccaccaCGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAacac                                     >  2:2025468/1‑100 (MQ=255)
                                                           ccaccaCGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAacac                                     <  1:1123876/100‑1 (MQ=255)
                                                            caccacGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACt                                    <  2:777937/100‑1 (MQ=255)
                                                            caccacGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACt                                    >  2:2574816/1‑100 (MQ=255)
                                                              ccacGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACttt                                  >  1:732964/1‑100 (MQ=255)
                                                               cacGGCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTg                                 <  2:2152462/100‑1 (MQ=255)
                                                                  ggCATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCg                              >  1:739375/1‑100 (MQ=255)
                                                                    cATTAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGtt                            <  2:1808591/100‑1 (MQ=255)
                                                                      ttAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTtgt                          >  2:318693/1‑100 (MQ=255)
                                                                       tAAAGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTtgtg                         <  2:259083/100‑1 (MQ=255)
                                                                         aaGATGCGCGCACTGTCAACATGCCGCACCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTAGCCGTTGTggg                       >  1:3598587/1‑100 (MQ=255)
                                                                         aaGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTggg                       <  1:1106714/100‑1 (MQ=255)
                                                                         aaGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTggg                       <  2:479208/100‑1 (MQ=255)
                                                                         aaGATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTggg                       <  2:3374572/100‑1 (MQ=255)
                                                                           gATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTCTTCCCGCGGCAACACTTTGCCGTTGTGGgtg                     >  2:3862184/1‑100 (MQ=255)
                                                                           gATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGgtg                     >  2:708329/1‑100 (MQ=255)
                                                                           gATGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGgtg                     >  2:3875042/1‑100 (MQ=255)
                                                                           caTGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGgtg                     <  1:4080858/99‑1 (MQ=255)
                                                                             tGCGCGCACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGtt                   <  2:1024539/100‑1 (MQ=255)
                                                                               cgcgcACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGtttt                 >  2:1417334/1‑100 (MQ=255)
                                                                                gcgcACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTc                <  2:786478/100‑1 (MQ=255)
                                                                                   cACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAg             <  2:133727/100‑1 (MQ=255)
                                                                                   cACCGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAg             <  2:1345746/100‑1 (MQ=255)
                                                                                     ccGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGAc           <  1:2449331/100‑1 (MQ=255)
                                                                                      cGTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACt          >  2:2342002/1‑100 (MQ=255)
                                                                                       gTCGACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATTCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTg         <  1:2886954/100‑1 (MQ=255)
                                                                                       gTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTg         <  1:1847871/100‑1 (MQ=255)
                                                                                       gTCAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTg         >  1:3612606/1‑100 (MQ=255)
                                                                                         cAACATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAg       <  1:3381991/100‑1 (MQ=255)
                                                                                          aaCATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAGt      >  1:2664197/1‑100 (MQ=255)
                                                                                           aCATGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACCTTGCCGTTGTGGGTGTTTTCCAGACTGAGTa     <  2:2751207/100‑1 (MQ=255)
                                                                                            cATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTaa    >  1:529818/1‑100 (MQ=255)
                                                                                             aTGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAAt   <  1:2272597/100‑1 (MQ=255)
                                                                                              tGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAAtt  >  1:3811824/1‑100 (MQ=255)
                                                                                              tGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAAtt  >  1:3452955/1‑100 (MQ=255)
                                                                                              tGCCGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAAtt  >  2:27352/1‑100 (MQ=255)
                                                                                                 |                                                                                                
CAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATT  >  NC_000913/908710‑908903

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: ‑