Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 124,866 A→C 5.0% D617A (GAC→GCC)  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_000913124,8660AC 5.0% 414.8 / 8.6 140D617A (GAC→GCC) aceEpyruvate dehydrogenase, decarboxylase component E1, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base A (78/55);  new base C (3/4);  total (81/59)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

Only 100 of 144 total aligned reads displayed.

GGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAA  >  NC_000913/124717‑125014
                                                                                                                                                     |                                                                                                                                                    
ggAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGa                                                                                                                                                      >  2:58333/1‑150 (MQ=255)
ggAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGa                                                                                                                                                      <  1:1034840/150‑1 (MQ=255)
ggAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGa                                                                                                                                                      <  2:512359/150‑1 (MQ=255)
ggAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGa                                                                                                                                                      >  1:1123538/1‑150 (MQ=255)
  aaGGGATCAACGAGCTGGGCTCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGAGAcc                                                                                                                                                    <  2:671279/150‑1 (MQ=255)
  aaGGGATCAACGAGCGGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCGCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTCTTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGAcc                                                                                                                                                    <  1:975649/150‑1 (MQ=255)
   aGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACCGCACCCACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCa                                                                                                                                                   >  1:1336909/1‑150 (MQ=255)
   aGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGTGATCTGTGCTGGGCGGCTGGCGACCa                                                                                                                                                   >  2:283662/1‑150 (MQ=255)
    gggATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGGTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGAGCTGTGCTGGTCGCCAGCCGCCCAg                                                                                                                                                  <  1:1155275/150‑1 (MQ=255)
     ggATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGCTGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGTCGCCAGCCGCCCAGc                                                                                                                                                 <  1:1528138/150‑1 (MQ=255)
         cAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGTCGCCAGCCGCCCAGCacag                                                                                                                                             <  1:873828/150‑4 (MQ=255)
         cAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACCCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAgc                                                                                                                                             <  2:695/150‑1 (MQ=255)
          aaCGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAgcg                                                                                                                                            >  1:570418/1‑150 (MQ=255)
            cGAGCTGGGCGCAGGCTGGTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGcgcg                                                                                                                                          <  1:1276695/150‑1 (MQ=255)
             gAGCTGGGCGCAGGGTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACATTCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCTATCTGTGCTGGGCGGCTTGCGACCAGCAAGCGCGt                                                                                                                                         <  2:1581524/150‑1 (MQ=255)
                cTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGc                                                                                                                                      <  2:1902154/150‑1 (MQ=255)
                   ggCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGCTCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTc                                                                                                                                   >  2:1333625/1‑150 (MQ=255)
                     cgcAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGACTTCCt                                                                                                                                 >  2:1221408/1‑150 (MQ=255)
                              ttCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATcggcgg                                                                                                                        >  1:1655538/1‑150 (MQ=255)
                               tCCTGTCTGGCAGCGGCGAGCTCTTACAGCACCAACATTCTGCCGATGATCCCGTTCTACATCTATTAGTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGt                                                                                                                       <  1:615871/150‑1 (MQ=255)
                               tCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACCATCTGCCGATGATCCCGTTCTACATGTTTTACTCGATGTTCGGGTTCGAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCTAGCGCGTGGCTTCCTGATCGGCGGt                                                                                                                       <  1:537713/150‑1 (MQ=255)
                               tCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCCGCAAGCGCGTGGCTTCCTGATCGGCGGt                                                                                                                       >  2:355289/1‑150 (MQ=255)
                               tCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGt                                                                                                                       >  2:925209/1‑150 (MQ=255)
                                   ggctggcAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACtt                                                                                                                   >  1:812076/1‑150 (MQ=255)
                                     ctggcAGCGGCGCCCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTcc                                                                                                                 >  1:960940/1‑150 (MQ=255)
                                        gcAGCGGCGACCTCTTACAGCCCCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTCGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGt                                                                                                              <  1:594897/150‑1 (MQ=255)
                                        gcAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCggg                                                                                                              >  2:708254/1‑149 (MQ=255)
                                         cAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTc                                                                                                             >  1:1532894/1‑150 (MQ=255)
                                            cggcgACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTa                                                                                                          <  2:1547182/150‑1 (MQ=255)
                                             ggcgACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGTGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTac                                                                                                         <  2:1410225/150‑1 (MQ=255)
                                             ggcgACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTac                                                                                                         >  2:1378115/1‑150 (MQ=255)
                                             ggcgACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTac                                                                                                         <  1:864520/150‑1 (MQ=255)
                                               cgACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTacca                                                                                                       >  2:1419419/1‑150 (MQ=255)
                                                gACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTaccac                                                                                                      <  2:439470/150‑1 (MQ=255)
                                                 aCCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTaccacc                                                                                                     >  1:1789563/1‑150 (MQ=255)
                                                  ccTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTAGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCAccc                                                                                                    >  1:1184840/1‑150 (MQ=255)
                                                    tctTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCAcccag                                                                                                  >  2:781648/1‑148 (MQ=255)
                                                    tctTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTg                                                                                                  <  2:8764/150‑1 (MQ=255)
                                                    tctTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTg                                                                                                  >  1:772976/1‑150 (MQ=255)
                                                     ctTACAGCTCCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGCCCAGCACGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGa                                                                                                 >  2:118955/1‑150 (MQ=255)
                                                       tACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAAc                                                                                               <  1:652473/150‑1 (MQ=255)
                                                       tACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAAc                                                                                               >  2:1428882/1‑150 (MQ=255)
                                                           gCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCg                                                                                           >  2:1290204/1‑150 (MQ=255)
                                                                        gCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTgcagca                                                                              >  2:1010521/1‑150 (MQ=255)
                                                                        gCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCCTATTGGCGATCTGTGCTGGGCGGCTGGAGACCAGCAAGCGCGTCGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCAgcagca                                                                              >  1:142033/1‑150 (MQ=255)
                                                                            atgatCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGaa                                                                          <  1:1422095/150‑1 (MQ=255)
                                                                             tgatCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAg                                                                         >  1:1110927/1‑150 (MQ=255)
                                                                              tatCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGa                                                                        <  1:457843/149‑1 (MQ=255)
                                                                                tCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATg                                                                      >  1:668544/1‑150 (MQ=255)
                                                                                tCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATg                                                                      >  1:1845746/1‑150 (MQ=255)
                                                                                 cccGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATgg                                                                     <  1:1909492/150‑1 (MQ=255)
                                                                                 cccGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCGGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATgg                                                                     <  1:1240698/150‑1 (MQ=255)
                                                                                  ccGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGt                                                                    >  2:1089348/1‑150 (MQ=255)
                                                                                     ttCTACATCTATTTCTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGGGGCTTCCTGATCGGCGGTACTTCCGGGCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTccc                                                                 <  2:1659976/150‑3 (MQ=255)
                                                                                        tACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGc                                                              >  1:1362068/1‑150 (MQ=255)
                                                                                        tACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGGGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGACCGGCGAAGGTCTGCAGCACGAAGATGGTCACAGc                                                              <  2:1142564/150‑1 (MQ=255)
                                                                                             cTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCCAGAAGATGGTCACAGCCACAt                                                         >  2:1039670/1‑150 (MQ=255)
                                                                                              tATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACAtt                                                        >  1:1049876/1‑150 (MQ=255)
                                                                                              gaTTACTCGATGTTCGGCGTCCAGCGTATTGGCGATCTGTGCTGGGAGGCTGGCGACCATCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGACGTACCCCCCTGAACGGAGAAGGTCTGCAGCACGAAGATGGTCACAGCCACAtt                                                        >  1:933682/2‑150 (MQ=255)
                                                                                                ttACTCGATGTTCTGCTTCCAGCGTATTGGCGATCTGTGATGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACAACCCTGAACTTCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACTTTCa                                                      <  1:630541/150‑1 (MQ=255)
                                                                                                ttACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCa                                                      <  1:849191/150‑1 (MQ=255)
                                                                                                 tACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAg                                                     >  2:55967/1‑150 (MQ=255)
                                                                                                 tACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAg                                                     >  2:894141/1‑150 (MQ=255)
                                                                                                  aCTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCCGt                                                    >  2:1403383/1‑150 (MQ=255)
                                                                                                  aCTCGATGTTCGGCTTCCAGCGTATTGGCGAACTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACCGCCACATTCAGt                                                    >  1:347755/1‑150 (MQ=255)
                                                                                                  aCTCGATGTTCGGCTTCCAGAGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGt                                                    >  1:304449/1‑150 (MQ=255)
                                                                                                   cTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTc                                                   >  2:549109/1‑150 (MQ=255)
                                                                                                     cGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGc                                                 <  1:798964/150‑1 (MQ=255)
                                                                                                     cGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGc                                                 >  1:753994/1‑150 (MQ=255)
                                                                                                      gATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCt                                                >  2:1809877/1‑150 (MQ=255)
                                                                                                           tCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTa                                           >  1:1181870/1‑150 (MQ=255)
                                                                                                           tCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTa                                           >  1:1107474/1‑150 (MQ=255)
                                                                                                              gCTTCCAGCGGATTGGCGAGCTGTGCGGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATcc                                        <  1:52548/150‑1 (MQ=255)
                                                                                                                ttCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCg                                      <  2:960940/150‑1 (MQ=255)
                                                                                                                       gTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTa                               <  2:427284/150‑1 (MQ=255)
                                                                                                                        tATTGGCGATCTGCGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTAt                              >  1:446791/1‑150 (MQ=255)
                                                                                                                         aTTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTAtc                             >  1:1059540/1‑150 (MQ=255)
                                                                                                                          ttGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTAtct                            <  2:1655538/150‑1 (MQ=255)
                                                                                                                               gATCTGTTCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGCAGATGGTCACAGCCCCATTCAGTCGCTGCCTATCCCGCACTGTATCTCTGAc                       >  2:1877416/1‑150 (MQ=255)
                                                                                                                               gATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCCCCTTCAGTCGAGGACTATCCCGAACTGTATCTCGTAc                       >  2:1884722/1‑150 (MQ=255)
                                                                                                                               gATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTAc                       >  2:1070616/1‑150 (MQ=255)
                                                                                                                               gATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTAc                       >  1:1803921/1‑150 (MQ=255)
                                                                                                                                aTCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACg                      <  2:1442752/150‑1 (MQ=255)
                                                                                                                                  gcttgCTGGTCGCCAGCCGCCCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGAc                    >  2:1901568/4‑150 (MQ=255)
                                                                                                                                   tgtgCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGCACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGAcc                   >  2:40900/1‑150 (MQ=255)
                                                                                                                                   ggtgCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTGCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGGCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGAcc                   <  1:150799/149‑1 (MQ=255)
                                                                                                                                      gCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCgg                >  2:1848156/1‑150 (MQ=255)
                                                                                                                                      gCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCgg                <  2:1123538/150‑1 (MQ=255)
                                                                                                                                       ccgcgcggcTGGCGACCAGCAAGCGCGCGGCGCCCCTCCCCCCCGCACTTCCCGCCGCACCACCCTCAACCCCCCAGGGCCGCCCCACGACGATCCTCACAGCCACATTCCGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgc               <  2:1789563/146‑1 (MQ=255)
                                                                                                                                        tGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgct              <  1:1118261/150‑1 (MQ=255)
                                                                                                                                          ggcggcTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgctta            <  2:611502/150‑1 (MQ=255)
                                                                                                                                          ggcggcTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgctta            <  2:1251359/150‑1 (MQ=255)
                                                                                                                                           gcggcTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgcttac           <  2:1450148/150‑1 (MQ=255)
                                                                                                                                           gcggcTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgcttac           >  2:404270/1‑150 (MQ=255)
                                                                                                                                            cggcTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACCTTCAGTCGCTGACTATCCCGAACTGTATCTCGTACGACCCGgcttacg          >  1:347316/1‑150 (MQ=255)
                                                                                                                                            cggcTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgcttacg          >  2:674813/1‑150 (MQ=255)
                                                                                                                                             ggcTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgcttacgc         <  1:1716786/150‑1 (MQ=255)
                                                                                                                                                   cGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTGCCACCCTGATCGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGa   >  1:97798/1‑150 (MQ=255)
                                                                                                                                                   cGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGa   <  1:925209/150‑1 (MQ=255)
                                                                                                                                                    gACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGaa  <  1:58333/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAA  >  NC_000913/124717‑125014

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