Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,550,656 T→A 5.4% Y66N (TAC→AAC) ‡ 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,6560TA 5.4% 285.1 / 3.7 93Y66N (TAC→AAC) ‡fdnHformate dehydrogenase‑N, Fe‑S (beta) subunit, nitrate‑inducible
Reads supporting (aligned to +/- strand):  ref base T (48/40);  new base A (2/3);  total (50/43)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

CCATCACGCCGCCTTCTCAGGTGCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCC  >  NC_000913/1550507‑1550797
                                                                                                                                                      |                                                                                                                                             
ccATCACGCCGCCTTCTCAGGTGCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTTTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGttt                                                                                                                                               <  2:1088559/150‑1 (MQ=255)
 cATCACGCCGCCTTCTCAGGTGCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTa                                                                                                                                              >  1:441401/1‑150 (MQ=255)
 cATAACGCCGCCTTCTCAGGTGCGGGATTAAAAAGCAGAAGTCGCAAAACTTATCGACGTTTACACCTGTATCGGCTGTAAAGCCTGTGAGGTGGCGTGTTAGGAGGGGAACGACCTCCGTGATGAAGTGGGGCACTGCGTC‑GGGGGTTa                                                                                                                                              >  2:988307/1‑150 (MQ=37)
  aTCAGGCCGCCTTCTGAGGTGCGTGATTACAAAGCTGATGTCGCAAAACTTATCGACGTTTCCACCTGTATCGTCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGCC‑GGGGTTTAc                                                                                                                                             <  2:680855/150‑1 (MQ=37)
  aTCACGCCGCCTTCTCAGGTGCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTAc                                                                                                                                             >  2:1107546/1‑150 (MQ=255)
    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                                                                                             >  1:521171/1‑150 (MQ=255)
                                                  cTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCttt                                                                                             <  2:1681139/150‑1 (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:1004289/1‑150 (MQ=255)
                                                                         cGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCaa                                                                      >  2:1016127/1‑150 (MQ=255)
                                                                         cGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCaa                                                                      >  1:1004705/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: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                                                            <  2:649720/150‑1 (MQ=255)
                                                                                   gCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTgg                                                            >  1:668895/1‑150 (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                                   >  2:1284591/1‑150 (MQ=255)
                                                                                                            ggAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGc                                   <  1:1016127/150‑1 (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)
                                                                                                                                                      |                                                                                                                                             
CCATCACGCCGCCTTCTCAGGTGCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCC  >  NC_000913/1550507‑1550797

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