Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,550,740 C→A 5.6% L94M (CTG→ATG)  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,7400CA 5.6% 242.8 / 3.6 71L94M (CTG→ATG) fdnHformate dehydrogenase‑N, Fe‑S (beta) subunit, nitrate‑inducible
Reads supporting (aligned to +/- strand):  ref base C (32/35);  new base A (2/2);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.10e-01

CTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGTCGATTTCCAGTCGGAAAACTGCATCGGCTGTGGTTACTGCATTGCCGGGTG  >  NC_000913/1550592‑1550888
                                                                                                                                                    |                                                                                                                                                    
cTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGc                                                                                                                                                      >  2:372902/1‑149 (MQ=255)
  gTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCt                                                                                                                                                     >  7:165873/1‑148 (MQ=255)
   tCAGGTGTCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTg                                                                                                                                                    <  8:89133/148‑1 (MQ=255)
        tGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGt                                                                                                                                              <  8:25960/149‑1 (MQ=255)
          gCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTa                                                                                                                                             >  8:49373/1‑148 (MQ=255)
          gCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTa                                                                                                                                             >  7:120711/1‑148 (MQ=255)
           cGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTaaa                                                                                                                                           <  3:202908/149‑1 (MQ=255)
            gtgtTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTa                                                                                                                                             >  4:42984/1‑146 (MQ=255)
               ttCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGAc                                                                                                                                        <  2:79839/148‑1 (MQ=255)
                  ggAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCt                                                                                                                                    >  7:12502/1‑149 (MQ=255)
                  ggAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCt                                                                                                                                    <  1:350189/149‑1 (MQ=255)
                  ggAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTGTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCt                                                                                                                                    >  4:90723/1‑149 (MQ=255)
                     gTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTa                                                                                                                                 <  1:70103/149‑1 (MQ=255)
                             aCATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACtgt                                                                                                                         <  3:42984/149‑1 (MQ=255)
                                tCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACtgtg                                                                                                                        >  1:161822/1‑147 (MQ=255)
                                    tgatgaAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATc                                                                                                                  <  4:204896/149‑1 (MQ=255)
                                      atgaAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACtg                                                                                                                          >  4:212516/1‑139 (MQ=255)
                                      atgaAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACtg                                                                                                                          <  3:212516/139‑1 (MQ=255)
                                       tgaAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCg                                                                                                               >  8:61149/1‑149 (MQ=255)
                                       tgaAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCg                                                                                                               <  1:395511/149‑1 (MQ=255)
                                           gTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGATGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTg                                                                                                           >  8:102662/1‑149 (MQ=255)
                                            tGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGc                                                                                                          >  7:153751/1‑149 (MQ=255)
                                              gggCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCt                                                                                                        >  3:340950/1‑149 (MQ=255)
                                                 cACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGaa                                                                                                     <  8:120711/149‑1 (MQ=255)
                                                      cGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGt                                                                                                <  5:19350/149‑1 (MQ=255)
                                                                 aCGATAACCCCGCCGATCTGAGCGCCAAGTACTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAATACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTctgct                                                                                     >  4:285796/1‑149 (MQ=255)
                                                                        ccccGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGt                                                                                          >  2:284962/1‑137 (MQ=255)
                                                                        ccccGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGt                                                                                          <  1:284962/137‑1 (MQ=255)
                                                                        ccccGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCGGGTGCaa                                                                              >  6:281290/1‑149 (MQ=255)
                                                                         cccGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGATGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAAt                                                                             >  7:59135/1‑149 (MQ=255)
                                                                         cccGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACAGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAAt                                                                             <  6:317292/149‑1 (MQ=255)
                                                                             ccgATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCAtt                                                                         <  7:376452/149‑1 (MQ=255)
                                                                             ccgATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCAtt                                                                         <  4:256568/149‑1 (MQ=255)
                                                                                 tCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAg                                                                      >  1:199563/1‑148 (MQ=255)
                                                                                     aGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGc                                                                 <  7:49373/149‑1 (MQ=255)
                                                                                     aGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGATGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTc                                                                        <  7:102662/142‑1 (MQ=255)
                                                                                         ccAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCGCTGTGAAGATCCCGGCTGCCt                                                                                                        >  4:9437/1‑106 (MQ=255)
                                                                                         ccAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCGCTGTGAAGATCCCGGCTGCCt                                                                                                        <  3:9437/106‑1 (MQ=255)
                                                                                          cAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACg                                                            >  1:116650/1‑149 (MQ=255)
                                                                                           aaGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACgg                                                           >  5:349381/1‑149 (MQ=255)
                                                                                            aGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACggg                                                          <  5:364263/149‑1 (MQ=255)
                                                                                                cTGGACGGTGATGCGCTTTAGCGAAACCGAACAGATCGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTg                                                      <  5:306777/149‑1 (MQ=255)
                                                                                                   gACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGTc                                                    <  8:165873/148‑1 (MQ=255)
                                                                                                        tGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGTCGATTTc                                              >  1:122315/1‑149 (MQ=255)
                                                                                                             cgcTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGTCGATTTCCAGTc                                         <  7:242600/149‑1 (MQ=255)
                                                                                                                tttAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGATGTATGCACTGTGAAGATCCCGCCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGACGATTTCCAGTCgg                                       >  6:313831/1‑148 (MQ=255)
                                                                                                                   atcGAAACCGAACAGAACGGCAAGCTGCAGTCGCTGATCCGTAAAGACGGCAGTATGAACTGTGAAGAAACCGGCAGCCAGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACCCTAACGCGATTGTCAATTTCAAGTCTGaaa                                    >  4:250245/3‑148 (MQ=255)
                                                                                                                    gCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGTCGATTTCCAGTCGGAAAAc                                  >  4:177993/1‑149 (MQ=255)
                                                                                                                      gAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGTCGATTTCCAGTCGGAAAACTg                                >  8:203065/1‑149 (MQ=255)
                                                                                                                       aaaCCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGTCGATTTCCAGTCGGAAAACTGc                               <  3:90723/149‑1 (MQ=255)
                                                                                                                           cGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGTCGATTTCCAGTCGGAAAACTGCATCg                           >  7:15166/1‑149 (MQ=255)
                                                                                                                            gAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGTCGATTTCCAGTCGGAAAACTGCATCgg                          <  8:394546/149‑1 (MQ=255)
                                                                                                                            gAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGTCGATTTCCAGTCGGAAAACTGCATCgg                          <  3:285796/149‑1 (MQ=255)
                                                                                                                             aaCAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGTCGATTTCCAGTCGGAAAACTGCATCGGc                         <  2:116650/149‑1 (MQ=255)
                                                                                                                                 gAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGTCGATTTCCAGTCGGAAAACTGCATCGGCtgtg                     <  7:225351/149‑1 (MQ=255)
                                                                                                                                 gAACGGCAAGCTGGAGTGGATGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGTCGATTTCCAGTCGGAAAACTGCATCGGCtgtg                     <  8:59135/149‑1 (MQ=255)
                                                                                                                                   aCGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGTCGATTTCCAGTCGGAAAACTGCATCGGCTGTGGt                   <  5:281290/149‑1 (MQ=255)
                                                                                                                                     ggCAAGCTGGAGTGGCTTATCCGTAAAGACTTCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGTCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGTATTGTCGATTTCCAGTCGGAAAACTGCATCGGCTGTGGTTa                 <  3:250245/149‑1 (MQ=255)
                                                                                                                                      gCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGTCGATTTCCAGTCGGAAAACTGCATCGGCTGTGGTTAc                <  6:271303/149‑1 (MQ=255)
                                                                                                                                      gCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGTCGATTTCCAGTCGGAAAACTGCATCGGCTGTGGTTAc                >  6:56944/1‑149 (MQ=255)
                                                                                                                                      gCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGTCGATTTCCAGTCGGAAAACTGCATCGGCTGTGGTTAc                >  4:108618/1‑149 (MQ=255)
                                                                                                                                      gCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGGCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGTCGATTTCCAGTCGGAAAACTGCATCGGCTGTGGTTAc                <  2:161822/149‑1 (MQ=255)
                                                                                                                                         aGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGTCGATTTCCAGTCGGAAAACTGCATCGGCTGTGGTTACTGc             <  1:372902/149‑1 (MQ=255)
                                                                                                                                          gCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGTCGATTTCCAGTCGGAAAACTGCATCGGCTGTGGTTACTGCa            >  4:214619/1‑149 (MQ=255)
                                                                                                                                           cTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGTCGATTTCCAGTCGGAAAACTGc                               >  8:216383/1‑129 (MQ=255)
                                                                                                                                           cTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGTCGATTTCCAGTCGGAAAACTGc                               <  7:216383/129‑1 (MQ=255)
                                                                                                                                           cTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGTCGATTTCCAGTCGGAAAACTGCATCGGCTGTGGTTACTGCAt           <  6:349381/149‑1 (MQ=255)
                                                                                                                                           cTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGTCGATTTCCAGTCGGAAAACTGCATCGGCTGTGGTTACTGCAt           <  4:185756/149‑1 (MQ=255)
                                                                                                                                                gTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGTCGATTTCCAGTCGGAAAACTGCATCGGCTGTGGTTACTGCATTGCCg      >  7:202761/1‑149 (MQ=255)
                                                                                                                                                 tGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCTATCATTCAGTACGCTAACGGGATTGTCGATTTCCAGTCGGAAAACTGCATCGGCTGAGGTTACTGCATTGCCgg     >  1:620/1‑149 (MQ=255)
                                                                                                                                                   gCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTATCGGGATTGTCGATTTCCAGTCGGAAAACTGCATCGGCTGTGGTTACTGCATTGCCGGgt   >  4:151737/1‑149 (MQ=255)
                                                                                                                                                   gCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGTCGATTTCCAGTCGGAAAACTGCATCGGCTGTGGTTACTGCATTGCCGGgt   >  1:310642/1‑149 (MQ=255)
                                                                                                                                                    cTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGTCGATTTCCAGTCGGAAAACTGCATCGGCTGTGGTTACTGCATTGCCGGgtg  <  2:350488/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCTGCTGGTGCAATCATTCAGTACGCTAACGGGATTGTCGATTTCCAGTCGGAAAACTGCATCGGCTGTGGTTACTGCATTGCCGGGTG  >  NC_000913/1550592‑1550888

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: