Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 126,405 A→G 6.6% S237S (TCA→TCG aceF → pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913126,4050AG 6.6% 265.1 / 6.7 93S237S (TCA→TCGaceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base A (53/31);  new base G (3/3);  total (58/35)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCAGCGCCTGCGGC  >  NC_000913/126256‑126554
                                                                                                                                                     |                                                                                                                                                     
ccGGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCa                                                                                                                                                       >  1:1719181/1‑150 (MQ=255)
ccGGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCa                                                                                                                                                       >  2:1100244/1‑150 (MQ=255)
    ccGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATCGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGa                                                                                                                                                   <  1:373858/150‑1 (MQ=255)
     cGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGAGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGAt                                                                                                                                                  >  2:24053/1‑150 (MQ=255)
       cTAAACAGGAAGCAGCTCCTGTAGCGGCCACTGCACAAGCGGCTGGCGTGAAAGAAGTTACAGTTCCGGATATCGGCGGTGACGAAGTAGAAGTTACTGAAGTGATGGTGAAAGTGGGAGACAAAGTTGCCGCTGAACAGTCACTGATCa                                                                                                                                                <  1:770398/150‑1 (MQ=255)
        tAAACAGGAAGCAGCTCCGGCAGCGGCCCCAGCACCAGCAGCAGGCGTGAAAGAAGTTAACGTTCCGGATAACGGCGGTGACGAAGTTGAAGTGACAGAAGAGAAGTTGAAACTGGGCGACAAACATGCCGCTGAACAGTCACTGATCAc                                                                                                                                               >  2:78286/1‑150 (MQ=255)
          aaCAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCg                                                                                                                                             <  1:1362187/150‑1 (MQ=255)
          aaCAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCg                                                                                                                                             >  1:1052687/1‑150 (MQ=255)
          aaCAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGATGCCGCTGAACAGTCACTGATCACCg                                                                                                                                             >  1:1293636/1‑150 (MQ=255)
           aCAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGt                                                                                                                                            >  2:1771326/1‑150 (MQ=255)
            cAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTa                                                                                                                                           <  2:476934/150‑1 (MQ=255)
              ggAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGa                                                                                                                                         >  2:418442/1‑150 (MQ=255)
                aaGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAg                                                                                                                                       <  1:1426998/150‑1 (MQ=255)
                     gcTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCCCTGATCACCGTAGAAGGCGAc                                                                                                                                  >  2:246663/1‑150 (MQ=255)
                          ggCAGCGGCCCCTGCACCAGCGGCTGGCGGGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGc                                                                                                                             <  1:154812/150‑1 (MQ=255)
                             aGCGGCCCCTGCACCCGCGCCTGGCGTGAAAGAAGTTAACGTTCCGGATAACGGCGGTGACTAAGTTGAAGTGACTGAAGTGATGTTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTc                                                                                                                          <  1:797361/150‑1 (MQ=25)
                              gCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCt                                                                                                                         >  1:1160544/1‑150 (MQ=255)
                                ggCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGATGTGACGAAGTTGATGTGTCTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTAt                                                                                                                       >  2:1815086/1‑150 (MQ=255)
                                 gCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAACGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATg                                                                                                                      <  2:1593125/150‑1 (MQ=255)
                                   cccTGCACCAGCGGCTGGCGTGAAAGTAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGa                                                                                                                    <  1:784560/150‑1 (MQ=255)
                                    ccTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGaa                                                                                                                   >  1:184914/1‑150 (MQ=255)
                                     cTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAg                                                                                                                  <  1:618850/150‑1 (MQ=255)
                                        cACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTc                                                                                                               >  2:1007546/1‑150 (MQ=255)
                                         aCCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTcc                                                                                                              >  2:831339/1‑150 (MQ=255)
                                         aCCAGCGGCTGGCGTGAAAGAAGTGAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGGCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTcc                                                                                                              <  1:305371/150‑1 (MQ=255)
                                            aGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGgc                                                                                                           >  1:1302669/1‑150 (MQ=255)
                                            aGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGgc                                                                                                           >  2:1738666/1‑150 (MQ=255)
                                             gCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGgcg                                                                                                          >  2:1361803/1‑150 (MQ=255)
                                             gCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGGGAAAGTGGGCGACAAAGTTGCCGCGGAACAGTCACGGATCACCGTAGAAGGCGACAAAGCTGCTAGGGAAGTTCCGgcg                                                                                                          >  2:1368229/1‑150 (MQ=21)
                                                gCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCg                                                                                                       >  2:284881/1‑150 (MQ=255)
                                                  tGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGtt                                                                                                     >  1:1680363/1‑150 (MQ=255)
                                                  tGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGGCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGtt                                                                                                     <  2:1355832/150‑1 (MQ=255)
                                                     cGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGc                                                                                                  >  2:1619023/1‑150 (MQ=255)
                                                      gTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCa                                                                                                 <  2:1719181/150‑1 (MQ=255)
                                                         aaaGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGc                                                                                              >  1:386135/1‑150 (MQ=255)
                                                          aagaagGTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGcg                                                                                             <  2:526195/150‑1 (MQ=255)
                                                            gaagTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGcgtc                                                                                           >  1:362324/1‑150 (MQ=255)
                                                               gTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTAGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTg                                                                                        >  2:265179/1‑150 (MQ=255)
                                                                   aCGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAgg                                                                                    >  2:1097878/1‑150 (MQ=255)
                                                                    cGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGa                                                                                   >  2:987998/1‑150 (MQ=255)
                                                                    cGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGa                                                                                   >  1:1039934/1‑150 (MQ=255)
                                                                     gTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGACGTGATGGTGAACGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGaa                                                                                  <  2:1856471/150‑1 (MQ=255)
                                                                     gTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGaa                                                                                  >  2:994756/1‑150 (MQ=255)
                                                                     gTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGaa                                                                                  >  1:113342/1‑150 (MQ=255)
                                                                     gTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGaa                                                                                  >  1:1795212/1‑150 (MQ=255)
                                                                     catCCGGATATCGGCGGTGACGAAGTTGAAGTGACTCAAGACATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCCTTTGCAGCCGTCGTGAAGGaa                                                                                  >  1:664951/3‑150 (MQ=21)
                                                                        ccGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTg                                                                               >  1:1576241/1‑150 (MQ=255)
                                                                            atatCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAg                                                                           <  1:9582/150‑1 (MQ=255)
                                                                            atatCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGGGCGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAg                                                                           <  1:66859/150‑1 (MQ=255)
                                                                               tCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAATCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGCAACTGAAAGTCa                                                                        >  1:778264/1‑150 (MQ=255)
                                                                                 ggcggTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAAc                                                                      >  2:705750/1‑150 (MQ=255)
                                                                                   cggTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGt                                                                    <  1:1760388/150‑1 (MQ=255)
                                                                                       gACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGc                                                                <  1:1635474/150‑1 (MQ=255)
                                                                                       gACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGc                                                                >  2:799288/1‑150 (MQ=255)
                                                                                       gACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGc                                                                <  2:1892630/150‑1 (MQ=255)
                                                                                       gACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGc                                                                <  1:833281/150‑1 (MQ=255)
                                                                                          gAAGTTGAAGTGACTGAATTGATGGTGAATGTTGGCGACAAAGTTTCAGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGAt                                                             <  2:818455/150‑1 (MQ=21)
                                                                                            aGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGACCTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATaa                                                           >  1:136232/1‑150 (MQ=255)
                                                                                             gTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATaaa                                                          <  1:1100244/150‑1 (MQ=255)
                                                                                               tGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGt                                                        >  1:1540144/1‑150 (MQ=255)
                                                                                                 aaGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGa                                                      <  2:688248/150‑1 (MQ=255)
                                                                                                  aGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGaa                                                     <  1:265179/150‑1 (MQ=255)
                                                                                                  aGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTAGGCGATAAAGTGaa                                                     >  2:232581/1‑150 (MQ=255)
                                                                                                     gACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCCCTGATCTCCGTAGAAGGCGACAACGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAAc                                                  >  1:106074/1‑150 (MQ=255)
                                                                                                         gAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGc                                              >  2:1276048/1‑150 (MQ=255)
                                                                                                         gAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGc                                              >  1:975503/1‑150 (MQ=255)
                                                                                                          aaGAGATGGTGAAAGTCGTCGTCAAACTTGCCTCTGAACACTCTCTGATCACCGTATATGGCGACAATGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCt                                             <  2:895067/150‑1 (MQ=255)
                                                                                                            gTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCg                                           >  1:344776/1‑150 (MQ=255)
                                                                                                            gTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCg                                           <  2:1807398/150‑1 (MQ=255)
                                                                                                                 ggTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGAt                                      >  1:562835/1‑150 (MQ=255)
                                                                                                                         tGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGAGTATTATCt                              >  2:258690/1‑150 (MQ=255)
                                                                                                                         tGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCCTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGCAAACTGGCTCGCTGATTATGATCt                              >  2:348013/1‑150 (MQ=255)
                                                                                                                           ggCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAATTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTc                            >  2:77357/1‑150 (MQ=255)
                                                                                                                            gCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCg                           >  1:784418/1‑150 (MQ=255)
                                                                                                                             cGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGa                          >  2:1380901/1‑150 (MQ=255)
                                                                                                                               aCAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAg                        <  1:1703252/150‑1 (MQ=255)
                                                                                                                                   aGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGACTATGATCTTCGAAGTTGa                    >  2:1909081/1‑150 (MQ=255)
                                                                                                                                   aGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGc                                                                >  2:144312/1‑106 (MQ=17)
                                                                                                                                   aGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATa                                                            >  1:1266179/1‑110 (MQ=17)
                                                                                                                                   aGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGc                                                                <  1:144312/106‑1 (MQ=255)
                                                                                                                                   aGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATa                                                            <  2:1266179/110‑1 (MQ=255)
                                                                                                                                     ttGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAg                  <  1:1007546/150‑1 (MQ=255)
                                                                                                                                      tGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAgg                 <  1:705750/150‑1 (MQ=255)
                                                                                                                                        ccGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAGGTTCCGGCGCCGTTTGAATGCGTCGTGAAGGAACTGAACGTCAACGGGGGCGATAAAGTGATAACTGGGTCGCTGCTTATGATATTCCAAGTTGAAGgcg               >  2:315599/1‑150 (MQ=25)
                                                                                                                                          gCTGAACAGTCGCTGATCACCGTAGAAGGCGACAATGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACt                                                                                <  2:14744/83‑1 (MQ=255)
                                                                                                                                          gCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACt                                                                                >  1:14744/1‑83 (MQ=255)
                                                                                                                                          gCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCa             >  2:1642388/1‑150 (MQ=255)
                                                                                                                                            tGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCAgc           >  1:420095/1‑150 (MQ=255)
                                                                                                                                             gAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCAgcg          >  2:965360/1‑150 (MQ=255)
                                                                                                                                                  gTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCAGCGCCTGc     >  2:522758/1‑150 (MQ=255)
                                                                                                                                                  gTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCAGCGCCTGc     >  2:1477332/1‑150 (MQ=255)
                                                                                                                                                  gTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCAGCGCCTGc     <  2:1160544/150‑1 (MQ=255)
                                                                                                                                                    cACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCAGCGCCTGCgg   >  1:880796/1‑150 (MQ=255)
                                                                                                                                                    cACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGATCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAACCTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCAGCGCCTGCgg   >  1:432954/1‑150 (MQ=255)
                                                                                                                                                    cACTGATCACCGTAGAAGGCGACAAAGCATCTATGGAAGTTCCGGCGCCGTTTGCAGGCGGCGGGATGGAACTGAAAGTCAACGTGCGCGATAAAGGGAAAACTGGCGCGCTGATTATGATGTTCGAAGTTGAAGGCGCAGCGCCTGCgg   <  1:908793/150‑1 (MQ=25)
                                                                                                                                                     aCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCAGCGCCTGCGGc  <  1:1619023/150‑1 (MQ=255)
                                                                                                                                                     aCTGAGCACCGTAGACGGCGCCAAAGCTTATATGCAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGGGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCAGCGCCTGCGGc  <  1:1165581/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCAGCGCCTGCGGC  >  NC_000913/126256‑126554

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: