Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 126,392 C→A 6.6% A233D (GCC→GAC)  aceF → pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913126,3920CA 6.6% 377.8 / 8.2 127A233D (GCC→GAC) aceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base C (81/37);  new base A (6/3);  total (87/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 133 total aligned reads displayed.

AGGCGCGGCAGCTCCGGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCG  >  NC_000913/126243‑126541
                                                                                                                                                     |                                                                                                                                                     
aGGCGCGGCAGCTCCGGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGAGAAAGTTGc                                                                                                                                                       >  1:1830685/1‑150 (MQ=255)
 ggCCCGGCAGCTCCGGCCGCTAAACAGGCAGCAGCCCCGGCAGCGGCCCTTGCACCATCTGCTGGCGTGAAAGAAGTTAACGACCCGTATAGCGGCTTTGACTAAGCTTAAGTGACCGAAGTTATGGTGAAAGTGGGCGACAAAGTTGcc                                                                                                                                                      <  2:2087744/150‑1 (MQ=255)
  gcgcGGCAGCTCCGGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCg                                                                                                                                                     <  2:747991/150‑1 (MQ=255)
   cgcgGCAGCTCCGGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAGAGAAGTTAACGGTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGc                                                                                                                                                    >  1:570853/1‑150 (MQ=255)
     cgGCAGCTCCGGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACAAGCTGATGCCGAGAAAGAAGTTCAAGGTCATGATCTCGGCGGTGACGAAGTTGAAGATACTGAAGTGCTGGTGAATGTCGGGGAGACAGTTGCCGcac                                                                                                                                                  >  2:182595/1‑148 (MQ=255)
           cTCCGGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGt                                                                                                                                            >  2:525695/1‑150 (MQ=255)
           cTCCGGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGt                                                                                                                                            >  1:1751503/1‑150 (MQ=255)
           cTCCGGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGt                                                                                                                                            >  2:710140/1‑150 (MQ=255)
               ggCCGCTAAACAGGATGCACATCCGTCATCGGCCTCTCCACCACCGCCGGGCTTGCTATCTGTTAACGTTCCGGATATCGGCTGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACt                                                                                                                                        <  2:1611135/150‑1 (MQ=18)
                gCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTg                                                                                                                                       <  1:1249856/150‑1 (MQ=255)
                gCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTg                                                                                                                                       >  1:660686/1‑150 (MQ=255)
                 ccGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGa                                                                                                                                      >  1:1657237/1‑150 (MQ=255)
                 ccGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGGGATGGGGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTgc                                                                                                                                      >  2:1311742/1‑149 (MQ=255)
                   gCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATc                                                                                                                                    >  1:480343/1‑150 (MQ=255)
                     tAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCAc                                                                                                                                  <  1:126871/150‑1 (MQ=255)
                     tAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCAGGTTATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAATGTTGCCGCTGAACAGTCACTGATCAc                                                                                                                                  <  2:1251830/150‑1 (MQ=255)
                       aaaaGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGCCTGAAGTGATGGTGAATGTGGGCGACAAAGTTGCCGCTGACCAGTCACTGAGAACCg                                                                                                                                >  2:1621716/4‑150 (MQ=255)
                       aaCAGGAAGCAGCTCCGGCAGCGGCCCCAGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATGTCGGCGGTGACGAAGTTGCAGTGACTGAAGTGATTGTGAAAGTGGGCGACAAAGTTGCCGCTGAACCGGCACTGCTCACCg                                                                                                                                >  1:372749/1‑150 (MQ=255)
                        aCAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGGTCCGGATATCGGCGGTGACGAAGTTGAAGTGCCTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGCTCACCGt                                                                                                                               >  1:153001/1‑150 (MQ=255)
                         cAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTa                                                                                                                              >  2:1152512/1‑150 (MQ=255)
                           ggAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTACCGTTCCGGATAACGGCGGTGACGAAGTTGAAGTGACTGAAGAGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGa                                                                                                                            >  1:243937/1‑150 (MQ=255)
                           ggAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAATTTGGCGACAAAGTTTCCGCTGAACAGTCACTGATCACCGTAGa                                                                                                                            >  2:370848/1‑150 (MQ=255)
                           ggAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTGGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGa                                                                                                                            <  1:1886908/150‑1 (MQ=255)
                                 agcTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGa                                                                                                                      >  1:1658089/1‑150 (MQ=255)
                                  gcTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGAATGAAGTGATGGTGAAAGTGGGCTACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGAc                                                                                                                     >  2:1963795/1‑150 (MQ=255)
                                  gcTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGATGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGAc                                                                                                                     >  1:1129230/1‑150 (MQ=255)
                                   cTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGAATGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACa                                                                                                                    >  2:276991/1‑150 (MQ=255)
                                       ggCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGc                                                                                                                >  2:1942152/1‑150 (MQ=255)
                                       ggCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGc                                                                                                                >  1:839282/1‑150 (MQ=255)
                                          aGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTc                                                                                                             >  2:931918/1‑150 (MQ=255)
                                           gCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCt                                                                                                            >  1:1693257/1‑150 (MQ=255)
                                               ccccTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATgg                                                                                                        >  2:1736427/1‑150 (MQ=255)
                                                         aGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGgc                                                                                              >  1:1918763/1‑150 (MQ=255)
                                                                ggCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGttt                                                                                       >  1:1375547/1‑150 (MQ=255)
                                                                  cGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGc                                                                                     >  1:2104656/1‑150 (MQ=255)
                                                                     taaaGACGTTAACGTTCCGGATATCGGCGGCGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAgg                                                                                  <  1:1567912/149‑1 (MQ=255)
                                                                     gAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTTTGGAAGTTCCGGCGCCGTTTGCAgg                                                                                  >  1:195937/1‑150 (MQ=255)
                                                                      aaaGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGc                                                                                 >  2:521773/1‑150 (MQ=255)
                                                                       tagaagTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAGAGATGGGGAAAGTGGGCGACAAAGGTGCCGCTGCACAGTCACTGATCACCGTTGAAGGCGAGAAAGCTTCTATGGAAGTTCCGGCTCCGTTTGCAGGcg                                                                                <  1:300397/149‑1 (MQ=14)
                                                                       aagaagTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGcg                                                                                <  1:601508/150‑1 (MQ=255)
                                                                        agacgTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACAGAAGTCATGGTGAAAGTGGGCGACAAAGTTGAAGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCTGCTCCGTTTGCAGGcga                                                                               >  2:300397/1‑149 (MQ=11)
                                                                        agaagTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGcgt                                                                               >  2:1630741/1‑150 (MQ=255)
                                                                           agTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAATTTCCGGCGCCGTTTGCAGGcgtcgt                                                                            >  2:514059/1‑150 (MQ=255)
                                                                           agTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGcgtcgt                                                                            >  1:889673/1‑150 (MQ=255)
                                                                           agTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGcgtcgt                                                                            >  1:1811187/1‑150 (MQ=255)
                                                                            gTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTg                                                                           >  1:881140/1‑150 (MQ=255)
                                                                                accctCCGGAAACCGGCGGAGACGAAGTTGAAGTGACTTAAGTTATCTTGAAAGTCGTCCACATAGTTGCCGCTGAACATTCACTGATCACCTTATAATGCGACAACGCATCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAgg                                                                       <  2:1175415/146‑1 (MQ=255)
                                                                                 cGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGa                                                                      >  1:2052938/1‑150 (MQ=255)
                                                                                     ccGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTg                                                                  >  1:1991928/1‑150 (MQ=255)
                                                                                       ggATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGGGAAGGAACTGaa                                                                >  1:1120735/1‑150 (MQ=255)
                                                                                         atatCGGCTGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCTCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAg                                                              <  2:1230396/150‑1 (MQ=255)
                                                                                         atatAGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGTGCGACAAATTTTCCGCGGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAg                                                              <  2:348732/150‑1 (MQ=25)
                                                                                           atCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGATTCTATGGAAGTTACGGCGCCGTTTGCAGGCGTCGTGAATGAACTGAAAGta                                                            >  1:365081/1‑149 (MQ=255)
                                                                                            tCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACTAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCa                                                           >  1:1398623/1‑150 (MQ=255)
                                                                                                  gTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTg                                                     >  2:149370/1‑150 (MQ=255)
                                                                                                  gTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTg                                                     >  2:1753409/1‑150 (MQ=255)
                                                                                                   tGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTgg                                                    >  2:1723968/1‑150 (MQ=255)
                                                                                                   tGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTgg                                                    >  1:495135/1‑150 (MQ=255)
                                                                                                   cgACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTgg                                                    >  1:558194/2‑150 (MQ=255)
                                                                                                           ttGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAg                                            >  1:437483/1‑150 (MQ=255)
                                                                                                           ttGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTCCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAg                                            >  1:108182/1‑150 (MQ=255)
                                                                                                              aaGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGa                                         <  2:2051494/150‑1 (MQ=255)
                                                                                                               aGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCAGCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGaa                                        <  2:2085407/150‑1 (MQ=255)
                                                                                                                  gACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAAc                                     >  2:1434589/1‑150 (MQ=255)
                                                                                                                  gACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAAc                                     >  1:230059/1‑150 (MQ=255)
                                                                                                                   aCTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACt                                    <  2:1658089/150‑1 (MQ=255)
                                                                                                                    cTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTg                                   >  2:630955/1‑150 (MQ=255)
                                                                                                                    cTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTg                                   >  1:19148/1‑150 (MQ=255)
                                                                                                                    cTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTg                                   >  1:587346/1‑150 (MQ=255)
                                                                                                                    cTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTg                                   >  1:1741197/1‑150 (MQ=255)
                                                                                                                     tGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTgg                                  <  2:997400/150‑1 (MQ=255)
                                                                                                                      gAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGc                                 >  1:1826489/1‑150 (MQ=255)
                                                                                                                         gTGATTGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAATCTTCTATGGAAGTTCCTGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCg                              >  2:752989/1‑150 (MQ=255)
                                                                                                                         gTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGATCg                              >  1:726612/1‑150 (MQ=255)
                                                                                                                          tGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGc                             >  2:201147/1‑150 (MQ=255)
                                                                                                                          tGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGc                             <  1:931918/150‑1 (MQ=255)
                                                                                                                           gATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCt                            <  2:664561/150‑1 (MQ=255)
                                                                                                                            aTGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTg                           <  2:881140/150‑1 (MQ=255)
                                                                                                                               gTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGAtt                        <  1:131547/150‑1 (MQ=255)
                                                                                                                                tGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTa                       >  1:82592/1‑150 (MQ=255)
                                                                                                                                 gAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTAt                      >  2:2008530/1‑150 (MQ=255)
                                                                                                                                       gggCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCtt                >  2:1873090/1‑150 (MQ=255)
                                                                                                                                        ggCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGATAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCttg               >  2:1410386/1‑149 (MQ=255)
                                                                                                                                        ggCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTc               >  1:671571/1‑150 (MQ=255)
                                                                                                                                         gCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCg              <  2:1294397/150‑1 (MQ=255)
                                                                                                                                         cccaCAAAGTTGCCGCTGAACACTCACTGATCACCCTCCCAGGCCACAAAGCTTCTATGCAAGTTCCGTCGCCGTTTGCAGTCCGAGTGAAGGAACTGAAAGTCAACGTGGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCg              <  2:1782932/147‑1 (MQ=255)
                                                                                                                                          cGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTg                                                                  >  1:147616/1‑97 (MQ=11)
                                                                                                                                          cGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTg                                                                  <  2:147616/97‑1 (MQ=11)
                                                                                                                                          cGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATa                                               >  1:809935/1‑116 (MQ=17)
                                                                                                                                          cGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATa                                               <  2:809935/116‑1 (MQ=17)
                                                                                                                                          cGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGa             >  2:846526/1‑150 (MQ=255)
                                                                                                                                           gACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGaa            >  2:1467803/1‑150 (MQ=255)
                                                                                                                                            aCAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAg           <  1:1189332/150‑1 (MQ=255)
                                                                                                                                             taaaGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGGCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGt          >  2:106171/2‑150 (MQ=18)
                                                                                                                                              aaaGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGtt         >  1:1052689/1‑150 (MQ=255)
                                                                                                                                               aaGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTg        >  2:679485/1‑150 (MQ=255)
                                                                                                                                                   tGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAgg    <  1:759916/150‑1 (MQ=255)
                                                                                                                                                    gCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGgc   >  2:935984/1‑150 (MQ=255)
                                                                                                                                                     ccGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAATTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGgcg  >  1:1642126/1‑150 (MQ=255)
                                                                                                                                                     ccGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGgcg  >  1:728218/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGGCGCGGCAGCTCCGGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCG  >  NC_000913/126243‑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: