Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,550,660 A→T 5.6% D67V (GAT→GTT) ‡ 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,6600AT 5.6% 257.4 / 4.4 89D67V (GAT→GTT) ‡fdnHformate dehydrogenase‑N, Fe‑S (beta) subunit, nitrate‑inducible
Reads supporting (aligned to +/- strand):  ref base A (43/40);  new base T (2/3);  total (46/43)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

CACGCCGCCTTCTCAGGTGCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTG  >  NC_000913/1550511‑1550809
                                                                                                                                                      |                                                                                                                                                     
cACGCCGCCTTCTCAGGTGCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGa                                                                                                                                                       >  2:627370/1‑150 (MQ=255)
 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                                                                                  >  1:1004705/1‑150 (MQ=255)
                                                                     cGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCaa                                                                                  >  2:1016127/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:1084717/1‑150 (MQ=255)
                                                                       gCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGc                                                                                >  2:759764/1‑150 (MQ=255)
                                                                         tGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTg                                                                              >  2:769383/1‑150 (MQ=255)
                                                                          gTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTgg                                                                             >  2:955891/1‑150 (MQ=255)
                                                                          gTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTgg                                                                             <  2:307639/150‑1 (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:954539/150‑1 (MQ=255)
                                                                                      aGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATcc                                                                 <  1:17548/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                           >  1:1634289/1‑150 (MQ=255)
                                                                                                                            gTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGc                           >  2:185575/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)
                                                                                                                                                      |                                                                                                                                                     
CACGCCGCCTTCTCAGGTGCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTG  >  NC_000913/1550511‑1550809

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