Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 125,233 C→A 7.5% V739V (GTC→GTA) ‡ 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,2330CA 7.5% 989.6 / 47.6 278V739V (GTC→GTA) ‡aceEpyruvate dehydrogenase, decarboxylase component E1, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base C (129/128);  new base A (9/12);  total (138/140)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 283 total aligned reads displayed.

CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATA  >  NC_000913/125086‑125381
                                                                                                                                                   |                                                                                                                                                    
cACTACGCTGAACGAAAACTACAACATGCCGGCAATGCCGCAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTcc                                                                                                                                                     >  2:684760/1‑149 (MQ=255)
   tACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTg                                                                                                                                                  >  3:290204/1‑149 (MQ=255)
       cTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAagc                                                                                                                                              <  2:1049503/149‑1 (MQ=255)
           aCGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCt                                                                                                                                          <  1:171965/149‑1 (MQ=255)
           aCGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCt                                                                                                                                          <  2:1258361/149‑1 (MQ=255)
             gAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGCTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCTTCACGGACGTGAcgcag                                                                                                                                            >  2:1613947/1‑145 (MQ=255)
                aaCTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGAt                                                                                                                                     <  4:1829413/149‑1 (MQ=255)
                  cTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATcc                                                                                                                                   <  1:1626640/149‑1 (MQ=255)
                   tACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCt                                                                                                                                  >  3:495222/1‑149 (MQ=255)
                   tACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCt                                                                                                                                  <  1:1464915/149‑1 (MQ=255)
                       acaTGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCg                                                                                                                              >  1:483026/1‑149 (MQ=255)
                          tGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGaaa                                                                                                                           <  3:618755/149‑1 (MQ=255)
                          tGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGaaa                                                                                                                           <  3:1459127/149‑1 (MQ=255)
                           gCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGaaa                                                                                                                           <  4:1122470/148‑1 (MQ=255)
                            ccGGCAATGCCGGTAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGa                                                                                                                         <  4:1609510/149‑1 (MQ=255)
                              ggCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTg                                                                                                                                 <  4:862187/139‑1 (MQ=255)
                               gCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTa                                                                                                                      >  4:1728663/1‑149 (MQ=255)
                               gCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTa                                                                                                                      <  2:483026/149‑1 (MQ=255)
                                 aaTGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAagcagg                                                                                                                                           >  2:621523/1‑125 (MQ=255)
                                 aaTGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACg                                                                                                                    <  1:100320/149‑1 (MQ=255)
                                  aTGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACgg                                                                                                                   >  1:337132/1‑149 (MQ=255)
                                  aTGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACgg                                                                                                                   >  3:1648424/1‑149 (MQ=255)
                                    gCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTc                                                                                                                                                      <  2:1389142/112‑1 (MQ=255)
                                             tGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCTTCACGGACGTGAcgca                                                                                                                                             <  2:360118/112‑1 (MQ=255)
                                             tGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAg                                                                                                              >  2:319493/1‑143 (MQ=255)
                                             tGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGtt                                                                                                           <  2:1514646/146‑1 (MQ=255)
                                              gCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTc                                                                                                          >  1:1592726/1‑146 (MQ=255)
                                                tGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACg                                                                                                     >  2:57170/1‑149 (MQ=255)
                                                 gAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCTTCACGGACGTGAcgcag                                                                                                                                            <  1:636189/109‑1 (MQ=255)
                                                    gAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGa                                                                                                       >  4:1480477/1‑143 (MQ=255)
                                                         tATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCTTCACGGACGTGAcgca                                                                                                                                             >  2:1157648/1‑100 (MQ=255)
                                                          aTCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATc                                                                                                                                    >  3:8786/1‑108 (MQ=255)
                                                          aTCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGt                                                                                           <  2:1865982/149‑1 (MQ=255)
                                                             cGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTg                                                                                                                                 >  2:916113/1‑108 (MQ=255)
                                                             cGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTg                                                                                                                                 <  1:916113/108‑1 (MQ=255)
                                                              gTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGAtt                                                                                                                       >  4:1378103/1‑117 (MQ=255)
                                                              gTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGAcc                                                                                       >  2:941344/1‑149 (MQ=255)
                                                               tAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCTTCACGGACGTGAcgcag                                                                                                                                            >  1:398456/1‑95 (MQ=255)
                                                                   ggTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCtt                                                                                  >  4:662814/1‑149 (MQ=255)
                                                                     tATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCt                                                                                                                                          <  3:134396/91‑1 (MQ=255)
                                                                      aTCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCt                                                                                                                                          <  4:806539/90‑1 (MQ=255)
                                                                         tACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGa                                                                            >  3:181834/1‑149 (MQ=255)
                                                                            aaaCTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCt                                                                                                         >  1:1173464/1‑117 (MQ=255)
                                                                             aaCTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAagca                                                                                                                                             >  1:1856756/1‑80 (MQ=255)
                                                                                tCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTGACCTGCGTCACGTCCGTGAAGCAGCt                                                                                                                                          >  4:1305949/1‑80 (MQ=255)
                                                                                tCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTg                                                                                                                                 <  1:850086/89‑1 (MQ=255)
                                                                                tCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGgcg                                                                     >  4:1793291/1‑149 (MQ=255)
                                                                                  gAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCa                                                                                <  2:1104843/136‑1 (MQ=255)
                                                                                  gAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCa                                                                                >  1:1104843/1‑136 (MQ=255)
                                                                                   aaaCTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGt                                                                  <  3:1728663/149‑1 (MQ=255)
                                                                                    aaCTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTg                                                                 >  1:245225/1‑149 (MQ=255)
                                                                                       tATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATg                                                              <  4:1357720/149‑1 (MQ=255)
                                                                                         ttGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGATATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGc                                                                                             >  4:1002758/1‑116 (MQ=255)
                                                                                         ttGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGt                                                            >  1:1328401/1‑149 (MQ=255)
                                                                                         ttGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGt                                                            <  1:684624/149‑1 (MQ=255)
                                                                                          tGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTtata                                                                                               >  3:835354/1‑113 (MQ=255)
                                                                                                  gCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACgg                                                                                                                   <  3:964055/85‑1 (MQ=255)
                                                                                                  gCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATtgt                                                   <  1:941344/149‑1 (MQ=255)
                                                                                                  gCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATtg                                                    <  4:337980/148‑1 (MQ=255)
                                                                                                  gCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATtg                                                    >  1:1404833/1‑148 (MQ=255)
                                                                                                  gCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATtg                                                    >  1:18201/1‑148 (MQ=255)
                                                                                                  gCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATtg                                                    <  2:1404833/148‑1 (MQ=255)
                                                                                                   cAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCg                                                                                                                                                    >  3:1056833/1‑51 (MQ=255)
                                                                                                   cAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCg                                                                                                                                                    <  2:1759727/51‑1 (MQ=255)
                                                                                                     aaGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATgg                                                             >  4:193599/1‑136 (MQ=255)
                                                                                                     aaGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATgg                                                             <  3:193599/136‑1 (MQ=255)
                                                                                                     aaGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGaa                                                <  3:640371/149‑1 (MQ=255)
                                                                                                        gTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGc                                                                                             <  1:773364/101‑1 (MQ=255)
                                                                                                        gTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGc                                                                                             >  2:773364/1‑101 (MQ=255)
                                                                                                               tCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACa                                      <  2:40718/149‑1 (MQ=255)
                                                                                                                   ctgctCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACa                                      <  4:253952/145‑1 (MQ=255)
                                                                                                                    tgctCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACAtgctg                                 >  3:3621/1‑149 (MQ=255)
                                                                                                                       tCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGa                                                       >  3:1249294/1‑124 (MQ=255)
                                                                                                                       tCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCAc                              >  3:1377742/1‑149 (MQ=255)
                                                                                                                        cGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCAcc                             >  3:697996/1‑149 (MQ=255)
                                                                                                                            tCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCt                         >  1:1055031/1‑149 (MQ=255)
                                                                                                                            tCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCt                         >  1:555058/1‑149 (MQ=255)
                                                                                                                             ccGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTg                        <  1:599218/149‑1 (MQ=255)
                                                                                                                             ccGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTg                        <  1:708711/149‑1 (MQ=255)
                                                                                                                              cGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTgg                       <  4:575329/149‑1 (MQ=255)
                                                                                                                               ggTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGaa                                                >  4:1251540/1‑123 (MQ=255)
                                                                                                                               ggTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGa                      >  3:1778071/1‑149 (MQ=255)
                                                                                                                               ggTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGCTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGa                      >  3:333691/1‑149 (MQ=255)
                                                                                                                                 ttCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATtgt                                                   <  2:1825367/118‑1 (MQ=255)
                                                                                                                                 ttCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATtgt                                                   >  3:1266205/1‑118 (MQ=255)
                                                                                                                                 ttCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGaaa                    <  4:41318/149‑1 (MQ=255)
                                                                                                                                    tATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATtgt                                                   <  4:950061/115‑1 (MQ=255)
                                                                                                                                     cagcTGCTTCACGGACGTGACGCAGCTGAGATGCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGAcctcc                                                                                    <  1:1218714/78‑1 (MQ=37)
                                                                                                                                     cagcTGCTTCACGGACGTGACGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGAcctcc                                                                                    >  2:1218714/4‑81 (MQ=37)
                                                                                                                                      tCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTAATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCt                         >  4:1099220/1‑139 (MQ=255)
                                                                                                                                       ccTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGCTGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCcgc              <  2:653740/149‑1 (MQ=255)
                                                                                                                                        cTGCTTCACGGACGTGACGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCa                                                          <  3:1120312/104‑1 (MQ=255)
                                                                                                                                        cTGCATCACGGACGTGACGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGaaa                    <  2:1445933/142‑1 (MQ=255)
                                                                                                                                           cttCACGGACGTGACGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATtgt                                                   >  2:970523/3‑108 (MQ=255)
                                                                                                                                           cttCACGGACGTGACGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATtgt                                                   <  1:970523/106‑1 (MQ=255)
                                                                                                                                             tCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTg                        <  4:1539349/133‑1 (MQ=255)
                                                                                                                                             tCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTg                        >  3:1539349/1‑133 (MQ=255)
                                                                                                                                             tCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTg                        >  2:86046/1‑133 (MQ=255)
                                                                                                                                               aCGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATg                                                              >  4:303891/1‑93 (MQ=255)
                                                                                                                                                cGTCCGTGAAGCAGCTGAGATCCTGGCGAAGGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATAGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGt     <  1:1220351/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATA  >  NC_000913/125086‑125381

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: