Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 126,147 T→G 6.5% A151A (GCT→GCG aceF → pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913126,1470TG 6.5% 198.1 / 2.0 79A151A (GCT→GCGaceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base T (35/37);  new base G (2/4);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTG  >  NC_000913/126001‑126295
                                                                                                                                                  |                                                                                                                                                    
gcggcgGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAATTCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCg                                                                                                                                                   >  1:523074/1‑150 (MQ=255)
 cggcggCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGt                                                                                                                                                  <  2:1300589/150‑1 (MQ=255)
           aGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCAc                                                                                                                                        <  1:520130/150‑1 (MQ=255)
            gACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCAcc                                                                                                                                       >  1:790259/1‑150 (MQ=255)
              cGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAATTCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGACGCTGATCACCGTAGAAGGCGACAAGCCTTCTATGGAAGTTCCGGCTCCGTTTGCAGGctacgg                                                                                                                                     >  2:1015613/1‑145 (MQ=17)
               gTTAACGTTCCGGATATCGGCTGCGACTTAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTg                                                                                                                                    <  2:514458/150‑1 (MQ=255)
               gTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTg                                                                                                                                    >  1:1092180/1‑150 (MQ=255)
               gTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTg                                                                                                                                    <  2:1144694/150‑1 (MQ=255)
                  aaCGTTCCGGATATCGGCAGCTACGAATTTGAAGTTACCGTAAACCTGGTGTAAGTAGGCGTTAAAGCTGAAGCTGAACAGTCGCTGATCACCGTATAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGaaa                                                                                                                                 <  2:616335/150‑1 (MQ=25)
                  aaCGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGaaa                                                                                                                                 >  2:1261815/1‑150 (MQ=255)
                        ccGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATc                                                                                                                           >  2:689385/1‑150 (MQ=255)
                        aagttcatCCGCCGCGAAGAAGTTGAATTCACCGACATCCTGGTCAAAGTTGGCGAAAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCCGCTCCGTTTGCTGGCACCGTGAAAGAGATc                                                                                                                           <  1:531833/144‑1 (MQ=17)
                             tatCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGt                                                                                                                      >  2:855760/1‑150 (MQ=255)
                               tCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGa                                                                                                                    >  2:623431/1‑150 (MQ=255)
                               tCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGa                                                                                                                    >  1:1219230/1‑150 (MQ=255)
                                  gCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCTCCGTTTGCAGGCGACGTGAAAGACCTGAAAGTCAACg                                                                                                                 >  2:633056/1‑150 (MQ=14)
                                    cgCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTTGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTg                                                                                                               <  2:668971/149‑1 (MQ=255)
                                           aaGTTGAAGTGACCGAAATCCTGGTGATAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCTCAGTAGAAGGAGACAAGGCTTCTATGGAAGTTCCGGCTGCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACa                                                                                                        <  1:218764/150‑1 (MQ=25)
                                           aaGTTGAAATCACCGAGATCCTGGTTAAAGTGGGCGACAATGTTGAAGCCGAACAGTCGCTGATCACCGTACAAGGCGACAAGGCATCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACa                                                                                                        <  1:237862/150‑1 (MQ=17)
                                            aGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAg                                                                                                                       <  2:1134053/134‑1 (MQ=11)
                                            aGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAg                                                                                                                       >  1:1134053/1‑134 (MQ=11)
                                                gAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAATGAGCTGAAAGTGAACGTGGGTGAAAAAgtg                                                                                                   >  1:798034/1‑150 (MQ=255)
                                                  aatGACCGAGATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCCGAACAGTAGCTGATCACCGTAGAAGGCGACAATGATTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTc                                                                                                 <  1:882019/148‑1 (MQ=21)
                                                   gagACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTTAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGATCCGTTTGCTGGCGGCGTTATAGATATCAAATTAAACTTGGGAGAAAAAgtgaca                                                                                                >  2:992254/3‑147 (MQ=18)
                                                   gTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCAGGCGACGTGAATGAGCTCAAAGTGAACGTGGGAGAGAAAgtgaca                                                                                                >  2:423705/1‑147 (MQ=18)
                                                     gACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGATCAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTAc                                                                                              >  1:1233041/1‑150 (MQ=255)
                                                      cccGAAATCCTGGTGAAAGTTGGCGAAAAAGTTGAAGCTGAACAGTCGCTGATCACGGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTAcc                                                                                             <  2:457217/149‑1 (MQ=255)
                                                               cTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTg                                                                                    >  1:796141/1‑150 (MQ=255)
                                                                 ggTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGAt                                                                                  >  1:609380/1‑150 (MQ=255)
                                                                 ggTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGAt                                                                                  >  2:837824/1‑150 (MQ=255)
                                                                    gAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTAt                                                                               >  1:1321831/1‑150 (MQ=255)
                                                                    caaaGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAACGCGTCTATGGAAGTTCCGTCTCCGGTTGCTGGAAGCGTGAAAGAGATCAAAGTGAAAGTTGGAGACAAAGAGACTACCGGCTCGCTGATTAt                                                                               >  2:1092904/2‑150 (MQ=17)
                                                                      aaGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAATGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGa                                                                                                                                   <  2:919100/96‑1 (MQ=11)
                                                                      aaGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGa                                                                                                                                   >  1:919100/1‑96 (MQ=11)
                                                                       aGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCATAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAAAGAGATCAAAGTGAACGTTGGAGAAAAAGTGACTACAGGCTCGCTGATTATGGt                                                                            >  1:876008/1‑150 (MQ=18)
                                                                        gatGGCGATAAAGCTGAAGCGGAACTGTCTCTGATCACCGTTGCAGGCGACTAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCATATTGAACGTGGGTGACAAAGTGTCTACCGGAACGCTGATTATGGTc                                                                           <  1:245983/148‑1 (MQ=25)
                                                                          tGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCtt                                                                         >  2:839301/1‑150 (MQ=255)
                                                                             cGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGa                                                                      >  1:1264658/1‑150 (MQ=255)
                                                                               aTAAAGTTGAAGCTGATCAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAg                                                                    <  1:1290355/150‑1 (MQ=255)
                                                                                    gTTGAAGCTGAACAGTCGCTGATCACCTTAGAACGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTcgcg                                                               <  1:779824/150‑1 (MQ=255)
                                                                                     ttGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCgg                                                              <  1:766150/150‑1 (MQ=255)
                                                                                      tGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCggg                                                             <  1:719070/150‑1 (MQ=255)
                                                                                       gAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCCGATTATGGTCTTCGAAGTCGCGGGt                                                            <  1:9826/150‑1 (MQ=255)
                                                                                             gAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAACGTGTCTCCCGGCTCGCTGATTAGGGTCTTCGAAGTCGCGGGTGAAGCa                                                      >  2:200080/1‑150 (MQ=255)
                                                                                             gAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCAACGTGAAAGAGATCAAAGTCAACGTGGGAGAAAAAGTGACAACAGGCTCGCTGATTATGGAGATCGAAGTCGCGGCTgcatcc                                                      >  2:1023264/1‑145 (MQ=17)
                                                                                               aCAGTCTCTGATCACCGTAGAAGGCGACAAGGCGTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAgg                                                    <  2:969468/150‑1 (MQ=255)
                                                                                               aCAGTCGCTGATCACCGTAGACGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGGGAACGTGGGTGACAAAGTGGCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAgg                                                    <  1:490904/150‑1 (MQ=255)
                                                                                                cAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGgc                                                   <  1:1286270/150‑1 (MQ=255)
                                                                                                  gTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGCCGCCGCTGGCTGGCCCCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGgcgc                                                 <  2:239327/150‑1 (MQ=25)
                                                                                                   ccGCTGATCACCGTAGAAGCCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGcgcg                                                <  2:367929/149‑1 (MQ=255)
                                                                                                    cGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCgg                                               >  1:754625/1‑150 (MQ=255)
                                                                                                       tGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAg                                            <  2:1153478/150‑1 (MQ=255)
                                                                                                        gATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGc                                           >  1:685825/1‑150 (MQ=255)
                                                                                                        gATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGc                                           <  1:1028600/150‑1 (MQ=255)
                                                                                                           cACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTcc                                        <  1:623431/150‑1 (MQ=255)
                                                                                                           cACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTcc                                        >  2:746986/1‑150 (MQ=255)
                                                                                                            aCCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGCGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCg                                       <  1:1319671/150‑1 (MQ=255)
                                                                                                             ccGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCgg                                      >  1:259410/1‑150 (MQ=255)
                                                                                                               gaaGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGcc                                    <  1:851598/148‑1 (MQ=255)
                                                                                                               gTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCGGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGcc                                    <  1:1092904/150‑1 (MQ=255)
                                                                                                               gTAGAAGGCGACAAGGCTTCTATGGAAGGTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGcc                                    <  1:809930/150‑1 (MQ=255)
                                                                                                                     ggCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTaaa                              >  2:406923/1‑150 (MQ=255)
                                                                                                                       ctaCAAGGCTTCTATGGAAGTTCCGGCTCCGGTTGCTGGGACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATCGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACa                            <  1:975075/148‑1 (MQ=255)
                                                                                                                       cGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACa                            >  1:943384/1‑150 (MQ=255)
                                                                                                                       cGACAAGGCTTCTAGGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACa                            <  1:497696/150‑1 (MQ=255)
                                                                                                                        gACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAg                           >  1:616282/1‑150 (MQ=255)
                                                                                                                               cTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAagcag                    <  1:855760/150‑1 (MQ=255)
                                                                                                                                ttCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAagcagc                   <  2:1092180/150‑1 (MQ=255)
                                                                                                                                  cTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAACCAGGAAGCAGCTc                 >  1:1104023/1‑150 (MQ=255)
                                                                                                                                     tGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAAGCAGCTCCgg              <  1:20121/150‑1 (MQ=255)
                                                                                                                                     tGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAAGCAGCTCCgg              >  1:775166/1‑150 (MQ=255)
                                                                                                                                      ggAAGTTCCGGCGCAGTTTGCTGGCACCGTGAAAGATATCAAAGTGAACGTGTGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCTGCCGCTAAACAGGAAGCAGCTCCGGc             <  2:1021361/150‑1 (MQ=255)
                                                                                                                                      ggAAGCTCCGGCCCCGTTTGCTGGCACCGAGAAAGAGATCACACTGAACGTGTCTGACAAAGTGTCTACCGTCTCTCTGATTATGGTCTGCGAAGTCGCGGGTGAAGCAGGCGCTGCAGCTCCTGCAGCTAATCAGGAAGCAGCTCCGGc             <  2:910740/150‑1 (MQ=255)
                                                                                                                                       gAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAAGCAGCTCCggac            >  2:226289/1‑148 (MQ=255)
                                                                                                                                        aaGTTCCGGCTCCGGTTGCTCGCACCGTGAAAGAGATCAAAGTCAACGTGGGTGACAACGCGTCTACCCGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAAGCAGCTCCGGCAg           <  2:671379/150‑1 (MQ=255)
                                                                                                                                         aGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAAGCAGCTCCGGCAGc          <  1:689385/150‑1 (MQ=255)
                                                                                                                                          gTTCCTGCTCCGGTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAAGCAGCTCCGGCAGCg         <  2:772281/150‑1 (MQ=255)
                                                                                                                                          gTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAAGCAGCTCCGGCAGCg         >  1:1283535/1‑150 (MQ=255)
                                                                                                                                                 cTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTg  <  2:685825/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTG  >  NC_000913/126001‑126295

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: