Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 125,226 G→T 6.5% R737L (CGT→CTT)  aceE → pyruvate dehydrogenase, decarboxylase component E1, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913125,2260GT 6.5% 986.5 / 41.6 275R737L (CGT→CTT) aceEpyruvate dehydrogenase, decarboxylase component E1, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base G (131/125);  new base T (8/10);  total (139/136)
Fisher's exact test for biased strand distribution p-value = 6.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01

Only 100 of 276 total aligned reads displayed.

TACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGC  >  NC_000913/125078‑125369
                                                                                                                                                    |                                                                                                                                               
 actacATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGt                                                                                                                                                >  2:501311/1‑149 (MQ=255)
  ctacATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTc                                                                                                                                               <  4:1346755/149‑1 (MQ=255)
      aTCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGt                                                                                                                                           >  3:239626/1‑149 (MQ=255)
        cACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTcc                                                                                                                                         >  4:1140891/1‑149 (MQ=255)
        cACTACGCTGAACGAAAACTACAACATGCCGGCAATGCCGCAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTcc                                                                                                                                         >  2:684760/1‑149 (MQ=255)
           tACGCTGAACGAAAACTACCACATGCCGGCAATTCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTg                                                                                                                                      >  3:1609510/1‑149 (MQ=255)
           tACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTg                                                                                                                                      >  3:1357720/1‑149 (MQ=255)
           tACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTg                                                                                                                                      >  3:290204/1‑149 (MQ=255)
                 gAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCTTCACGGACGTGAcgcag                                                                                                                                <  1:1218865/149‑1 (MQ=255)
                  aaCGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCTTCACGGACGTGAcgcagg                                                                                                                               <  3:1691260/149‑2 (MQ=255)
                   aCGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCt                                                                                                                              <  2:1258361/149‑1 (MQ=255)
                   aCGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCt                                                                                                                              <  1:171965/149‑1 (MQ=255)
                     gAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTga                                                                                                                            >  1:910757/1‑149 (MQ=255)
                          cTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATcc                                                                                                                       <  1:1626640/149‑1 (MQ=255)
                            aCCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTg                                                                                                                     <  3:1812221/149‑1 (MQ=255)
                            aCCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTg                                                                                                                     <  1:501311/149‑1 (MQ=255)
                                 aTGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGaa                                                                                                                >  3:411717/1‑149 (MQ=255)
                                  tGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGaaa                                                                                                               <  3:628596/149‑1 (MQ=255)
                                   gCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGaaa                                                                                                               >  3:1122470/1‑148 (MQ=255)
                                    ccGGCAATGCCGGTAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGa                                                                                                             <  4:1609510/149‑1 (MQ=255)
                                    ccGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGa                                                                                                             <  1:1421532/149‑1 (MQ=255)
                                      ggCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTg                                                                                                                     <  4:862187/139‑1 (MQ=255)
                                      ggCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTg                                                                                                                     >  3:862187/1‑139 (MQ=255)
                                      ggCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGAtt                                                                                                           <  1:582592/149‑1 (MQ=255)
                                       gCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGt                                                                                                                                       >  2:640686/1‑120 (MQ=255)
                                       gCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTa                                                                                                          >  4:1728663/1‑149 (MQ=255)
                                         aaTGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACg                                                                                                        <  1:100320/149‑1 (MQ=255)
                                          aTGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACgg                                                                                                       >  1:337132/1‑149 (MQ=255)
                                            gCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTc                                                                                                                                          <  2:1389142/112‑1 (MQ=255)
                                               ggAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTc                                                                                                                                               <  2:1610838/104‑1 (MQ=255)
                                                     tGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTgaga                                                                                                                          >  1:1402522/1‑119 (MQ=255)
                                                     tGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGtt                                                                                               >  1:1514646/1‑146 (MQ=255)
                                                     tGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGtt                                                                                               <  2:1514646/146‑1 (MQ=255)
                                                        tGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACg                                                                                         >  2:57170/1‑149 (MQ=255)
                                                         gAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCTTCACGGACGTGAcgcag                                                                                                                                <  1:636189/109‑1 (MQ=255)
                                                          aGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAACGATTACGGCGTAGGTTCTGTCGtt                                                                                       <  1:1597481/149‑1 (MQ=255)
                                                           ggAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGttt                                                                                      <  4:1419275/149‑1 (MQ=255)
                                                            gAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGa                                                                                           <  3:1480477/143‑1 (MQ=255)
                                                               ggTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTg                                                                                            >  3:1767814/1‑139 (MQ=255)
                                                                gTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGc                                                                                 <  4:239626/149‑1 (MQ=255)
                                                                 tATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCTTCACGGACGTGAcgca                                                                                                                                 <  3:657304/100‑1 (MQ=255)
                                                                  aTCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCt                                                                                             <  2:1679957/135‑1 (MQ=255)
                                                                   tCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGc                                                                                 >  2:396736/1‑146 (MQ=255)
                                                                     cGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTg                                                                                                                     >  2:916113/1‑108 (MQ=255)
                                                                      gTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGAtt                                                                                                           <  3:1378103/117‑1 (MQ=255)
                                                                      gTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGAtt                                                                                                           >  4:1378103/1‑117 (MQ=255)
                                                                      gTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGAcc                                                                           >  2:941344/1‑149 (MQ=255)
                                                                       tAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCTTCACGGACGTGAcgcag                                                                                                                                >  1:398456/1‑95 (MQ=255)
                                                                             tATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCa                                                                    >  4:184403/1‑149 (MQ=255)
                                                                              aTCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCt                                                                                                                              >  3:806539/1‑90 (MQ=255)
                                                                              aTCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCt                                                                                                                              <  4:806539/90‑1 (MQ=255)
                                                                              aTCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGa                                                                             >  4:1119459/1‑139 (MQ=255)
                                                                              aTCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGATCTGACGTTTATAGCGTGa                                                                             <  3:1119459/139‑1 (MQ=255)
                                                                                 tACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGa                                                                >  3:1024803/1‑149 (MQ=255)
                                                                                   cAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGaaa                                                                                                               <  1:121661/100‑1 (MQ=255)
                                                                                    aaaCTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCt                                                                                             <  2:1173464/117‑1 (MQ=255)
                                                                                     aaCTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAagca                                                                                                                                 >  1:1856756/1‑80 (MQ=255)
                                                                                      aCTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGaaa                                                                                                               <  2:1827327/97‑1 (MQ=255)
                                                                                        tCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTGACCTGCGTCACGTCCGTGAAGCAGCt                                                                                                                              >  4:1305949/1‑80 (MQ=255)
                                                                                        tCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCt                                                                                                                              <  3:1305949/80‑1 (MQ=255)
                                                                                        tCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGa                                                                             <  1:929052/129‑1 (MQ=255)
                                                                                        tCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGgcg                                                         >  3:575329/1‑149 (MQ=255)
                                                                                        tCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGgcg                                                         >  3:917754/1‑149 (MQ=255)
                                                                                          gAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGAcctcc                                                                        >  1:942392/1‑132 (MQ=255)
                                                                                          gAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCa                                                                    <  2:1104843/136‑1 (MQ=255)
                                                                                            aaCTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTg                                                     >  1:245225/1‑149 (MQ=255)
                                                                                            aaCTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTg                                                     >  1:1261001/1‑149 (MQ=255)
                                                                                               tATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCTTCACGGACGTGacgc                                                                                                                                  >  4:264274/1‑66 (MQ=38)
                                                                                               tATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATg                                                  >  2:1885692/1‑149 (MQ=255)
                                                                                                aTTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATgg                                                 >  1:40718/1‑149 (MQ=255)
                                                                                                  tGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTc                                               <  3:184403/149‑1 (MQ=255)
                                                                                                     aGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGa                                                    >  4:1524647/1‑141 (MQ=255)
                                                                                                        tAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGAtt                                         <  4:181834/149‑1 (MQ=255)
                                                                                                          gCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACgg                                                                                                       <  3:964055/85‑1 (MQ=255)
                                                                                                          gCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACgg                                                                                                       >  2:1023841/1‑85 (MQ=255)
                                                                                                          gCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGc                                                                                 <  1:1137919/107‑1 (MQ=255)
                                                                                                          gCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATtgt                                       <  4:1024803/149‑1 (MQ=255)
                                                                                                          gCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATtg                                        <  4:337980/148‑1 (MQ=255)
                                                                                                          gCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATtg                                        >  3:337980/1‑148 (MQ=255)
                                                                                                           cAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATtgtg                                      <  2:337132/149‑1 (MQ=255)
                                                                                                             aaGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGaa                                    <  2:1768544/149‑1 (MQ=255)
                                                                                                             aaGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGaa                                    >  1:653740/1‑149 (MQ=255)
                                                                                                                gTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGc                                                                                 <  1:773364/101‑1 (MQ=255)
                                                                                                                gTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGc                                                                                 >  2:773364/1‑101 (MQ=255)
                                                                                                                         aGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACAtg                        <  2:1415559/149‑1 (MQ=255)
                                                                                                                            tgctCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATtgt                                       <  4:177211/131‑1 (MQ=255)
                                                                                                                                    tCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCt             >  1:1055031/1‑149 (MQ=255)
                                                                                                                                    tCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCt             >  2:890439/1‑149 (MQ=255)
                                                                                                                                    tCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCt             >  2:444411/1‑149 (MQ=255)
                                                                                                                                     ccGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTg            <  1:708711/149‑1 (MQ=255)
                                                                                                                                      cGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTgg           <  4:575329/149‑1 (MQ=255)
                                                                                                                                       ggTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGaa                                    >  4:1251540/1‑123 (MQ=255)
                                                                                                                                       ggTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGaa                                    <  3:1251540/123‑1 (MQ=255)
                                                                                                                                         ttCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATtgt                                       >  1:1825367/1‑118 (MQ=255)
                                                                                                                                         ttCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGaaa        <  4:41318/149‑1 (MQ=255)
                                                                                                                                            tATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATtgt                                       <  4:950061/115‑1 (MQ=255)
                                                                                                                                            tATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTg                               >  4:636411/1‑123 (MQ=255)
                                                                                                                                             cagcTGCTTCACGGACGTGACGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGAcctcc                                                                        >  2:1218714/4‑81 (MQ=37)
                                                                                                                                               ccTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCcgc  >  2:609141/1‑149 (MQ=255)
                                                                                                                                                cTGCATCACGGACGTGACGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGaaa        <  2:1445933/142‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGC  >  NC_000913/125078‑125369

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: