Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 126,405 A→G 8.8% 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 8.8% 391.6 / 20.9 139S237S (TCA→TCGaceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base A (83/42);  new base G (6/6);  total (89/50)
Fisher's exact test for biased strand distribution p-value = 3.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 141 total aligned reads displayed.

GGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCAGCGCCTGCGGC  >  NC_000913/126258‑126554
                                                                                                                                                   |                                                                                                                                                     
ggCCGCTAAACAGGATGCACATCCGTCATCGGCCTCTCCACCACCGCCGGGCTTGCTATCTGTTAACGTTCCGGATATCGGCTGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACt                                                                                                                                                     <  2:1611135/150‑1 (MQ=18)
 gCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTg                                                                                                                                                    >  1:660686/1‑150 (MQ=255)
 gCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTg                                                                                                                                                    <  1:1249856/150‑1 (MQ=255)
  ccGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGa                                                                                                                                                   >  1:1657237/1‑150 (MQ=255)
  ccGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGGGATGGGGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTgc                                                                                                                                                   >  2:1311742/1‑149 (MQ=255)
  ccGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGAAATCGGCGGTGCCGAAGGTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGa                                                                                                                                                   >  1:59902/1‑150 (MQ=255)
    gCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATc                                                                                                                                                 >  1:480343/1‑150 (MQ=255)
      tAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCAc                                                                                                                                               <  2:908192/150‑1 (MQ=255)
      tAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCAc                                                                                                                                               <  1:126871/150‑1 (MQ=255)
      tAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCAGGTTATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAATGTTGCCGCTGAACAGTCACTGATCAc                                                                                                                                               <  2:1251830/150‑1 (MQ=255)
        aaaaGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGCCTGAAGTGATGGTGAATGTGGGCGACAAAGTTGCCGCTGACCAGTCACTGAGAACCg                                                                                                                                             >  2:1621716/4‑150 (MQ=255)
        aaCAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCg                                                                                                                                             <  1:435407/150‑1 (MQ=255)
        aaCAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCg                                                                                                                                             >  1:1663762/1‑150 (MQ=255)
        aaCAGGAAGCAGCTCCGGCAGCGGCCCCAGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATGTCGGCGGTGACGAAGTTGCAGTGACTGAAGTGATTGTGAAAGTGGGCGACAAAGTTGCCGCTGAACCGGCACTGCTCACCg                                                                                                                                             >  1:372749/1‑150 (MQ=255)
         aCAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGGTCCGGATATCGGCGGTGACGAAGTTGAAGTGCCTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGCTCACCGt                                                                                                                                            >  1:153001/1‑150 (MQ=255)
          cAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTa                                                                                                                                           >  2:1152512/1‑150 (MQ=255)
          aagcaagaagCACCGGCAGCCGCCCCAGCACCAGCTGCTGGCGCGAAAGAAGTTAACGTTCAGTATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTTATTGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTa                                                                                                                                           <  2:1645134/146‑1 (MQ=255)
            ggAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAATTTGGCGACAAAGTTTCCGCTGAACAGTCACTGATCACCGTAGa                                                                                                                                         >  2:370848/1‑150 (MQ=255)
            ggAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTGGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGa                                                                                                                                         <  1:1886908/150‑1 (MQ=255)
              aaGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAg                                                                                                                                       >  1:29996/1‑150 (MQ=255)
                  agcTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGCCTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGa                                                                                                                                   >  2:1870042/1‑150 (MQ=255)
                   gcTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGAATGAAGTGATGGTGAAAGTGGGCTACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGAc                                                                                                                                  >  2:1963795/1‑150 (MQ=255)
                   gcTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGATGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGAc                                                                                                                                  >  1:1129230/1‑150 (MQ=255)
                    cTCCGGCAGCGGCCCCTGCACCATCGGCTTGCGTGAAAGAACTTTACGTTCCGGAATTCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACa                                                                                                                                 <  1:1810257/150‑1 (MQ=25)
                      ccGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACaaa                                                                                                                               >  1:1423353/1‑150 (MQ=255)
                        ggCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGc                                                                                                                             >  2:1942152/1‑150 (MQ=255)
                        ggCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGc                                                                                                                             >  1:839282/1‑150 (MQ=255)
                         gCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCt                                                                                                                            >  2:591349/1‑150 (MQ=255)
                            gCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCt                                                                                                                         >  1:1693257/1‑150 (MQ=255)
                                ccccTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATgg                                                                                                                     >  2:1736427/1‑150 (MQ=255)
                                 cccTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGa                                                                                                                    >  1:1870500/1‑150 (MQ=255)
                                       aCCAGCGGCTGGGGTTAAAGAAGTTAACGTTCCGGATATCGGCGGGGACGAAGTTGAAGTGACTGAATTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTcc                                                                                                              <  2:1553227/150‑1 (MQ=255)
                                          aGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGgc                                                                                                           >  1:1918763/1‑150 (MQ=255)
                                                   cGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGc                                                                                                  >  1:2104656/1‑150 (MQ=255)
                                                      taaaGACGTTAACGTTCCGGATATCGGCGGCGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAgg                                                                                               <  1:1567912/149‑1 (MQ=255)
                                                      gAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAgg                                                                                               >  2:1022589/1‑150 (MQ=255)
                                                      gAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAgg                                                                                               >  1:495478/1‑150 (MQ=255)
                                                        tagaagTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAGAGATGGGGAAAGTGGGCGACAAAGGTGCCGCTGCACAGTCACTGATCACCGTTGAAGGCGAGAAAGCTTCTATGGAAGTTCCGGCTCCGTTTGCAGGcg                                                                                             <  1:300397/149‑1 (MQ=14)
                                                         agacgTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACAGAAGTCATGGTGAAAGTGGGCGACAAAGTTGAAGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCTGCTCCGTTTGCAGGcga                                                                                            >  2:300397/1‑149 (MQ=11)
                                                            agTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAATTTCCGGCGCCGTTTGCAGGcgtcgt                                                                                         >  2:514059/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)
                                                                          atatCGGCTGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCTCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAg                                                                           <  2:1230396/150‑1 (MQ=255)
                                                                          atatCGGCTGCGACGAAGTTGAAGTGACTGAAATGATGGTGAAAGTGGGCGAAAAAGTTGAAGCTGAACAGTCACTGATCACCGTAGAAGGCGACAACGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAg                                                                           <  2:1304094/150‑1 (MQ=17)
                                                                          atatAGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGTGCGACAAATTTTCCGCGGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAg                                                                           <  2:348732/150‑1 (MQ=25)
                                                                             tCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACTAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCa                                                                        >  1:1398623/1‑150 (MQ=255)
                                                                                   gTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTg                                                                  >  2:1753409/1‑150 (MQ=255)
                                                                                    tGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTgg                                                                 >  1:495135/1‑150 (MQ=255)
                                                                                    tGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTgg                                                                 >  2:781534/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                                                  >  1:230059/1‑150 (MQ=255)
                                                                                                   gACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAAc                                                  >  2:1434589/1‑150 (MQ=255)
                                                                                                    aCTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACt                                                 <  2:1658089/150‑1 (MQ=255)
                                                                                                     cTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTg                                                >  2:630955/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)
                                                                                                      tGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTgg                                               <  1:1152512/150‑1 (MQ=255)
                                                                                                       gAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGc                                              >  1:2067120/1‑150 (MQ=255)
                                                                                                       gAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGc                                              >  1:1826489/1‑150 (MQ=255)
                                                                                                         aGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTc                                            <  2:1751503/150‑1 (MQ=255)
                                                                                                          gTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGATCg                                           >  1:726612/1‑150 (MQ=255)
                                                                                                           tGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGc                                          >  2:201147/1‑150 (MQ=255)
                                                                                                           tGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGc                                          <  1:931918/150‑1 (MQ=255)
                                                                                                             aTGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTg                                        <  2:29996/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)
                                                                                                                 tGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTa                                    <  1:1942152/150‑1 (MQ=255)
                                                                                                                      gTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATc                               <  2:1965963/150‑1 (MQ=255)
                                                                                                                         ggCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGATAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCttg                            >  2:1410386/1‑149 (MQ=255)
                                                                                                                         ggCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTc                            >  1:671571/1‑150 (MQ=255)
                                                                                                                           cGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTg                                                                               >  1:147616/1‑97 (MQ=11)
                                                                                                                           cGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTg                                                                               <  2:147616/97‑1 (MQ=11)
                                                                                                                           cGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATa                                                            <  2:809935/116‑1 (MQ=17)
                                                                                                                            gACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGaa                         >  1:1616282/1‑150 (MQ=255)
                                                                                                                            gACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGaa                         >  2:1467803/1‑150 (MQ=255)
                                                                                                                             aCAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAg                        <  1:1189332/150‑1 (MQ=255)
                                                                                                                               aaaGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGtt                      >  1:1052689/1‑150 (MQ=255)
                                                                                                                                aaGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTg                     >  2:679485/1‑150 (MQ=255)
                                                                                                                                    tGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAgg                 <  1:759916/150‑1 (MQ=255)
                                                                                                                                     gaagCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAg                                                                           <  1:349235/88‑1 (MQ=255)
                                                                                                                                     gCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGgc                >  2:935984/1‑150 (MQ=255)
                                                                                                                                      ccGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAATTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGgcg               >  1:1642126/1‑150 (MQ=255)
                                                                                                                                      ccGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGgcg               >  1:728218/1‑150 (MQ=255)
                                                                                                                                        gCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCa             >  2:273261/1‑150 (MQ=255)
                                                                                                                                        gCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCa             <  1:686498/150‑1 (MQ=255)
                                                                                                                                          tGAACAGTCGCTGATCACCGTAGAAGGCGACAATGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCg                                                               <  2:72732/98‑1 (MQ=255)
                                                                                                                                           gAACAGTCTCTGATCACCGTAGAAGGCGACAATGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCAgcg          <  2:1826087/150‑1 (MQ=255)
                                                                                                                                           gAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCAgcg          <  2:660686/150‑1 (MQ=255)
                                                                                                                                             aCAGTCGCTGATCACCGTAGAAGGCGACAATGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAgg                                                                                    <  1:1106955/74‑1 (MQ=255)
                                                                                                                                             aCAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAgg                                                                                    >  2:1106955/1‑74 (MQ=255)
                                                                                                                                              cAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTTTGATCTTCGAAGTTGAAGGCGCAGCGCCt       >  2:78871/1‑150 (MQ=255)
                                                                                                                                                gTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCAGCGCCTGc     <  1:1311742/150‑1 (MQ=255)
                                                                                                                                                 tCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCAGCGCCTGCg    <  1:710140/150‑1 (MQ=255)
                                                                                                                                                   aCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCAGCGCCTGCGGc  <  2:1423353/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTGCACCAGCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCAGCGCCTGCGGC  >  NC_000913/126258‑126554

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: