Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,347,285 C→G 5.4% M14I (ATG→ATC)  fumB ← anaerobic class I fumarate hydratase (fumarase B)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,347,2850CG 5.4% 283.6 / 6.7 92M14I (ATG→ATC) fumBanaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base C (43/44);  new base G (3/2);  total (46/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

Only 100 of 101 total aligned reads displayed.

GGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATC  >  NC_000913/4347140‑4347434
                                                                                                                                                 |                                                                                                                                                     
ggAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTACGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCg                                                                                                                                                   <  1:624544/150‑1 (MQ=255)
 gAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTTGATAGTATTCGGTATTGTCTTTCCCCATCgg                                                                                                                                                  >  1:1545290/1‑150 (MQ=255)
 gAGCATAAAAGAAGCGTCGTGAAAGGCATGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAATTCGGCAACGCTAACGTAATCGGAATTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCgg                                                                                                                                                  >  2:563495/1‑150 (MQ=255)
      aaaaaGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAg                                                                                                                                             <  2:184724/149‑1 (MQ=255)
       aaaaGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCGATGGGGAAAga                                                                                                                                            <  1:1316049/150‑2 (MQ=255)
         acgaacCGTCGTGCAAGGCTTGCTGCGCCAGCAGGGGCAGGGCTCCTGGTTCCACTTTCAGGATCGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGAGAGCAGATAGTATTCGGCATTGCCTTTCCCGCTGGGGAAAgaca                                                                                                                                          <  2:1621844/148‑4 (MQ=255)
          agaagCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGc                                                                                                                                         >  1:911745/1‑150 (MQ=255)
             agCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTg                                                                                                                                      >  1:51920/1‑150 (MQ=255)
             agCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTg                                                                                                                                      >  1:1021466/1‑150 (MQ=255)
               cgtcgtGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGt                                                                                                                                    <  2:1358826/150‑1 (MQ=255)
                gtcgtGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCATCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGAATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTa                                                                                                                                   >  1:692226/1‑150 (MQ=255)
                ctcgtGAAAGGCTTGCGGCGCCAGCAGGGGCAGGGCGTCTGGTTCCACTTTCAGGATGGGTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGACTATATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTa                                                                                                                                   <  1:459179/149‑1 (MQ=255)
                 tcgtGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAg                                                                                                                                  >  1:727173/1‑150 (MQ=255)
                  cgtGAAAGGGTTGGTGCGCCTGCAGGGGAAGGGGTTCTGGTTCCTCTTTGAGTATGGTGTCGCGGTCGTACTCGGCAAATCTAACGTATTCGTATGTGAGTTGATTGTATTCGGTTTTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGa                                                                                                                                 <  2:678173/150‑1 (MQ=255)
                  cgtGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGa                                                                                                                                 >  1:39813/1‑150 (MQ=255)
                    tGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTGCTGGGTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATa                                                                                                                               <  1:1288322/150‑1 (MQ=255)
                     gAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATaa                                                                                                                              >  2:717342/1‑150 (MQ=255)
                       aaGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAg                                                                                                                            <  2:500490/150‑1 (MQ=255)
                        aGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAgg                                                                                                                           >  2:1085298/1‑150 (MQ=255)
                         ggCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCAGTCGAAGTCGGCAACGCTAACGTAATCGGATGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAggg                                                                                                                          >  1:715013/1‑150 (MQ=255)
                              gctgcGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGATGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGgtttgt                                                                                                                     >  1:236879/1‑150 (MQ=255)
                                cgcgcCAGCAGGGGCAGGGCGTCTGGTTCCACTTTCATGATGGTTTCGCCGTCGACGTCGTCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGgtttgttt                                                                                                                   <  2:163846/149‑1 (MQ=255)
                                    ccAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACa                                                                                                               <  2:237008/150‑1 (MQ=255)
                                     cagcagGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACAt                                                                                                              <  2:58739/150‑1 (MQ=255)
                                      agcagGGTCAGGTCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAAAGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATa                                                                                                             <  1:522829/150‑1 (MQ=255)
                                          gggTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCtt                                                                                                         <  2:169993/150‑1 (MQ=255)
                                            gTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTTTTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTcc                                                                                                       <  2:197389/150‑1 (MQ=255)
                                              cAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAg                                                                                                     >  2:1444401/1‑150 (MQ=255)
                                              cAGGGCTTCTGGTTCCACTTTCAGGATGGCTTCGCCGTCGAAGTCGGCAACGCTAACTTAATCGTAAGTGAGTAGATATTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAg                                                                                                     <  1:598252/150‑1 (MQ=255)
                                               aGGGATTCTTGTTCCACTTTCAGGATGGTTTCGCTGTCTAAGTCTGCAACGCTAACGTATTGGGAAGTGAGTAGATAGTATTCGGTATTGTATTTCCCCATCGTGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACTTAGCTTCCAGc                                                                                                    <  1:287286/150‑1 (MQ=255)
                                                gggCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGTAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGcc                                                                                                   <  1:627499/150‑1 (MQ=255)
                                                gggCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGcc                                                                                                   >  2:1083007/1‑150 (MQ=255)
                                                gggATTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCGCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGcc                                                                                                   <  2:1724611/150‑1 (MQ=255)
                                                  gCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTg                                                                                                 <  1:1685045/150‑1 (MQ=255)
                                                  gCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTg                                                                                                 >  2:294386/1‑150 (MQ=255)
                                                       tGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAAct                                                                                            >  2:1225532/1‑150 (MQ=255)
                                                        ggTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAActc                                                                                           >  2:1356660/1‑150 (MQ=255)
                                                         gTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAActct                                                                                          >  1:1175642/1‑150 (MQ=255)
                                                          ttCCAGTTTGAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTg                                                                                         <  2:190035/150‑1 (MQ=255)
                                                          ttCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTg                                                                                         <  2:1144671/150‑1 (MQ=255)
                                                             cACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTAt                                                                                      >  1:794069/1‑150 (MQ=255)
                                                              aCTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTAtt                                                                                     <  1:604796/150‑1 (MQ=255)
                                                                tttCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTg                                                                                   <  2:817772/150‑1 (MQ=255)
                                                                tttCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTg                                                                                   >  2:538300/1‑150 (MQ=255)
                                                                 ttCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGt                                                                                  <  2:684589/150‑1 (MQ=255)
                                                                 ttCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGCAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGt                                                                                  <  2:619501/150‑1 (MQ=255)
                                                                 tgaaGGATGGTTTCGCCTTCGAATTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGTTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGt                                                                                  <  2:1653579/147‑1 (MQ=255)
                                                                  tCAGGATGGTTTCGCCGTCGAAGTGGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGtt                                                                                 <  2:1215486/150‑1 (MQ=255)
                                                                   cAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCACCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTa                                                                                >  2:283934/1‑150 (MQ=255)
                                                                      gATGGTTTCGCCGTCTAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTc                                                                             <  1:446452/150‑1 (MQ=255)
                                                                          gTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGaaa                                                                         >  1:786118/1‑150 (MQ=255)
                                                                          gTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGaaa                                                                         >  2:803891/1‑150 (MQ=255)
                                                                             tCGCCGTCGAAGTCGGCAACGCTAACGTAATCGTAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATg                                                                      <  1:126729/150‑1 (MQ=255)
                                                                                ccGTCGAAGTCGGCAACGCTAACGTAATCGTATGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCg                                                                   <  1:504238/150‑1 (MQ=255)
                                                                                 cGTCGAAGTCGGCAACGCTAACGTAATCGTAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGt                                                                  >  1:443320/1‑150 (MQ=255)
                                                                                   gcGAAGTCGGCAACGCTCACGTAATCGGACGTGAGTAGATAGTATTCGGGATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGGTTGTGTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGGGc                                                                <  1:1009829/149‑1 (MQ=255)
                                                                                     gAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCcg                                                              <  2:1021466/150‑1 (MQ=255)
                                                                                        gTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGt                                                           <  1:1291310/150‑1 (MQ=255)
                                                                                           ggCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTACAGCCTGTAACTCTGTATTTGTTATTCGAAAATGTCGTGCCGCGTAGt                                                        >  2:1517025/1‑150 (MQ=255)
                                                                                              aaCGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCAACGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATCGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCAc                                                     >  2:1221601/1‑150 (MQ=255)
                                                                                               aCGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGACTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACg                                                    >  1:1392462/1‑150 (MQ=255)
                                                                                                  cTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGggc                                                 <  1:583133/150‑1 (MQ=255)
                                                                                                   tAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGggcg                                                >  1:919645/1‑150 (MQ=255)
                                                                                                    cttccgaTCTGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGggcgg                                               >  2:1309310/7‑150 (MQ=255)
                                                                                                    aaCGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGggcgg                                               <  1:634660/150‑1 (MQ=255)
                                                                                                         aaTCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGc                                          >  2:1516128/1‑150 (MQ=255)
                                                                                                         aaTCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGT‑ACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCa                                         >  2:937475/1‑150 (MQ=255)
                                                                                                             ggAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGAAAAg                                      >  2:1640732/1‑150 (MQ=255)
                                                                                                                aGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGc                                   <  2:638722/150‑1 (MQ=255)
                                                                                                                aGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGc                                   >  1:1186323/1‑150 (MQ=255)
                                                                                                                aGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGAGTAGTCACGGGCGGCACGCAAAGTGc                                   >  2:1188032/1‑150 (MQ=255)
                                                                                                                aGTGAGTAGATAGTATTCGGTAATGTCATTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTTACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGc                                   <  2:861226/150‑1 (MQ=255)
                                                                                                                   gAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCAtt                                >  2:697900/1‑150 (MQ=255)
                                                                                                                      tAGATAGTATTCGGGATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTATGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTtat                             >  1:149742/1‑150 (MQ=255)
                                                                                                                        gATAGTATTCGGTATTCTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAACGGGTTTGTTTGACATACCATCCAGCCTGTCACGCTGTATTTGTTATTCTATAATGGCGTGCCGCGTAGTCAAGGGCAGCACGCAAAGTGCATTTATaa                           >  2:100598/1‑150 (MQ=255)
                                                                                                                          tAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATTAAGGGTTTGTTTGACATAGCTTCCAGACTGTAACTCTGTATTTGTTATTAGAATATGGCGTGCCGCGTAGTCAAGGGCGGCACGCATAGTGCATTTATAAGa                         >  1:332854/1‑150 (MQ=255)
                                                                                                                            gTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGTGTTTGTTTGACATCGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAAc                       >  2:225523/1‑150 (MQ=255)
                                                                                                                            gTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAAc                       <  2:911745/150‑1 (MQ=255)
                                                                                                                             tcttCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTTTATCTCTGTATTTGTTATTCTAAAATGGCGTGCCGCGTAGTCACGGGCTGCACGCAAAGTGCATTTATAAGAAcc                      <  1:1663512/148‑1 (MQ=255)
                                                                                                                             tATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAAcc                      <  1:717342/150‑1 (MQ=255)
                                                                                                                             tATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGGCACGGGCGGCACGCATAGTGCATTTATAAGAAcc                      <  1:1258858/150‑1 (MQ=255)
                                                                                                                              aTTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAAccc                     >  2:789929/1‑150 (MQ=255)
                                                                                                                              aTTAGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCATGTAACTCTGTATATGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGAGGCACGCAATGTGAATTTATAAGAAccg                     >  1:1240290/1‑149 (MQ=255)
                                                                                                                                   gtattgtattTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACa                <  1:703842/150‑1 (MQ=255)
                                                                                                                                      ttGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATcg             <  1:803891/150‑1 (MQ=255)
                                                                                                                                       tGTCTTTCCCCCTCGGGAAAGGTGCCGGGTAGATAAAGAGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATcgc            >  1:1111841/1‑150 (MQ=255)
                                                                                                                                       caccTTTCCCGATGGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATcgc            >  2:541023/4‑150 (MQ=255)
                                                                                                                                        accTTTCCCGATGGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATcgcg           >  2:1730919/3‑150 (MQ=255)
                                                                                                                                        accTTTCCCGATGGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATcgcg           >  2:717108/3‑150 (MQ=255)
                                                                                                                                          cattCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCAGTAACTCTGACGTTGTTATTCCAACATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCgcggc         <  2:241582/148‑1 (MQ=255)
                                                                                                                                          cTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCgcggc         <  2:570535/150‑1 (MQ=255)
                                                                                                                                              ccccATCTGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAg     <  1:563495/150‑1 (MQ=255)
                                                                                                                                              ccccATCGGGAATGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTAATCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAg     <  2:57772/150‑1 (MQ=255)
                                                                                                                                              ccccATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAg     >  1:1390644/1‑150 (MQ=255)
                                                                                                                                               cccATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGa    >  1:1696802/1‑150 (MQ=255)
                                                                                                                                               cccATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGa    >  2:1282572/1‑150 (MQ=255)
                                                                                                                                               cccATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGa    >  2:93347/1‑150 (MQ=255)
                                                                                                                                                ccATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGAt   <  1:283934/150‑1 (MQ=255)
                                                                                                                                                 cATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATc  >  1:122539/1‑150 (MQ=255)
                                                                                                                                                 cATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATc  <  2:1044175/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GGAGCATAAAAGAAGCGTCGTGAAAGGCTTGCTGCGCCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATC  >  NC_000913/4347140‑4347434

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: ‑