Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 126,069 G→C 7.1% V125V (GTG→GTC aceF → pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913126,0690GC 7.1% 251.8 / 2.6 87V125V (GTG→GTCaceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base G (44/34);  new base C (2/5);  total (46/41)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGG  >  NC_000913/125920‑126206
                                                                                                                                                     |                                                                                                                                         
tCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCTGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTg                                                                                                                                           <  1:965480/150‑1 (MQ=255)
 cctccgACGGTGCACCAGACGCTGAACATGCTCAGGCAGAAGCGAAGACAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCTGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGa                                                                                                                                          <  2:603370/147‑1 (MQ=38)
 ccgccgACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCAGCGGCGGCAAAAGACGATAACGTACCGGATATCCGCAGCGAAGAAGATGAAGATACGGAACTAATGGTGa                                                                                                                                          >  2:602284/1‑150 (MQ=255)
 ccgccgACGGCGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCTGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGa                                                                                                                                          <  2:597845/150‑1 (MQ=255)
  cgccgACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGaa                                                                                                                                         >  1:1095722/1‑150 (MQ=255)
  cgccgACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGaa                                                                                                                                         >  1:1101681/1‑150 (MQ=255)
  cgccgACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGaa                                                                                                                                         <  1:422322/150‑1 (MQ=255)
   gccgACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGaaa                                                                                                                                        >  1:180427/1‑150 (MQ=255)
    ccgACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGCCGAAGTTGAAGTGACCGAAATCCTGGTGAAAg                                                                                                                                       >  1:899386/1‑150 (MQ=255)
       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                                                                                                        <  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          <  1:908321/150‑1 (MQ=255)
                                                                                                                                 gAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAgtg          <  2:1101681/150‑1 (MQ=255)
                                                                                                                                    gTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCt       <  1:337478/150‑1 (MQ=255)
                                                                                                                                    gTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCAGGCGACGTGAAGGAGATCAAAGTCAACGTTGGAGAAAAAgtgaaa       >  1:327003/1‑147 (MQ=17)
                                                                                                                                     tGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTa      >  1:250116/1‑150 (MQ=255)
                                                                                                                                     tGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTa      >  1:214064/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                                                  >  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)
                                                                                                                                                     |                                                                                                                                         
TCCGCCGACGGTGCAGCAGACGCTGCACCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGG  >  NC_000913/125920‑126206

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: