Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 126,126 G→A 6.8% K144K (AAG→AAA aceF → pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913126,1260GA 6.8% 391.6 / 6.3 129K144K (AAG→AAAaceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base G (47/72);  new base A (6/3);  total (53/76)
Fisher's exact test for biased strand distribution p-value = 1.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 133 total aligned reads displayed.

GCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAAGC  >  NC_000913/125977‑126275
                                                                                                                                                     |                                                                                                                                                     
gcTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACCTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAg                                                                                                                                                       <  2:81095/150‑1 (MQ=255)
 cTCCGGCAGCAGCACCAGCGGTTGCGGCGGCAAAAGTCGTTAACGTTCCGGATATCGGCAGCGACGAAGGTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCGGAAGAGGCGCTGATCACCGTAGAAGGCGACAAgg                                                                                                                                                      <  1:1101211/150‑1 (MQ=25)
     ggCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTc                                                                                                                                                  >  2:1773907/1‑150 (MQ=255)
     ggCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTc                                                                                                                                                  >  2:1543500/1‑150 (MQ=255)
         gcagcaCCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATg                                                                                                                                              <  2:1901605/150‑1 (MQ=255)
          cagcaCCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATgg                                                                                                                                             >  2:903080/1‑150 (MQ=255)
            gcaCCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGaa                                                                                                                                           >  2:591454/1‑150 (MQ=255)
                cAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTc                                                                                                                                       <  2:1887211/150‑1 (MQ=255)
                cAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTc                                                                                                                                       >  2:698084/1‑150 (MQ=255)
                    gggtGCGGCGGCAAAAGACGGTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGc                                                                                                                                   <  2:100105/147‑1 (MQ=255)
                        gcggcgGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCg                                                                                                                               >  1:662399/1‑150 (MQ=255)
                          ggcggcAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGtt                                                                                                                             <  1:361613/150‑1 (MQ=255)
                          ggcggcAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGtt                                                                                                                             <  1:1876768/150‑1 (MQ=255)
                               cAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTg                                                                                                                        <  1:1788178/150‑1 (MQ=255)
                                      cGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGt                                                                                                                 >  1:791125/1‑150 (MQ=255)
                                      cGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCCTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGt                                                                                                                 >  2:364377/1‑150 (MQ=255)
                                       gTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTg                                                                                                                <  2:1665668/150‑1 (MQ=255)
                                           aCGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAag                                                                                                            <  1:639492/150‑1 (MQ=255)
                                           aCGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAag                                                                                                            <  1:1501713/150‑1 (MQ=255)
                                             gTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAagag                                                                                                          >  2:1050116/1‑150 (MQ=255)
                                               tCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAATGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGAt                                                                                                        >  1:1567960/1‑150 (MQ=255)
                                                ccGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTTGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGTCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATc                                                                                                       <  2:16004/150‑1 (MQ=255)
                                                ccGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTATGCTGGCACCGTGAAAGAGATc                                                                                                       >  2:5209/1‑150 (MQ=255)
                                                 cGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCa                                                                                                      >  1:1937683/1‑150 (MQ=255)
                                                      atCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAg                                                                                                                                          <  2:1583932/109‑1 (MQ=255)
                                                      atCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAg                                                                                                                                          >  1:1583932/1‑109 (MQ=255)
                                                       tCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGa                                                                                                <  1:763370/150‑1 (MQ=255)
                                                          acAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACg                                                                                             <  2:1491941/149‑1 (MQ=255)
                                                            aGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTg                                                                                           <  2:1140245/150‑1 (MQ=255)
                                                             gCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCTCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTg                                                                                                                <  2:946556/128‑1 (MQ=14)
                                                              cgacgaAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTggg                                                                                         <  1:903080/150‑1 (MQ=255)
                                                                acgaAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTg                                                                                       <  1:431797/150‑1 (MQ=255)
                                                                    aGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACaa                                                                                   <  1:1165948/150‑1 (MQ=255)
                                                                      ttGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAg                                                                                 <  1:1078463/150‑1 (MQ=255)
                                                                        gAAGTGACCGAAATCCTGGTGAAAGTTTGCGATAAAGTTGAAGCTGAAGAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTGGCTGGCACCGTGAAAGAGAGCAAAGTGAACGTGGGTGACAAAgtg                                                                               <  1:297520/150‑1 (MQ=255)
                                                                        gAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGaa                                                                                                                                           >  2:1885836/1‑90 (MQ=255)
                                                                        gAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGaa                                                                                                                                           <  1:1885836/90‑1 (MQ=255)
                                                                        gAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGGTGAAGCTCAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAgtg                                                                               <  2:263886/150‑1 (MQ=255)
                                                                              aCCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTAcc                                                                         <  2:1709955/150‑1 (MQ=255)
                                                                               ccGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCg                                                                        >  2:1305869/1‑150 (MQ=255)
                                                                                  aaaTCCTGGCGTACGTGCGCGATATATTCTAAGCAGAACAGTCTTTGATCACCGTAGATGTTGAGATGTCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCt                                                                     <  2:2002196/150‑1 (MQ=255)
                                                                                       cTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTg                                                                >  1:561569/1‑150 (MQ=255)
                                                                                            gAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTAt                                                           >  1:982357/1‑150 (MQ=255)
                                                                                            gAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTAt                                                           <  2:1833266/150‑1 (MQ=255)
                                                                                              aaGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGATGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGCTCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATgg                                                         >  2:304401/1‑150 (MQ=255)
                                                                                               aGTTGGCGATAAAGTTGAAGCTGATGAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCTCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGt                                                        <  1:838370/150‑1 (MQ=255)
                                                                                                gtgggCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGAt                                                                                                        >  2:1284922/4‑101 (MQ=14)
                                                                                                  tGGCGATAAAGTTGACGCTGAACAGTCGCTGATCACCGTAGAAGGCGCCACGGCTTCTAGGGAAGTTCCGGCGCCGTTTGCAGGCGCCGGGAATGAGCTCAAAGTAAACGTTGGGGAAAAAGTGTCTACAGGCGCGCTGATTATGATCta                                                     >  1:1344871/1‑149 (MQ=11)
                                                                                                    gCGATAAACTTGCAGCTGAACAGTCGCTGAGCACCGTAGAATGCGACAATGATTCTATGGAAGTTCAGGCTCCGATTGCTGTCACCGTGAAAGAGATCAAAGTGAACGTGTGTGACAAAGTGTCTACCGGCGCTCTGATTATGGTCTTCg                                                   <  2:340548/150‑1 (MQ=255)
                                                                                                     cGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGa                                                  >  2:663996/1‑150 (MQ=255)
                                                                                                     cGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGa                                                  <  2:1430881/150‑1 (MQ=255)
                                                                                                      gATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACAGTCTCGCTGATTATGGTCTTCGaa                                                 >  2:107432/1‑150 (MQ=255)
                                                                                                       aTAAAGTTGAAGCTGAACAGTCGCTGATCCCCGTAGAAGGCGACAAGGCTTCTAAGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTTACAAAGTGTCTACCGGCTCGGTGATTATGGTCTTCGAAg                                                >  1:1243346/1‑150 (MQ=255)
                                                                                                       aTAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAg                                                <  1:1765829/150‑1 (MQ=255)
                                                                                                           aGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTcgc                                            >  1:1441945/1‑150 (MQ=255)
                                                                                                              tGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCggg                                         <  2:1059866/150‑1 (MQ=255)
                                                                                                              tGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCggg                                         <  2:1518919/150‑1 (MQ=255)
                                                                                                              ggAAGCTGAACAGGCGCTGAGCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCggg                                         <  1:1250824/149‑1 (MQ=255)
                                                                                                                   cTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAg                                    <  2:1187049/150‑1 (MQ=255)
                                                                                                                   cTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAg                                    >  1:1365357/1‑150 (MQ=255)
                                                                                                                   cTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAg                                    <  1:1729599/150‑1 (MQ=255)
                                                                                                                    tGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGGCGTGAAGGAGCTCAAAGTGAACGTGGGCGAGAAAGTGTCTACGGGCTCGCTGATTATGTTCTTCGAAGTGGCGGGCGCAGc                                   >  1:850926/1‑150 (MQ=14)
                                                                                                                     gAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCa                                  <  1:753048/150‑1 (MQ=255)
                                                                                                                     gAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCa                                  <  1:727601/150‑1 (MQ=255)
                                                                                                                     gAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGGTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCa                                  <  1:1773907/150‑1 (MQ=255)
                                                                                                                     gAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCGTCTATGGAAGTTCCGGCTCCGAATGCTGGCACCGTGAAAGAGATCAAAGTGAAAGTGGGTGACAAAGACTCTACCGGCTCGCTGATTATGGTGTTCGAAGACGCCGGAGAAGCa                                  >  2:747129/1‑150 (MQ=21)
                                                                                                                      aaCAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAg                                 <  1:1543500/150‑1 (MQ=255)
                                                                                                                        cAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGt                                                                                                  <  1:1956429/83‑1 (MQ=255)
                                                                                                                        cAGTCGCTGATCACCGTAGAAGGCGACAAAGCATCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGt                                                                                                  >  2:1956429/1‑83 (MQ=255)
                                                                                                                         aGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGgcg                              >  2:1001682/1‑150 (MQ=255)
                                                                                                                         aGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGGCTTCGAAGTCGCGGGTGAAGCAGgcg                              <  2:544631/150‑1 (MQ=255)
                                                                                                                          gTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGgcgc                             >  1:1012311/1‑150 (MQ=255)
                                                                                                                          gTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGgcgc                             <  2:1915281/150‑1 (MQ=255)
                                                                                                                           gcGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGcgcg                            <  2:1989662/149‑1 (MQ=255)
                                                                                                                            cGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCg                                                                                                                               <  1:48791/50‑1 (MQ=255)
                                                                                                                            cGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCg                                                                                                                               <  1:1892026/50‑1 (MQ=255)
                                                                                                                            cGCTGATCACCGTAGAAGGCGACAAGGCGTCTATGGAAGTTCCGGCTCCg                                                                                                                               >  2:1892026/1‑50 (MQ=255)
                                                                                                                            cGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCTCCGGTGGCGGGGAACGTGAAAGAGATCAAAGTGACGGTGGGGGAGAAAGAGGCGACCGGCTCTCTGATTATGGTCTTCGAATACGCCGCAGAAGCAGgcgccg                           >  1:536884/1‑148 (MQ=9)
                                                                                                                             gCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGc                          <  2:1854328/150‑1 (MQ=255)
                                                                                                                             gCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGc                          <  2:2001046/150‑1 (MQ=255)
                                                                                                                               tGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGAGCAAAGTGAACGTGGGGGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAg                        <  1:1645050/150‑1 (MQ=255)
                                                                                                                                gATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGc                       <  1:747129/150‑1 (MQ=255)
                                                                                                                                    aCCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCg                   <  1:856317/150‑1 (MQ=255)
                                                                                                                                     ccGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCgg                  >  2:840183/1‑150 (MQ=255)
                                                                                                                                     ccGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCgg                  <  2:536884/150‑1 (MQ=255)
                                                                                                                                     ccGTAGAAGGCGACAAGGCATCTATGGAAGTTCCGGCTCCGGTGGCTGGCAACGTTAAAGAGATCAAAGTGTACGTTGGAGACAAAGTGTCTACCGGCTCGCTGATTATGGTGTTCGAATTCGCGGCGGAAGCAGGAGCCGCAGCTCCgg                  >  2:1218085/1‑150 (MQ=255)
                                                                                                                                      cGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGc                 <  2:1567960/150‑1 (MQ=255)
                                                                                                                                       gTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTATGCTGCCACCGTGAAAGAGAGCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGcc                <  1:166846/150‑1 (MQ=255)
                                                                                                                                       gTAGAACGCGACAAGGCTTCTATGGAAGTTCCGGCTGCGTTTGCTGGCACCGTTAAAGACATCAAAGTCAACGTGGGTGACATAGTGTCTACCGGCTCGCTGATTATTGTCTTCGAAGTCGCTGGTGAAGCAGGCGCGGCAGCTCCGGcc                <  2:263269/150‑1 (MQ=37)
                                                                                                                                           aaGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGAGAAGCAGGCGCGGCAGCTCCGGCCGCTa            <  2:955521/150‑1 (MQ=255)
                                                                                                                                              gCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGAGCAAAGTGAACGTGGGTGACAAAGTGTCTACCGTCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAAc         <  1:1321760/150‑1 (MQ=255)
                                                                                                                                              gCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGAGCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAAc         <  1:525356/150‑1 (MQ=255)
                                                                                                                                               cGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACa        <  1:1439831/150‑1 (MQ=255)
                                                                                                                                                gACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAg       <  1:1218085/150‑1 (MQ=255)
                                                                                                                                                 aCAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAgg      >  2:569825/1‑150 (MQ=255)
                                                                                                                                                  cAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGa     <  2:558232/150‑1 (MQ=255)
                                                                                                                                                    aGGCTTCTATTGTAGTTCCGTCTCCGTTTGCTGGCACCGTTAAAGATATTAAAGTGAACGTGTGTGACAAAGTGTCTACCGTATCGCTTATTATTGTCTTCGAAGTCGCGGGTGATTCTGTCGATGCAGCTCCGGCCGCTAAACAGGAag   <  2:232236/150‑1 (MQ=255)
                                                                                                                                                    aGGCTTCTATGGAAGATCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGCACTTGGGAGACAAATTGACTACCGTCTCGCTTATTATGGTCTTCGAAGTCGCGGTTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAag   <  2:183184/150‑1 (MQ=255)
                                                                                                                                                     ggCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAagc  <  2:84296/150‑1 (MQ=255)
                                                                                                                                                     ggCTTCTATGGAAGGTCCGGCGCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAagc  <  1:199191/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAAGC  >  NC_000913/125977‑126275

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