Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 126,393 C→A 5.0% A233A (GCC→GCA aceF → pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913126,3930CA 5.0% 279.5 / 2.7 92A233A (GCC→GCAaceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base C (52/35);  new base A (2/3);  total (54/38)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGGCAGCTCCGGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCG  >  NC_000913/126246‑126541
                                                                                                                                                   |                                                                                                                                                    
ctcgtcAGCTCAGGCAGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACAAGCGGGTGGCGTGAAAGAAGTTAACGTTGCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGc                                                                                                                                                    <  1:351188/145‑1 (MQ=255)
 ccgGCAGCTCCGGCCGCTAAACAGGACGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCt                                                                                                                                                   <  2:1463043/149‑1 (MQ=255)
   ggCAGCTCCGGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGa                                                                                                                                                 <  1:1278081/150‑1 (MQ=255)
     cAGCTCCGGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAAc                                                                                                                                               <  2:908378/150‑1 (MQ=255)
     cAGCTCCGGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGGTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGCCAAAGTTGCCGCTGAAc                                                                                                                                               >  2:1541186/1‑150 (MQ=255)
      accTCCGGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGAGATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACa                                                                                                                                              >  1:952176/3‑150 (MQ=255)
         tCCGGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTc                                                                                                                                           <  1:204746/150‑1 (MQ=255)
          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:1052687/1‑150 (MQ=255)
                    aaCAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCg                                                                                                                                <  1:1362187/150‑1 (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                                                                                              >  2:1738666/1‑150 (MQ=255)
                                                      aGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGgc                                                                                              >  1:1302669/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)
                                                                                  ccGGATATCGGCAGCGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTg                                                                                                                                                  <  2:300727/70‑1 (MQ=255)
                                                                                  ccGGATATCGGCAGCGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTg                                                                                                                                                  >  1:300727/1‑70 (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:1892630/150‑1 (MQ=255)
                                                                                                 gACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGc                                                   >  2:799288/1‑150 (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                                 >  1:975503/1‑150 (MQ=255)
                                                                                                                   gAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGc                                 >  2:1276048/1‑150 (MQ=255)
                                                                                                                    aaGAGATGGTGAAAGTCGTCGTCAAACTTGCCTCTGAACACTCTCTGATCACCGTATATGGCGACAATGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCt                                <  2:895067/150‑1 (MQ=255)
                                                                                                                      gTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCg                              <  2:1807398/150‑1 (MQ=255)
                                                                                                                      gTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCg                              >  1:344776/1‑150 (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)
                                                                                                                                                   |                                                                                                                                                    
CGCGGCAGCTCCGGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCG  >  NC_000913/126246‑126541

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: