Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,550,661 T→A 6.5% D67E (GAT→GAA) ‡ fdnH → formate dehydrogenase‑N, Fe‑S (beta) subunit, nitrate‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,550,6610TA 6.5% 259.7 / 4.3 89D67E (GAT→GAA) ‡fdnHformate dehydrogenase‑N, Fe‑S (beta) subunit, nitrate‑inducible
Reads supporting (aligned to +/- strand):  ref base T (44/39);  new base A (2/4);  total (46/43)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCCGCCTTCTCAGGTGCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGC  >  NC_000913/1550512‑1550810
                                                                                                                                                      |                                                                                                                                                     
aCGCCGCCTTCTCAGGTGCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTGTACGAt                                                                                                                                                       <  2:1621379/150‑1 (MQ=255)
    cgccTTCTCAGGTGCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAAcc                                                                                                                                                   >  1:1688763/1‑150 (MQ=255)
         tctcAGGTGCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCccgcc                                                                                                                                              >  2:690952/1‑150 (MQ=255)
          ctcAGGTGCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCccgccg                                                                                                                                             <  2:473189/150‑1 (MQ=255)
            cAGGTGCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGAt                                                                                                                                           >  1:649336/1‑150 (MQ=255)
               gTGCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGGTTACGATAACCCCGCCGATCTg                                                                                                                                        >  1:766889/1‑150 (MQ=255)
                 gCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAg                                                                                                                                      >  1:253851/1‑150 (MQ=255)
                 gCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGGCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGGC‑GGGGTTATCGTAAACCccgacgcagtgcc                                                                                                                                      <  2:520559/150‑8 (MQ=25)
                  cGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAgc                                                                                                                                     >  1:1610717/1‑150 (MQ=255)
                     gATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGcc                                                                                                                                  >  2:862078/1‑150 (MQ=255)
                      aTTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCa                                                                                                                                 >  1:1719440/1‑150 (MQ=255)
                          cAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTc                                                                                                                             >  2:954539/1‑150 (MQ=255)
                          cAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTATCGTAAACCccgacg                                                                                                                                             >  1:174712/1‑134 (MQ=25)
                          cAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGGGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTATCGTAAACCccgacg                                                                                                                                             <  2:174712/134‑1 (MQ=25)
                           aaaGCAGAAGTCGCAAAACTTATCGCCGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTCGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACTCTAACCCCGCCGATCTGAGCGCCAAGTcc                                                                                                                            >  1:1140815/1‑150 (MQ=255)
                              gCAGAAGTCGCAACACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTgg                                                                                                                         <  2:1226077/150‑1 (MQ=255)
                               cAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGCGCGCCAAGTCCTGGa                                                                                                                        >  1:1273017/1‑150 (MQ=255)
                               cAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGa                                                                                                                        >  1:514387/1‑150 (MQ=255)
                                aGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGAc                                                                                                                       >  2:617953/1‑150 (MQ=255)
                                   aGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGt                                                                                                                    >  2:467848/1‑150 (MQ=255)
                                     tCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGa                                                                                                                  <  2:799740/150‑1 (MQ=255)
                                      cGCAAAACTTATCGACGTTTCCACCTTTATCGGCTGTTAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGACGATCTGAGCGCCAAGTCCTGGACGGTGAt                                                                                                                 <  1:852051/150‑1 (MQ=255)
                                            aCTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCtt                                                                                                           >  2:1226036/1‑150 (MQ=255)
                                             cTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCttt                                                                                                          <  2:1681139/150‑1 (MQ=255)
                                             cTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCttt                                                                                                          >  1:521171/1‑150 (MQ=255)
                                                aTCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGc                                                                                                       <  2:1148399/150‑1 (MQ=255)
                                                 tCGACGTTTCCACCTGTATCGGCGGTAAAGCCTGGCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCg                                                                                                      <  2:992362/150‑1 (MQ=255)
                                                  cGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCTCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGa                                                                                                     >  2:136828/1‑150 (MQ=255)
                                                   gACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCATGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGaa                                                                                                    <  2:941172/150‑1 (MQ=255)
                                                        ttCCACCTGTATCGGCTGTACAGCCGGGCAGGTGGCGTGTTCGGCGTGGCACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGa                                                                                               <  1:869758/150‑1 (MQ=255)
                                                        ttCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGa                                                                                               <  2:1731305/150‑1 (MQ=255)
                                                                 tATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGGGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACgg                                                                                      <  1:1375666/150‑1 (MQ=255)
                                                                  aTCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATACCCCCTCCGATCTGAGCGCCAAGTCCTGGACGGTGATGAGCTTTAGCGAAACCGAACAGAACTGc                                                                                     >  2:162626/1‑150 (MQ=255)
                                                                  aTCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGc                                                                                     <  2:441401/150‑1 (MQ=255)
                                                                   tCGGCTGTAACGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCa                                                                                    <  1:969351/150‑1 (MQ=255)
                                                                    cGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGCACAGAACGGCaa                                                                                   >  2:1068844/1‑150 (MQ=255)
                                                                    cGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCaa                                                                                   >  2:1016127/1‑150 (MQ=255)
                                                                    cGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCaa                                                                                   >  1:1004705/1‑150 (MQ=255)
                                                                    cGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCaa                                                                                   >  1:1004289/1‑150 (MQ=255)
                                                                     ggCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATTTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAg                                                                                  >  2:1273345/1‑150 (MQ=255)
                                                                     ggCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAg                                                                                  >  1:1262567/1‑150 (MQ=255)
                                                                      gCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGc                                                                                 >  2:759764/1‑150 (MQ=255)
                                                                      gCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGc                                                                                 >  2:1084717/1‑150 (MQ=255)
                                                                        tGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTg                                                                               >  2:769383/1‑150 (MQ=255)
                                                                         gTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTgg                                                                              <  2:307639/150‑1 (MQ=255)
                                                                         gTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTgg                                                                              >  2:955891/1‑150 (MQ=255)
                                                                            aaGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCCCTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGCTGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGt                                                                           >  2:301907/1‑150 (MQ=255)
                                                                             aGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTg                                                                          >  2:357034/1‑150 (MQ=255)
                                                                             aGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTg                                                                          >  2:1365838/1‑150 (MQ=255)
                                                                              gCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTgg                                                                         >  1:668895/1‑150 (MQ=255)
                                                                              gCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTgg                                                                         <  2:649720/150‑1 (MQ=255)
                                                                               ccTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGc                                                                        <  1:1245583/150‑1 (MQ=255)
                                                                                cTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGTTTAAGTGGTGCACTGCGTC‑GTGGTTTACGATAACCGCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGAGGCt                                                                       <  1:1202635/150‑1 (MQ=255)
                                                                                cTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCt                                                                       <  2:1505469/150‑1 (MQ=255)
                                                                                cTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCt                                                                       <  2:1691913/150‑1 (MQ=255)
                                                                                cTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGCGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCt                                                                       <  1:1566076/150‑1 (MQ=255)
                                                                                   gcAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGAt                                                                    <  1:1120678/149‑1 (MQ=255)
                                                                                    cAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATc                                                                   <  1:849245/150‑1 (MQ=255)
                                                                                     aGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATcc                                                                  <  1:17548/150‑1 (MQ=255)
                                                                                     aGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATcc                                                                  <  1:954539/150‑1 (MQ=255)
                                                                                       gTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGt                                                                >  1:583697/1‑150 (MQ=255)
                                                                                           cGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGACTGATCCGTAAAg                                                            >  2:367830/1‑150 (MQ=255)
                                                                                             tgtTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGAc                                                          >  1:1447033/1‑150 (MQ=255)
                                                                                                tCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGc                                                       <  1:997630/150‑1 (MQ=255)
                                                                                                tCGGAGGGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGc                                                       <  1:93462/150‑1 (MQ=255)
                                                                                                 gggAGTGGAACGACATCCGGGATGAAGTGGGGCACTGCGGC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCATAAAGACGGCt                                                      <  2:87203/149‑1 (MQ=255)
                                                                                                  ggAGTGGAACGACATCCGTGATGACGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTg                                                     <  1:101250/150‑1 (MQ=255)
                                                                                                   gAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGt                                                    <  2:514387/150‑1 (MQ=255)
                                                                                                      cggAACGACATCCGTCCTGACCTGGCGCACTGCGGC‑GGGCTTTACGAAAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATg                                                 <  2:1176414/149‑1 (MQ=255)
                                                                                                       ggAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGc                                                <  1:1016127/150‑1 (MQ=255)
                                                                                                       ggAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGc                                                >  2:1284591/1‑150 (MQ=255)
                                                                                                        gAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCa                                               >  1:953230/1‑150 (MQ=255)
                                                                                                         aaCGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCAc                                              <  2:253851/150‑1 (MQ=255)
                                                                                                            gACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATATCCCCGCCGAGCTGAGCGCCAAGTCCTGGACGGTGATGGGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACtgt                                           <  2:1610717/150‑1 (MQ=255)
                                                                                                              cATCCGTGATGAAGGGGGGCACTGCGTCGGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACtgtg                                          <  1:584023/150‑1 (MQ=255)
                                                                                                                tCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAg                                       >  2:1725594/1‑150 (MQ=255)
                                                                                                                tCCGTGATGAAGTGGGGCACGGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAg                                       <  1:1471844/150‑1 (MQ=255)
                                                                                                                  cGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGAt                                     <  2:1688763/150‑1 (MQ=255)
                                                                                                                       tgaAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCgg                                <  1:1697558/150‑1 (MQ=255)
                                                                                                                          aGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTg                             >  1:1697250/1‑150 (MQ=255)
                                                                                                                           gTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGc                            >  2:185575/1‑150 (MQ=255)
                                                                                                                           gTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGc                            >  1:1634289/1‑150 (MQ=255)
                                                                                                                             ggggCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCt                          >  2:1620743/1‑150 (MQ=255)
                                                                                                                                   cTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGc                    >  1:722390/1‑150 (MQ=255)
                                                                                                                                       gTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGc                >  1:1580043/1‑150 (MQ=255)
                                                                                                                                        tC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGcc               <  1:467848/150‑1 (MQ=255)
                                                                                                                                        gc‑GGGGTTATCGTAAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAgg                     <  1:388852/143‑1 (MQ=255)
                                                                                                                                        gc‑GGGGTTATCGTAAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAgg                     >  2:388852/2‑144 (MQ=255)
                                                                                                                                                  aCGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTctgctg      <  1:185575/150‑1 (MQ=255)
                                                                                                                                                     aTAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTg   <  2:668895/150‑1 (MQ=255)
                                                                                                                                                      tAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGc  <  2:264639/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
ACGCCGCCTTCTCAGGTGCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGC  >  NC_000913/1550512‑1550810

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