Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,550,657 A→T 5.5% Y66F (TAC→TTC) ‡ 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,6570AT 5.5% 271.6 / 4.4 91Y66F (TAC→TTC) ‡fdnHformate dehydrogenase‑N, Fe‑S (beta) subunit, nitrate‑inducible
Reads supporting (aligned to +/- strand):  ref base A (46/40);  new base T (2/3);  total (48/43)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

CATCACGCCGCCTTCTCAGGTGCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTG  >  NC_000913/1550508‑1550806
                                                                                                                                                      |                                                                                                                                                     
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                                                                                                      <  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                                                                     <  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: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)
                                                                                                                                                      |                                                                                                                                                     
CATCACGCCGCCTTCTCAGGTGCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTC‑GGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTG  >  NC_000913/1550508‑1550806

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