Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 126,075 T→G 11.9% V127V (GTT→GTG aceF → pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913126,0750TG11.9% 192.7 / 8.2 83V127V (GTT→GTGaceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base T (37/34);  new base G (4/7);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 5.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGG  >  NC_000913/125927‑126220
                                                                                                                                                    |                                                                                                                                                 
aCGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTg                                                                                                                                                  <  2:878167/150‑1 (MQ=255)
  ggTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGc                                                                                                                                                <  1:1254634/150‑1 (MQ=255)
  ggTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTAAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGc                                                                                                                                                <  1:969222/150‑1 (MQ=255)
  ggTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGGTGACGTTCCGGATATCGGCGGCGAGGAAGTTGAAGTGACGGACATCCTGGTGAAAGTGGGc                                                                                                                                                >  1:1140956/1‑150 (MQ=255)
        gcagACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATaaa                                                                                                                                          >  1:1108970/1‑150 (MQ=255)
          agACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGt                                                                                                                                        >  2:936691/1‑150 (MQ=255)
          acaagCTGCACCTGCTCATCCACAATACAAGACATAAGCCGCTCCTGGAGCAGTAGTAGATGCTTCGTAGGATAAAGACTTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAATTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGt                                                                                                                                        <  2:210511/146‑1 (MQ=255)
              gCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCCCCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGCAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGaa                                                                                                                                    >  1:980316/1‑150 (MQ=255)
                 gCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAACTCCTGGTGAAAGTTGGCGAAAAAGTTGAAGCt                                                                                                                                 >  2:251070/1‑150 (MQ=255)
                 gCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCt                                                                                                                                 <  2:367442/150‑1 (MQ=255)
                  aaCCTGCTCAGGCAGAAGATTAGATAGAAGCAGCTGCGGCAGCACGAGCAGCGTCTGAGGTGGCAAAAGACCTTAACGTTCCGGACATCGGCAGCGACGAAGTTGAATTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTg                                                                                                                                <  2:1109101/149‑1 (MQ=255)
                     cTGCTCAGGCATCAGAGAAGAAACAAGCAGCTCCGGAATCATCACCAGCGGCTGCGGCGGCAAAAGACGTAAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAAc                                                                                                                             <  1:614293/150‑1 (MQ=255)
                        cTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGCAATCATGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGt                                                                                                                          >  2:105992/1‑150 (MQ=255)
                         tCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTc                                                                                                                         >  2:908321/1‑150 (MQ=255)
                            ggCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGAGGTTAACGTTCCGGATATCGGCAGCGACGAAGTAGAAGTGACCGAAATCCTGGTGATAGTTGGCGAGAATGTTGAAGCTGAACCGTCGCt                                                                                                                      >  2:502761/1‑150 (MQ=255)
                            ggCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCt                                                                                                                      <  1:946906/150‑1 (MQ=255)
                            ggCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCt                                                                                                                      <  2:539595/150‑1 (MQ=255)
                            ggCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTAGAAGCTGCACAGgcgcg                                                                                                                      >  1:159790/1‑149 (MQ=255)
                            ggCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGAACCAGCGGCTGCGGCGGCAAAAAACGTAAACGTTCCGGATATCGGCATCGACTAAGTTGAAGTGAGCGCAAGCGTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGcg                                                                                                                      >  2:190617/1‑149 (MQ=255)
                               agaagagaagAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATAACGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGAt                                                                                                                   >  1:1109944/1‑150 (MQ=255)
                                      aagaaagcAGCAGCTCCGGGATCAGCACCAGCGGCTGCTGCGGCAAAAGACGTTAACGTTCCGGACATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTa                                                                                                            <  1:442945/150‑1 (MQ=38)
                                       agaaagaaGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCTGCGACGAAGTTGATGTGACCGAAATGCTGGTGAAAGTTGGCGATAAAGTTGAGGCTGAACAGTCGCTGATCAGCGTAg                                                                                                           >  2:294379/1‑150 (MQ=25)
                                        gaaagaaGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTACGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCGTGGGGCAAGTAGGCGAGAAAGTTGTAGCTGACCAGACGCTGATCACCGTAGa                                                                                                          >  1:1191165/1‑150 (MQ=255)
                                            gaagCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGGTGAAGTGACGGAAGTCCTGGTGAAAGTTGGCGAGAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGc                                                                                                      >  1:94395/1‑150 (MQ=255)
                                             aagCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCg                                                                                                     <  1:239134/150‑1 (MQ=14)
                                                    tCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGc                                                                                              <  2:487892/150‑1 (MQ=255)
                                                     ccGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGAGAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCt                                                                                             >  1:1197442/1‑150 (MQ=255)
                                                          agcagcaCCAGCGGCGGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTTAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTAt                                                                                        <  1:1072667/150‑1 (MQ=255)
                                                           gcagcaCCAGCGGATGCGGCGGCAAAAGACGTTAAGGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACAGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACTGTCGCTGATCTCCGTAGAAGGGGACAAGGCTTCTATg                                                                                       >  1:52888/1‑150 (MQ=255)
                                                               caCCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAg                                                                                   <  2:553169/150‑1 (MQ=255)
                                                                 ccAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGtt                                                                                 >  1:733127/1‑150 (MQ=255)
                                                                   aGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAg                                                                                   >  1:84847/1‑146 (MQ=255)
                                                                   aGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAg                                                                                   <  2:84847/146‑1 (MQ=255)
                                                                   aGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGGGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTcc                                                                               <  2:180427/150‑1 (MQ=255)
                                                                    gCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCg                                                                              <  1:1311365/150‑1 (MQ=255)
                                                                     cGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCgg                                                                             <  1:1234664/150‑1 (MQ=255)
                                                                          gcggcgGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCg                                                                        >  1:1274822/1‑150 (MQ=255)
                                                                                gcAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCt                                                                  >  1:867035/1‑150 (MQ=255)
                                                                                gcAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCt                                                                  >  1:871128/1‑150 (MQ=255)
                                                                                gcAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACAGAAATCCTGGTGAAAGTTGGCGATAAAGTTGCAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGATCCGGCTCCGTTTGca                                                                  >  2:762722/1‑149 (MQ=255)
                                                                                 cAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTg                                                                 >  1:1264353/1‑150 (MQ=255)
                                                                                  aaacGACGTGAACGTTCCGGATATTGGCAGCGACGAACTTGAAGTTATCTAAATCATTGTGATACTTGGCTATAAATTTGAAGATTAACATTCTCTGATCAACGTAGAAGTAGAGAAGTCTTCTATGGAAGTTCCGGCTCCGTTTGCTgg                                                                <  2:14220/150‑1 (MQ=16)
                                                                                   aaaGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGc                                                               <  2:1108970/150‑1 (MQ=255)
                                                                                       aCGTTAACGTTCCGGATATCGGCAGCGCCGAAGATGAAGAGACCGAAATCCTGGTGAAAGTTGGAGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAATGCTTCTATGGAAGTTACGGCTACGTTTGCTGGCACCg                                                           >  2:77999/1‑150 (MQ=18)
                                                                                        cGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGCAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGt                                                          >  1:635061/1‑150 (MQ=255)
                                                                                               gTTGCGGATATCGGCAGTGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTTAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGTCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAagag                                                   <  2:38331/150‑1 (MQ=255)
                                                                                                ttCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAgaga                                                  <  2:899386/150‑1 (MQ=255)
                                                                                                 cccGGAAATCGCCAGCGACTAAGTTTAATTGTCCGAAAACCTGGTTAAAGTTGTCGATAAACTTGTATCTGACCAGTCGATGATCACCCTAGCAGGCGACAATGCATCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGAt                                                 <  2:391730/149‑1 (MQ=255)
                                                                                                  ccGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATc                                                >  1:417580/1‑150 (MQ=255)
                                                                                                       tatCGGCAGTGTCGAAGGGCAAGGTCCCGATATCCTTGCGAATGTGGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGt                                           <  1:434291/150‑1 (MQ=18)
                                                                                                           ggCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAAc                                       >  2:1010147/1‑150 (MQ=255)
                                                                                                            gCAGCGACGAAGTTGAAGTGTCCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACg                                      >  1:696539/1‑150 (MQ=255)
                                                                                                             gtcggatgaAGTTGAATTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGt                                     <  2:1000965/143‑1 (MQ=255)
                                                                                                                cgacgaAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCAGGCAGCGTGAATGATCTCAAAGTGAACGTTGGt                                 >  2:384682/1‑150 (MQ=18)
                                                                                                                 gacgaAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAACGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGACGTGAATGAGATCAAAGTCAACGTtggg                                 >  2:1072468/1‑149 (MQ=17)
                                                                                                                  acgaAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTg                                <  1:442239/150‑1 (MQ=255)
                                                                                                                    gaAGGTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGCAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGAc                              <  1:460989/150‑1 (MQ=255)
                                                                                                                     aaGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACa                             >  2:863028/1‑150 (MQ=255)
                                                                                                                      aGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGACGTGAATGAACTGAAAGTCAACGTTGGCGAGaa                            >  2:804713/1‑150 (MQ=14)
                                                                                                                        ttGAATTGACCGATATCCTGGTCAAAGTTGTCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAg                          <  1:1072095/150‑1 (MQ=255)
                                                                                                                        ttGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAg                          <  2:393486/150‑1 (MQ=17)
                                                                                                                        ttGAAATCACCGACATCCTGGTCAAAGTTGGCGAAAAAGTTGAAGCGGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAg                          <  2:951306/150‑1 (MQ=18)
                                                                                                                          gAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAgtg                        <  2:1101681/150‑1 (MQ=255)
                                                                                                                          gAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAgtg                        <  1:908321/150‑1 (MQ=255)
                                                                                                                             gTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCt                     <  1:337478/150‑1 (MQ=255)
                                                                                                                             gTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCAGGCGACGTGAAGGAGATCAAAGTCAACGTTGGAGAAAAAgtgaaa                     >  1:327003/1‑147 (MQ=17)
                                                                                                                              tGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTa                    >  1:214064/1‑150 (MQ=255)
                                                                                                                              tGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTa                    >  1:250116/1‑150 (MQ=255)
                                                                                                                                 ccGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTg                                                         <  1:1301545/110‑1 (MQ=14)
                                                                                                                                 ccGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTg                                                         >  2:1301545/1‑110 (MQ=14)
                                                                                                                                 ccGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCTCTGATCACCGTAGAAGGCGACAAAGc                                                                                              >  2:793583/1‑73 (MQ=255)
                                                                                                                                 ccGAAATCCTGGTGAAAGTGGGCGATAAAGTTGAAGCTGAACAGTCTCTGATCACCGTAGAAGGCGACAAAGc                                                                                              <  1:793583/73‑1 (MQ=14)
                                                                                                                                  cGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCgg                >  2:127380/1‑150 (MQ=255)
                                                                                                                                       tCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAgg                                                                <  1:83164/97‑1 (MQ=255)
                                                                                                                                       tCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAgg                                                                >  2:83164/1‑97 (MQ=255)
                                                                                                                                          tGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTg                                                                 >  1:696462/1‑93 (MQ=11)
                                                                                                                                          tGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTg                                                                 <  2:696462/93‑1 (MQ=11)
                                                                                                                                           ggTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTa                    >  2:834752/1‑137 (MQ=255)
                                                                                                                                           ggTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTa                    <  1:834752/137‑1 (MQ=255)
                                                                                                                                            gtcaaaGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACg                                      <  1:386735/115‑1 (MQ=17)
                                                                                                                                            gtcaaaGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACg                                      >  2:386735/4‑118 (MQ=17)
                                                                                                                                                aaGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCTCCGTAGAAGCCTACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGTAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATgg  <  2:52888/150‑1 (MQ=255)
                                                                                                                                                aaGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATgg  <  1:529296/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGG  >  NC_000913/125927‑126220

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: