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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913126,0930TC 7.2% 223.5 / 7.6 83A133A (GCT→GCCaceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base T (29/48);  new base C (3/3);  total (32/51)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

CCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAA  >  NC_000913/125947‑126240
                                                                                                                                                  |                                                                                                                                                   
ccTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCTGCAGCAGCACCAGCTGCTGCGGCGGCAAAAGACGTTAACGTTAAGGATATCGGCAGAGACGAAGTTGAAGTGCCCGAAATCCTGGTTATAGTTTGAGAATAAGTTGGTGCTGa                                                                                                                                                   >  1:1429956/1‑149 (MQ=255)
 cTGCTCAGGCAGAAGAGAAGACAGCAGCCGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGaa                                                                                                                                                  <  2:542420/149‑1 (MQ=255)
   gCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACa                                                                                                                                                <  1:1213011/149‑1 (MQ=255)
   gCTCAGGAAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGAGGCAAAAGACGTTAAAGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGAGATAAAGTTGAAGCTGAACa                                                                                                                                                <  2:300286/149‑1 (MQ=255)
       aGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAATTTGAAGCTGAACAGTCg                                                                                                                                            >  1:1074230/1‑149 (MQ=255)
        ggCAGAAGAGAAGAAAGAAGCAGCTCCGACAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGGTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGc                                                                                                                                           <  2:1370443/149‑1 (MQ=255)
          cAGAAGAGAAGAAAGAAGCAGCTCCGTCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTg                                                                                                                                         >  2:1428941/1‑149 (MQ=255)
            gtagagaagaAAGAAGCAGCTCGGGCAGCAGCACGAGCGGCTGCGGCGGGAAAAGACGTTAAGGTTCCGGATATCGGCAGCGACGAAGTTGAAGTTACCGAATTCCTGGTGAATGTTGGCGATAAAGTTGATGCTGAAAAGTCGCTGAt                                                                                                                                       <  2:487870/147‑1 (MQ=255)
             aagagaagaAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATc                                                                                                                                      >  2:498716/1‑149 (MQ=255)
                agaagaAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCAcc                                                                                                                                   <  2:998445/149‑1 (MQ=255)
                    gaaagaaGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGAAg                                                                                                                               >  2:346473/1‑149 (MQ=255)
                     aaagaaGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGa                                                                                                                              >  1:197108/1‑149 (MQ=255)
                           gcagcTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTTGTGAAAGTTGGCGATAAAGTTGCAGCTGAACAGTCGCTGATCACCGTAGAAGGCGa                                                                                                                        <  2:145429/149‑1 (MQ=255)
                           gcagcTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGa                                                                                                                        <  1:877324/149‑1 (MQ=255)
                               cTCCTGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTTGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAg                                                                                                                    <  1:11220/149‑1 (MQ=255)
                                tCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAgg                                                                                                                   >  1:65231/1‑149 (MQ=255)
                                  cGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGCCCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGAGCACCGTAGAAGGCGACAAGGCt                                                                                                                 >  2:1294481/1‑149 (MQ=255)
                                  cGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCt                                                                                                                 >  1:1246244/1‑149 (MQ=255)
                                      agcagcaCCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTa                                                                                                             >  1:1126016/1‑149 (MQ=255)
                                         atcaCCAGCGGCTGCTGCTGCAAAAGACGTTAACGTTCCGGATATCGGCATCGACGAAGTTGAAGTGATCGATATCCTTGTGAAAGTTGGCTATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATgg                                                                                                          <  2:307314/147‑1 (MQ=21)
                                         agcaCCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATgg                                                                                                          <  1:583996/149‑1 (MQ=255)
                                          gcaCCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGATAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGa                                                                                                         >  1:1363996/1‑149 (MQ=255)
                                          gcaCCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGa                                                                                                         <  2:119526/149‑1 (MQ=255)
                                          gcaCCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGGGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGACGCTGATCACCGTAGACGGCGAAAAGGCTTcgagtga                                                                                                         >  2:853188/1‑143 (MQ=25)
                                              cAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGCCCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGtt                                                                                                     >  1:386524/1‑149 (MQ=255)
                                               aGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTc                                                                                                    <  1:14707/149‑1 (MQ=255)
                                                gCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGAGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTcc                                                                                                   >  1:338936/1‑149 (MQ=255)
                                                  ggCTGCGGCGGCAAAAGAGGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGGTAAAGTTGAAGCTGAACAGGCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCgg                                                                                                 <  1:872480/149‑1 (MQ=255)
                                                  ggCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCgg                                                                                                 <  2:165274/149‑1 (MQ=255)
                                                      gcggcgGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTcc                                                                                             <  1:324090/149‑1 (MQ=255)
                                                            gcAAAATACCTTAACGTACCGCAAATCGGCAGCGACGAATTTGATTTTTATTAAATCCTGGTGAAATTTTTCGATAATGTTGAATCTTAACATTCGCTGATCACGGTAGATGGCGACAATGCTTCTATGGAAGTTCCGGCTCCGTTTGc                                                                                       <  2:1518707/149‑1 (MQ=255)
                                                              aaaaGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTg                                                                                     <  1:643056/149‑1 (MQ=255)
                                                              aaaaGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGAAGAAGGCGACATTGCTTCTATGGAAGTTCAGGCTCCGGTTGCTg                                                                                     >  2:1445402/1‑149 (MQ=18)
                                                               aaaGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGTCCGAAATCCTGGTGTAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTTGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTgg                                                                                    <  2:392072/149‑1 (MQ=255)
                                                                 aGACGCTAACGTTCCGGATCTCTGAAGCGACGAAGTTGAAGTGACCGAACTCCGGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCa                                                                                  <  2:1024390/149‑1 (MQ=18)
                                                                   aCGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCAcc                                                                                <  1:498716/149‑1 (MQ=255)
                                                                         aCGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCTTGaaa                                                                          <  1:702113/149‑1 (MQ=255)
                                                                            ttCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAagag                                                                       <  1:1247458/149‑1 (MQ=255)
                                                                             tCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCATTTGCTGGCACCGTGAAAgaga                                                                      <  1:1428941/149‑1 (MQ=255)
                                                                                    atCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCAGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTTGCACCGTGAAAGAGATCAAAGt                                                               >  1:849241/1‑149 (MQ=255)
                                                                                       ggCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGCACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGaa                                                            >  1:892877/1‑149 (MQ=255)
                                                                                        gCAGCGACGAAGCAGACGTGACCTAAATCCTGGTGTAATTTGGCGATAAAGTTGATGCTGAACAGTCGCTTATCACCGTAGAAGCCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGTAc                                                           <  2:597903/149‑1 (MQ=255)
                                                                                            cttcgtaGTTGAATTGACCGTTATCGTGGTGTAAGTTGGCTATAATGTTTAAGCTGAACTGTGGCTTTTCAGCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTgg                                                       <  1:145781/143‑1 (MQ=255)
                                                                                              acgaAGTTGAAGTGACCGAAATCCTGGGGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGt                                                     <  1:616373/149‑1 (MQ=255)
                                                                                                    ttGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACaaa                                               >  1:142090/1‑149 (MQ=255)
                                                                                                      gAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAgt                                             <  1:303602/149‑1 (MQ=255)
                                                                                                            aCCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACTTGGGTGACAAAGTGTCTAc                                       <  2:65231/149‑1 (MQ=255)
                                                                                                            aCCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTAc                                       <  1:1382288/149‑1 (MQ=255)
                                                                                                              cGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCTCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGcc                                                                                             <  1:1114666/93‑1 (MQ=255)
                                                                                                              cGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGcc                                                                                             >  2:1114666/1‑93 (MQ=255)
                                                                                                              cGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGGCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCg                                     <  1:820349/149‑1 (MQ=255)
                                                                                                                aaaTCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCTCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGc                                                                                   <  1:293775/101‑1 (MQ=14)
                                                                                                                aaaTCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGc                                   <  2:1239465/149‑1 (MQ=255)
                                                                                                                aaaTCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGc                                   >  2:17292/1‑149 (MQ=255)
                                                                                                                aaaTCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGc                                                                                   >  2:293775/1‑101 (MQ=255)
                                                                                                                aaaTCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCGTCTATGGAAGTTCCGGCTCCGGTTGCGGGCAGCGTGAAAGAGATCAAAGTGCACGTGGGTGACAAAGTGTCTACCGGc                                   <  2:422541/149‑1 (MQ=21)
                                                                                                                     cTGGGGAAAGTTGGCGATAAAGTTGCAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCATAGTGAACGTTGGTGACAAAGTGTCTACCGGCTCGCt                              <  2:165943/149‑1 (MQ=255)
                                                                                                                         tcaaaGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTgg                                                       <  2:427984/118‑1 (MQ=255)
                                                                                                                         tcaaaGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTgg                                                       >  1:427984/3‑120 (MQ=17)
                                                                                                                            acgTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGCCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATg                       <  1:937386/147‑1 (MQ=255)
                                                                                                                            aaGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGaaa                                                                          <  2:226374/98‑1 (MQ=255)
                                                                                                                            aaGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCATCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGaaa                                                                          >  1:226374/1‑98 (MQ=14)
                                                                                                                                 ggCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGaa                                                                                                        >  1:404833/1‑63 (MQ=255)
                                                                                                                                 ggCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGaa                                                                                                        <  2:404833/63‑1 (MQ=255)
                                                                                                                                 ggCGATAAAGTTGAAGCTGAACAGTCGCGGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCtt                  <  2:849241/149‑1 (MQ=255)
                                                                                                                                  gCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACAGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTc                 >  1:5746/1‑149 (MQ=255)
                                                                                                                                  gCGATAAAGGTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTc                 <  1:85374/149‑1 (MQ=255)
                                                                                                                                    gATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGa               >  2:1322550/1‑149 (MQ=255)
                                                                                                                                    gATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGa               >  1:1444524/1‑149 (MQ=255)
                                                                                                                                     aTAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGaa              <  2:1246244/149‑1 (MQ=255)
                                                                                                                                     aTAAAGTTGAAGCTGAACAGTCGCTGATCACCCTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCGGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGaa              <  2:792992/149‑1 (MQ=255)
                                                                                                                                       aaaGTTGAAGCTGAACAGTCGCTGATGACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAACGAGCTCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCGTCGACGt            >  1:1017830/1‑149 (MQ=255)
                                                                                                                                         aGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCCCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTGCCGGCGCGCGGATTATGGGCTTGGAAGTcg          >  1:118991/1‑149 (MQ=25)
                                                                                                                                           ttGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTcgcg        <  1:1087990/149‑1 (MQ=255)
                                                                                                                                           ttGAAGCTGAACAGTCGCTGAGCACCGTAGAAGGCGACAAGGCTGCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGGGAAAGAGCTCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTcgcg        <  1:167122/149‑1 (MQ=255)
                                                                                                                                           ttGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATc                                                                    <  2:1410897/89‑1 (MQ=255)
                                                                                                                                           ttGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATc                                                                    >  1:1410897/1‑89 (MQ=255)
                                                                                                                                            tGAAGCTGAACAGGCGCTGATCACCGTAGAAGGCGACAAGGCGTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCgg       <  1:719983/149‑1 (MQ=255)
                                                                                                                                              aaGCTGAACAGTCGCTGATCACCTTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGt     <  2:338936/149‑1 (MQ=255)
                                                                                                                                               aGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTg    >  2:744828/1‑149 (MQ=255)
                                                                                                                                               aGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTACTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGATATCAAAGTGAACGTTGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTg    <  1:588144/149‑1 (MQ=255)
                                                                                                                                                gCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCGttt                                                                                         >  1:295805/1‑63 (MQ=255)
                                                                                                                                                 cTGAACAGGCGCTGAGCACCCTAGAAGGCGAGAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGACCGTGGGTGACAATGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGaa  <  2:1363996/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CCTGCTCAGGCAGAAGAGAAGAAAGAAGCAGCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAA  >  NC_000913/125947‑126240

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: