Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 126,093 T→C 8.4% A133A (GCT→GCC aceF → pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913126,0930TC 8.4% 185.5 / 2.6 70A133A (GCT→GCCaceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base T (35/28);  new base C (2/4);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGT  >  NC_000913/125945‑126237
                                                                                                                                                    |                                                                                                                                                
cACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTg                                                                                                                                                 <  1:479325/150‑1 (MQ=255)
  ccTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGaa                                                                                                                                               >  1:585596/1‑150 (MQ=255)
        cAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCg                                                                                                                                         >  2:1290355/1‑150 (MQ=255)
        cAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCg                                                                                                                                         <  2:1180353/150‑1 (MQ=255)
              gaagagaagaAAGAAGCAGCTCCGGCAGCAGCCCCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACCGTCGCTGagc                                                                                                                                   >  2:851598/1‑148 (MQ=255)
              gaagagaagaAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTGCCTGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGATGAACAGTCGCTGCTc                                                                                                                                   >  2:809930/1‑150 (MQ=255)
               aagagaagaAAGAAGCAGCTCCGGCAGCTGCTCCAGCGGCTGCGGCGTCAAATGAGGTTAACTTTCCGGATATCGGCTGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCa                                                                                                                                  <  1:445517/150‑1 (MQ=255)
                   gaagaaAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGt                                                                                                                              <  1:230123/150‑1 (MQ=255)
                    aataaagaaGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTTGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTa                                                                                                                             <  2:746539/147‑1 (MQ=255)
                    aagaaagaAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTa                                                                                                                             <  2:1085930/150‑1 (MQ=255)
                      gaaagaaGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGa                                                                                                                           >  2:766150/1‑150 (MQ=255)
                      gaaagaaGCAGCTCCGGCAGAAGCACCAGCGGCTGCGGCGGCAAAAGACTTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGAACGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCCCCGTAGa                                                                                                                           >  2:20121/1‑150 (MQ=255)
                        aagaagCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAg                                                                                                                         <  2:595670/150‑1 (MQ=255)
                              cagcTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACa                                                                                                                   <  2:533745/150‑1 (MQ=255)
                              cagcTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACa                                                                                                                   <  1:179777/150‑1 (MQ=255)
                                     ggCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTc                                                                                                            >  2:9826/1‑150 (MQ=255)
                                          cagcaCCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCGGCGACGAAGTTGAAGTGACCGAAGTCCTGGTGAAAGTGGGCGAGAAAGTTGCCGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATgg                                                                                                       >  2:539517/1‑150 (MQ=17)
                                           agcaCCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGAGACGAAGTTGAAGTGACAGAAATGATGGTGAAAGTTGGAGAAAAAGTTGCAGCTGAACAGTCACTGATCACCTTAGAAGCCGAAAAAGCTTCTATGGa                                                                                                      >  2:1117064/1‑150 (MQ=14)
                                            gcaCCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGaa                                                                                                     >  1:772281/1‑150 (MQ=255)
                                             caCCAGCGGCGGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGATATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAg                                                                                                    <  2:34467/150‑1 (MQ=255)
                                                cAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTc                                                                                                 <  2:1241558/150‑1 (MQ=255)
                                                    ggCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGc                                                                                             <  2:585596/150‑1 (MQ=255)
                                                    ggCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGCTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGc                                                                                             >  2:885410/1‑150 (MQ=255)
                                                        gcggcgGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAATTCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCg                                                                                         >  1:523074/1‑150 (MQ=255)
                                                         cggcggCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGt                                                                                        <  2:1300589/150‑1 (MQ=255)
                                                                   aGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCAc                                                                              <  1:520130/150‑1 (MQ=255)
                                                                    gACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCAcc                                                                             >  1:790259/1‑150 (MQ=255)
                                                                      cGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAATTCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGACGCTGATCACCGTAGAAGGCGACAAGCCTTCTATGGAAGTTCCGGCTCCGTTTGCAGGctacgg                                                                           >  2:1015613/1‑145 (MQ=17)
                                                                       gTTAACGTTCCGGATATCGGCTGCGACTTAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTg                                                                          <  2:514458/150‑1 (MQ=255)
                                                                       gTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTg                                                                          >  1:1092180/1‑150 (MQ=255)
                                                                       gTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTg                                                                          <  2:1144694/150‑1 (MQ=255)
                                                                          aaCGTTCCGGATATCGGCAGCTACGAATTTGAAGTTACCGTAAACCTGGTGTAAGTAGGCGTTAAAGCTGAAGCTGAACAGTCGCTGATCACCGTATAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGaaa                                                                       <  2:616335/150‑1 (MQ=25)
                                                                          aaCGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGaaa                                                                       >  2:1261815/1‑150 (MQ=255)
                                                                                ccGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATc                                                                 >  2:689385/1‑150 (MQ=255)
                                                                                aagttcatCCGCCGCGAAGAAGTTGAATTCACCGACATCCTGGTCAAAGTTGGCGAAAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCCGCTCCGTTTGCTGGCACCGTGAAAGAGATc                                                                 <  1:531833/144‑1 (MQ=17)
                                                                                     tatCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGt                                                            >  2:855760/1‑150 (MQ=255)
                                                                                       tCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGa                                                          >  2:623431/1‑150 (MQ=255)
                                                                                       tCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGa                                                          >  1:1219230/1‑150 (MQ=255)
                                                                                          gCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCTCCGTTTGCAGGCGACGTGAAAGACCTGAAAGTCAACg                                                       >  2:633056/1‑150 (MQ=14)
                                                                                            cgCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTTGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTg                                                     <  2:668971/149‑1 (MQ=255)
                                                                                                   aaGTTGAAGTGACCGAAATCCTGGTGATAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCTCAGTAGAAGGAGACAAGGCTTCTATGGAAGTTCCGGCTGCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACa                                              <  1:218764/150‑1 (MQ=25)
                                                                                                   aaGTTGAAATCACCGAGATCCTGGTTAAAGTGGGCGACAATGTTGAAGCCGAACAGTCGCTGATCACCGTACAAGGCGACAAGGCATCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACa                                              <  1:237862/150‑1 (MQ=17)
                                                                                                    aGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAg                                                             <  2:1134053/134‑1 (MQ=11)
                                                                                                    aGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAg                                                             >  1:1134053/1‑134 (MQ=11)
                                                                                                     gTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGCCGCTGAACAGTCACTGATCa                                                                                                                                  <  2:865227/64‑1 (MQ=14)
                                                                                                     gTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCACTGATCa                                                                                                                                  >  1:865227/1‑64 (MQ=255)
                                                                                                        gAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAATGAGCTGAAAGTGAACGTGGGTGAAAAAgtg                                         >  1:798034/1‑150 (MQ=255)
                                                                                                          aatGACCGAGATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCCGAACAGTAGCTGATCACCGTAGAAGGCGACAATGATTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTc                                       <  1:882019/148‑1 (MQ=21)
                                                                                                           gagACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTTAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGATCCGTTTGCTGGCGGCGTTATAGATATCAAATTAAACTTGGGAGAAAAAgtgaca                                      >  2:992254/3‑147 (MQ=18)
                                                                                                           gTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCAGGCGACGTGAATGAGCTCAAAGTGAACGTGGGAGAGAAAgtgaca                                      >  2:423705/1‑147 (MQ=18)
                                                                                                             gACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGATCAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTAc                                    >  1:1233041/1‑150 (MQ=255)
                                                                                                              cccGAAATCCTGGTGAAAGTTGGCGAAAAAGTTGAAGCTGAACAGTCGCTGATCACGGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTAcc                                   <  2:457217/149‑1 (MQ=255)
                                                                                                                       cTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTg                          >  1:796141/1‑150 (MQ=255)
                                                                                                                         ggTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGAt                        >  2:837824/1‑150 (MQ=255)
                                                                                                                         ggTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGAt                        >  1:609380/1‑150 (MQ=255)
                                                                                                                         ggTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCgg                                                                                              >  2:354158/1‑80 (MQ=255)
                                                                                                                         ggTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCgg                                                                                              <  1:354158/80‑1 (MQ=255)
                                                                                                                            gAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTAt                     >  1:1321831/1‑150 (MQ=255)
                                                                                                                            caaaGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAACGCGTCTATGGAAGTTCCGTCTCCGGTTGCTGGAAGCGTGAAAGAGATCAAAGTGAAAGTTGGAGACAAAGAGACTACCGGCTCGCTGATTAt                     >  2:1092904/2‑150 (MQ=17)
                                                                                                                              aaGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAATGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGa                                                                         <  2:919100/96‑1 (MQ=11)
                                                                                                                              aaGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGa                                                                         >  1:919100/1‑96 (MQ=11)
                                                                                                                               aGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCATAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTTGGAGAAAAAGTGACTACAGGCTCGCTGATTATGGt                  >  1:876008/1‑150 (MQ=18)
                                                                                                                                gatGGCGATAAAGCTGAAGCGGAACTGTCTCTGATCACCGTTGCAGGCGACTAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCATATTGAACGTGGGTGACAAAGTGTCTACCGGAACGCTGATTATGGTc                 <  1:245983/148‑1 (MQ=25)
                                                                                                                                  tGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCtt               >  2:839301/1‑150 (MQ=255)
                                                                                                                                     cGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGa            >  1:1264658/1‑150 (MQ=255)
                                                                                                                                       aTAAAGTTGAAGCTGATCAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAg          <  1:1290355/150‑1 (MQ=255)
                                                                                                                                            gTTGAAGCTGAACAGTCGCTGATCACCTTAGAACGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTcgcg     <  1:779824/150‑1 (MQ=255)
                                                                                                                                             ttGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCgg    <  1:766150/150‑1 (MQ=255)
                                                                                                                                              tGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCggg   <  1:719070/150‑1 (MQ=255)
                                                                                                                                               gAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCCGATTATGGTCTTCGAAGTCGCGGGt  <  1:9826/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGT  >  NC_000913/125945‑126237

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: