Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 126,093 T→C 15.9% 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,0930TC15.9% 198.1 / 20.6 88A133A (GCT→GCCaceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base T (36/36);  new base C (7/7);  total (43/45)
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

GCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGA  >  NC_000913/125944‑126239
                                                                                                                                                     |                                                                                                                                                  
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                                                                                                      <  2:84847/146‑1 (MQ=255)
                                                  aGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAg                                                                                                      >  1:84847/1‑146 (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                                                                            >  2:1301545/1‑110 (MQ=14)
                                                                                                                ccGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTg                                                                            <  1:1301545/110‑1 (MQ=14)
                                                                                                                ccGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCTCTGATCACCGTAGAAGGCGACAAAGc                                                                                                                 >  2:793583/1‑73 (MQ=255)
                                                                                                                ccGAAATCCTGGTGAAAGTGGGCGATAAAGTTGAAGCTGAACAGTCTCTGATCACCGTAGAAGGCGACAAAGc                                                                                                                 <  1:793583/73‑1 (MQ=14)
                                                                                                                 cGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCgg                                   >  2:127380/1‑150 (MQ=255)
                                                                                                                      tCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAgg                                                                                   >  2:83164/1‑97 (MQ=255)
                                                                                                                      tCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAgg                                                                                   <  1:83164/97‑1 (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)
                                                                                                                                          aaaGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTc          <  1:1005351/150‑1 (MQ=255)
                                                                                                                                          aaaGTTGAAGCGGAACAGTCGCTGATGACCGTAGACGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTc          <  1:784935/150‑1 (MQ=255)
                                                                                                                                          aaaGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTgg                                                                                   <  2:625057/77‑1 (MQ=255)
                                                                                                                                          aaaGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTgg                                                                                   >  1:625057/1‑77 (MQ=255)
                                                                                                                                           aaGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCTCCGTGAATGAGCTCAAAGTGAACGTGGGTGACAAAGTGTATACCGGCTCGCTGATTATGATCTTCGAAGtag         >  2:1179146/1‑148 (MQ=18)
                                                                                                                                             gTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGTCGACAACGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCAACGTGAAAGAGCTCACAGTGAACGTGGGTGACAAAGTGGCTACCGGCTCGCTGATTTAGTTCTTCGAAGtagac       >  2:990236/1‑146 (MQ=25)
                                                                                                                                               tGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCggg     <  2:250116/150‑1 (MQ=255)
                                                                                                                                               tGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGGCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCggg     <  2:214064/150‑1 (MQ=255)
                                                                                                                                                  aGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGa  <  1:886078/150‑1 (MQ=255)
                                                                                                                                                  aGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAg                                                               >  1:444566/1‑89 (MQ=255)
                                                                                                                                                  aGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAg                                                               <  2:444566/89‑1 (MQ=255)
                                                                                                                                                  aGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTggg                                                     <  2:533835/99‑1 (MQ=255)
                                                                                                                                                  aGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTggg                                                     >  1:533835/1‑99 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGA  >  NC_000913/125944‑126239

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: