Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 126,126 G→A 9.0% K144K (AAG→AAA aceF → pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913126,1260GA 9.0% 371.7 / 12.1 131K144K (AAG→AAAaceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base G (55/64);  new base A (7/5);  total (62/69)
Fisher's exact test for biased strand distribution p-value = 5.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 137 total aligned reads displayed.

GCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAAGC  >  NC_000913/125977‑126275
                                                                                                                                                     |                                                                                                                                                     
gcTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAg                                                                                                                                                       >  2:1762634/1‑150 (MQ=255)
  tCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGc                                                                                                                                                     <  2:94997/150‑1 (MQ=255)
   ccGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCt                                                                                                                                                    >  2:867042/1‑150 (MQ=255)
   ccGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCAGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCt                                                                                                                                                    >  2:1767919/1‑150 (MQ=255)
      gcagcagcaCCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGATGTGACAGAACTCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCt                                                                                                                                                 >  1:921073/1‑150 (MQ=255)
      gcagcagcaCCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACAGAAATCCTGGTGAAAGTTGGCGAGAAAGTTGCCGCTGAACAGTCGCTGATCACCGAAGACGGCGGCGAggaagag                                                                                                                                                 >  1:1817610/1‑145 (MQ=18)
       cagcagcaCCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTa                                                                                                                                                <  2:1713734/150‑1 (MQ=255)
         gcagcaCCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGATCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATg                                                                                                                                              >  2:1930612/1‑150 (MQ=255)
         gcagcaCCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATg                                                                                                                                              >  1:49366/1‑150 (MQ=255)
         gcagcaCCAGCGGCTGCGGAGGCAAAAGACGTTAACGTACCGGATATCGGCAGCGACGCAGGTGAAGTGACCGAAATCCTGGTGAAAGGTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATg                                                                                                                                              >  1:228870/1‑150 (MQ=25)
          cagcaCCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGAGATCGGCAGCGACGAAGATGAAGTGACCGAAATCCTGGTGAAAGTTGGCGAGAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGTAGGCGACAAGGCTTCTATgt                                                                                                                                             >  2:359656/1‑149 (MQ=255)
            gcaCCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGaa                                                                                                                                           >  2:588640/1‑150 (MQ=255)
                 aGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTcc                                                                                                                                      <  2:1838764/150‑1 (MQ=255)
                   cGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAATTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACATTCGCTGATCTCCGTAGATGGCGACAAGGCTTCTATGGAAGTTCCgg                                                                                                                                    <  2:1094881/150‑1 (MQ=255)
                   cGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGGTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCgg                                                                                                                                    >  1:1261510/1‑150 (MQ=255)
                     gCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCt                                                                                                                                  >  1:837742/1‑150 (MQ=255)
                      cTGCGGCGGCAAAAGACGTTAACGTGCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTc                                                                                                                                 <  2:95890/150‑1 (MQ=255)
                               ctcaaGCCGTTAACGTTCCCGACTTGGGCAGCGACGAAGTTGAATTCACCGATATCCTGGTGAAAGTTGGCGACAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTg                                                                                                                        <  1:1858861/147‑1 (MQ=18)
                               cAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTg                                                                                                                        >  2:643320/1‑150 (MQ=255)
                                aaaaGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTgg                                                                                                                       <  1:1950610/150‑1 (MQ=255)
                                            cGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGATGGCACCGTGAAaga                                                                                                           >  2:143842/1‑150 (MQ=255)
                                            cGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGGCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAaga                                                                                                           <  1:1038679/150‑1 (MQ=255)
                                             gTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCGGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAagag                                                                                                          <  2:296958/150‑1 (MQ=255)
                                                 cGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCt                                                                                                                                                 >  1:1450898/1‑107 (MQ=255)
                                                  ggATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCaa                                                                                                     <  1:867043/150‑1 (MQ=255)
                                                    atttCGTCAGCGACGAAGTTGAAGGGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAg                                                                                                   <  1:1272704/147‑1 (MQ=255)
                                                    atatCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAg                                                                                                   <  2:1959778/150‑1 (MQ=255)
                                                    atatCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAg                                                                                                   <  2:1364296/150‑1 (MQ=255)
                                                        cGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGaa                                                                                               >  2:1900269/1‑150 (MQ=255)
                                                            aGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTAt                                                                                                                                               <  1:503748/98‑1 (MQ=255)
                                                            aGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTAt                                                                                                                                               >  2:503748/1‑98 (MQ=255)
                                                              cgacgaAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGAt                                                                                                        <  1:262653/135‑1 (MQ=255)
                                                              cgacgaAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGAt                                                                                                        >  2:262653/1‑135 (MQ=255)
                                                              cgacgaAGTTGAAGTGACCGAAATCCTGGGGAAAGTTGGCGATAAAGCTGAAGCTGCACACTCGCTGCTCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTggg                                                                                         <  2:372972/150‑1 (MQ=21)
                                                                acgaAGTTGAATTGTCCGCAATGCTGTTGAAAGGTGGAGATAAAGTTGCATCTGAACAGTCGCTGATCACCGTAGCTGGCGACAAGGCTTCTATGGAAGGTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTg                                                                                       <  1:339520/150‑1 (MQ=25)
                                                                acgaAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCTCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGtt                                                                                                                             <  2:1264745/112‑1 (MQ=255)
                                                                acgaAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGGGGGTg                                                                                       <  1:840650/150‑1 (MQ=255)
                                                                       tGAAGTGACCGAAATGCTGGTGAAAGATGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCCACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACATAgt                                                                                >  2:38387/1‑150 (MQ=255)
                                                                       tGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAgt                                                                                >  2:1738717/1‑150 (MQ=255)
                                                                       tGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATgg                                                                                                                                             <  2:1029949/89‑1 (MQ=255)
                                                                       tGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATgg                                                                                                                                             >  1:1029950/1‑89 (MQ=255)
                                                                        taaGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGAGCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAgtg                                                                               <  1:110531/149‑1 (MQ=255)
                                                                           gTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCt                                                                            <  2:1261509/150‑1 (MQ=255)
                                                                             gACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGTCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCGGGCCCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTAc                                                                          >  1:331184/1‑150 (MQ=255)
                                                                             gACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGAGAAAGTGTcggc                                                                          >  1:1926364/1‑147 (MQ=255)
                                                                               ccGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCg                                                                        >  1:1405726/1‑150 (MQ=255)
                                                                                   aaTCCTTGTGAAAGTTTTAGATAAAGTTGAAGCTGAACATTCGCGTATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTc                                                                    <  2:249119/150‑1 (MQ=37)
                                                                                    aTCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCg                                                                   <  1:1762635/150‑1 (MQ=255)
                                                                                      ccTGGTGAAAGTTGGCGATAATGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCt                                                                 <  1:1218350/150‑1 (MQ=255)
                                                                                       cTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGGAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTg                                                                >  1:1123448/1‑150 (MQ=255)
                                                                                       cTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTg                                                                >  1:1495070/1‑150 (MQ=255)
                                                                                        tgggtaaaGTTGTCGAAAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGa                                                               <  1:1853532/145‑1 (MQ=255)
                                                                                         ggTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGAt                                                              <  1:398799/150‑1 (MQ=255)
                                                                                          gTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGAtt                                                             >  2:505040/1‑150 (MQ=255)
                                                                                            gAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTAt                                                           <  1:965058/150‑1 (MQ=255)
                                                                                              aaGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCTCCGTGAAAAAAAACAAAGTGAGCTATCCAGTCAAATTGACTACCGGCTAGCTGAGTTTgg                                                         >  1:656389/1‑150 (MQ=17)
                                                                                              aaGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATgg                                                         <  2:1054891/150‑1 (MQ=255)
                                                                                               aGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGt                                                        <  1:588640/150‑1 (MQ=255)
                                                                                               aGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGt                                                        >  2:777192/1‑150 (MQ=255)
                                                                                                   ggCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGATCGCTGATTATGGTCTTc                                                    >  2:1626498/1‑150 (MQ=255)
                                                                                                        tAACGTTGTAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGt                                               <  2:1642000/150‑1 (MQ=255)
                                                                                                          aaGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCgg                                                                                                                                    >  2:1439520/1‑63 (MQ=255)
                                                                                                           aGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTcgc                                            >  1:645213/1‑150 (MQ=255)
                                                                                                           aGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTcgc                                            >  1:284267/1‑150 (MQ=255)
                                                                                                           aGTTGAAGCTGAACAGTCGCTGATCACAGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCGGATTATGGTCTTCGAAGTcgc                                            >  2:1662167/1‑150 (MQ=255)
                                                                                                              tGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCggg                                         <  2:617978/150‑1 (MQ=255)
                                                                                                                aaGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTg                                       >  1:1101315/1‑150 (MQ=255)
                                                                                                                 aGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGa                                      <  2:1193811/150‑1 (MQ=255)
                                                                                                                 aGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGa                                      <  1:1767920/150‑1 (MQ=255)
                                                                                                                  gCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGaa                                     <  1:1633485/150‑1 (MQ=255)
                                                                                                                   cTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGTGTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCTCCGTGAAAGAGATCAAAGTGAACGTGGTTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAg                                    <  1:139978/150‑1 (MQ=255)
                                                                                                                    tGAACAGTCGCTGATCACCGTAGAAGGCGACAAGTCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGc                                   <  1:1274433/150‑1 (MQ=255)
                                                                                                                    tGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGc                                   >  1:744366/1‑150 (MQ=255)
                                                                                                                    tGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCGTCTATGGAAGTTCCGGCTCCGGTTGCTGGCACCGTGACAGAGCTCTACGGGAACGTGGGTGACAAAGTTTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGc                                   <  1:987300/150‑1 (MQ=25)
                                                                                                                     gAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCa                                  >  1:416525/1‑150 (MQ=255)
                                                                                                                     gAACAGTCGCTGAGCACCGTAGAAGGCGAGAAGGATTCTATGTAAGTTCCGGCTGCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGCACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGGCTTCGAAGTCGCGGGTGAAGCa                                  <  1:359656/150‑1 (MQ=255)
                                                                                                                      aaCAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAg                                 <  2:1138604/150‑1 (MQ=255)
                                                                                                                        cAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGgc                               <  1:1900270/150‑1 (MQ=255)
                                                                                                                          gTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGgcgc                             >  1:490241/1‑150 (MQ=255)
                                                                                                                          gTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGgcgc                             <  1:1405370/150‑1 (MQ=255)
                                                                                                                           tCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGcgcg                            >  2:419081/1‑150 (MQ=255)
                                                                                                                           tCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGGTTGCTGGCACCGGGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGcgcg                            <  2:382690/150‑1 (MQ=255)
                                                                                                                             gCTGATCACCTTAGAAGGCGATAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGCGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGGCGCGGGTGAAGCAGGCGCGGc                          <  2:166377/150‑1 (MQ=255)
                                                                                                                              cTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCa                         <  1:1731501/150‑1 (MQ=255)
                                                                                                                              cTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCa                         <  1:1748558/150‑1 (MQ=255)
                                                                                                                              cTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCGTTTGCTGGCCCCGTGAAAGAGATCAAAGTGAACGGGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCa                         <  1:1008654/150‑1 (MQ=255)
                                                                                                                                    aCCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCCACGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGAGAAGCACGAGCGGCAGCTCCg                   >  2:1559159/1‑150 (MQ=255)
                                                                                                                                      cGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAAGGTGGTTGACAAAGCGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGc                 <  1:1548455/150‑1 (MQ=255)
                                                                                                                                      cGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCCGCAGCTCCGGc                 >  2:341461/1‑150 (MQ=255)
                                                                                                                                       gTAGAAGGCGACAAGGCTTCAATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGcc                >  2:805435/1‑150 (MQ=255)
                                                                                                                                        tAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCg               <  2:489736/150‑1 (MQ=255)
                                                                                                                                         aGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAGGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGc              <  2:103923/150‑1 (MQ=255)
                                                                                                                                          gAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGca             >  1:737765/1‑149 (MQ=255)
                                                                                                                                          gAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCt             <  2:49366/150‑1 (MQ=255)
                                                                                                                                          gAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCt             >  2:1908336/1‑150 (MQ=255)
                                                                                                                                               cGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAACCa        <  2:1956336/150‑1 (MQ=255)
                                                                                                                                                 aCAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAgg      <  1:2988/150‑1 (MQ=255)
                                                                                                                                                   aaGGCTTCTATGGAAGTTCCTGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGATCGGCAGCTCCGGCCGCTAAACATGaa    <  2:1697812/150‑1 (MQ=255)
                                                                                                                                                     ggCTTCTATTGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTTCCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAagc  <  1:1195635/150‑1 (MQ=255)
                                                                                                                                                     ggCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAagc  <  1:355321/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAAGC  >  NC_000913/125977‑126275

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: