Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 126,126 G→A 7.4% 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 7.4% 226.0 / 4.3 93K144K (AAG→AAAaceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base G (30/54);  new base A (4/3);  total (35/58)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAAG  >  NC_000913/125978‑126274
                                                                                                                                                    |                                                                                                                                                    
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:17292/1‑149 (MQ=255)
                                                                                 aaaTCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGc                                                                     <  2:1239465/149‑1 (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                                                                                                                                          <  2:404833/63‑1 (MQ=255)
                                                                                                  ggCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGaa                                                                                                                                          >  1:404833/1‑63 (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                                                                                                      >  1:1410897/1‑89 (MQ=255)
                                                                                                            ttGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATc                                                                                                      <  2:1410897/89‑1 (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)
                                                                                                                    tcaCAGTCGCGGCTCACCGTAGATGGCGACAAGGCTTCTATGGATGTTCCGGCTACGTTTGCTGGCACCGGGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGTTTATGTTCTTCGAAGTCGCGGGTGAAGc                                  <  1:853188/147‑1 (MQ=255)
                                                                                                                    gAACAGTCGCTTATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTATGCTGGCACCGTGAAAGAGCTCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGc                                  <  2:1477773/149‑1 (MQ=255)
                                                                                                                    gAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGc                                  >  2:17932/1‑149 (MQ=255)
                                                                                                                      aCAGGCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGGGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAg                                <  2:1443144/149‑1 (MQ=255)
                                                                                                                       cAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAgg                               >  2:1240072/1‑149 (MQ=255)
                                                                                                                       cAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGAGCAAAGTGAACGTGGGTGACAAAGTGGCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAgg                               <  2:625423/149‑1 (MQ=255)
                                                                                                                         gTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGACCGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGGGAAGCAGgcg                             <  1:1319133/149‑1 (MQ=255)
                                                                                                                          tCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGgcgc                            >  2:865528/1‑149 (MQ=255)
                                                                                                                          tCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGgcgc                            >  2:110086/1‑149 (MQ=255)
                                                                                                                            gCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCgg                          >  1:173220/1‑149 (MQ=255)
                                                                                                                             cTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTGCCGGCGCCGTTTGCTGGCACCGGGAAAGAGAGCAAAGTGAACGTGGGGGACAAAGTGTCTACCGGCGCGCTGATTATGGTCTTCGAAGTCGCGGGGGAAGCAGGCGCGGc                         <  1:699569/149‑1 (MQ=25)
                                                                                                                               gATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCac                       >  1:553316/1‑148 (MQ=255)
                                                                                                                               gATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAg                       <  1:822628/149‑1 (MQ=255)
                                                                                                                                   aCCGTAGAAGGCGACAAGCCTTCTATGGAAGTTCCGGCTCCGTTTGCTGCCACCGTGAAAGAGATCAACGTGAACGGGGCTGACAAAGTGTCTACCGGCGCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTcc                   <  2:197108/149‑1 (MQ=255)
                                                                                                                                    ccGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGTATCCGTTTGCTGGCACCGTGAAAGAGATCAAATTGACCGTGTGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCg                  <  1:1503876/149‑1 (MQ=255)
                                                                                                                                    ccGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAACTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCg                  >  1:533181/1‑149 (MQ=255)
                                                                                                                                      gTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTTGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGc                >  2:950481/1‑149 (MQ=255)
                                                                                                                                      ctAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAACGGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGc                <  1:1563102/148‑1 (MQ=255)
                                                                                                                                        ggatggCGACAAGGCTTCAAAGGACGATCCCCCTCCGGATGCGGGCACGGTCAACCACATCACAGGGCCCGGGGGTGACACACTGGCGACCGGCGCGCTGATTCTGGTCTTCGATGTCGCGGGAGACGCAGGCGCGGCAGCTCCGGCCg              <  2:1168604/145‑1 (MQ=255)
                                                                                                                                          aaGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGCAGCAGGCGCGGCAGCTCCGGCCGCt            >  2:1521262/1‑149 (MQ=255)
                                                                                                                                              cGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAAc        <  2:619634/149‑1 (MQ=255)
                                                                                                                                              cGACAAGGCTTCTATGGAAGCTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAAc        <  2:386524/149‑1 (MQ=255)
                                                                                                                                               gACAAGGCTTCTATTGAAGTTCCGCCGCCCTTTGCTGGCACCGTCAAAGACATCAAAGTGAACGGGGGTGACAACGTGCCTACCGGCTCTCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACa       <  1:1195264/149‑1 (MQ=255)
                                                                                                                                               gACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACAGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACa       >  2:777215/1‑149 (MQ=255)
                                                                                                                                                aCAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAg      <  2:1126016/149‑1 (MQ=255)
                                                                                                                                                  aaGGCATCTATGGAAGTTCCGGCTCCGTTTGATGGAACAGTGAAAGACCTCAAAGTCAACGTGGGTGACAAAGTTTCTAAAGTCTCGATGATTATGTTCTTATAAGTCGCGGTTGAAGCATTCTCGGCAGCTCCGGCCGCTAAACAGGa    <  2:1074230/149‑1 (MQ=255)
                                                                                                                                                   tggCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGTGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGaa   <  1:1445402/148‑1 (MQ=255)
                                                                                                                                                    ggCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAag  >  1:1302849/1‑149 (MQ=255)
                                                                                                                                                    ggCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAag  <  1:724101/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTCCGGCAGCAGCACCAGCGGCTGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAAG  >  NC_000913/125978‑126274

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: