Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 126,081 T→C 11.4% D129D (GAT→GAC aceF → pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913126,0810TC11.4% 182.3 / 9.8 81D129D (GAT→GACaceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base T (33/33);  new base C (4/5);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAG  >  NC_000913/125935‑126229
                                                                                                                                                  |                                                                                                                                                    
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                                                                                                                               <  2:539595/150‑1 (MQ=255)
                    ggCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCt                                                                                                                               <  1:946906/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                                               >  2:386735/4‑118 (MQ=17)
                                                                                                                                    gtcaaaGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACg                                               <  1:386735/115‑1 (MQ=17)
                                                                                                                                        aaGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCTCCGTAGAAGCCTACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGTAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATgg           <  2:52888/150‑1 (MQ=255)
                                                                                                                                        aaGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATgg           <  1:529296/150‑1 (MQ=255)
                                                                                                                                            gggCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCGTCTATGGAAGTTCCGGCTCCGTTTGCTGGCAACGTGAAAGAGATCAAAGTGAAAGTTGGTGACAAAGTGTCTACCGGCTAGCTGATTATGGTCtt       >  2:784935/2‑150 (MQ=21)
                                                                                                                                                 aTAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGGCTTCGAAg  <  1:304731/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAG  >  NC_000913/125935‑126229

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: