Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,287,037 C→T 13.0% intergenic (+511/+29) narI → / ← rttR nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,287,0370CT13.0% 132.7 / 9.8 62intergenic (+511/+29)narI/rttRnitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base C (25/28);  new base T (2/6);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAAC  >  NC_000913/1286785‑1287322
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                
tgcgtcggaagcctaagctgggTCGAAGAGTCTGCGCCTTCACCGTATGGGCAGAGTATTTGATTGGGGAGTCGGTTCAGGTTTCCGTTGCTTATGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTGGCCTGCGTCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggtg                                                                                                                                                                                                                                                              <  1:855270/269‑1 (MQ=14)
    gttcgaggctgaCCTTGGGTGAAGCTGGACAGCCTTCGCAGTGTGTGCAGAATGTAGGATTAGGGATTCGATGGAGAATTGCGGCGCTTGTGAAAGGGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGTGCTTCGCCCTTCGGGTCGTGGCCTGCGGCAACGCTCTCTCGCTGACGCTCGGGCCGAACCTTGGGGGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                                                                                                                                                                                          <  1:748616/278‑1 (MQ=14)
    gctcagtccGAACCTTGGTCGTAGCTTATCATCCGTGCCCGCTTGGGCAGAGTATGTGATTGCGTCTTCGTTAGAGATTTCCGGGGCTTTAGAAAGGGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCGTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGGGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATGCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                                                                                                                                                                                          <  2:359781/281‑1 (MQ=14)
    gcTCCAGTCGAACCTTGGTCGAAGCTTCTCCTCCTTCCCCGTTTGGGCATAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCGGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                                                                                                                                                                                          <  2:403200/289‑1 (MQ=255)
        aaGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAgg                                                                                                                                                                                                                                                      >  2:652275/1‑289 (MQ=255)
         aGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTCAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGa                                                                                                                                                                                                                                                     <  2:405164/289‑1 (MQ=255)
               aCCTTGGTCGAAGCTTCTCATCCTTCCCCGGTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACt                                                                                                                                                                                                                                                <  2:245962/288‑1 (MQ=255)
                               ctcATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTGGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgccgttg                                                                                                                                                                                                                               >  2:610044/1‑283 (MQ=255)
                                      ttCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGAGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCg                                                                                                                                                                                                                       >  2:692572/1‑289 (MQ=255)
                                       tCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGCATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGGGGTGGGGGAAGGATTACTCAGCGCTTCGCGCTTCGCCCTTctg                                                                                                                                                                                                                      >  2:540310/1‑287 (MQ=255)
                                       tCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgag                                                                                                                                                                                                                                                                                             <  2:60020/219‑1 (MQ=255)
                                       tCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgag                                                                                                                                                                                                                                                                                             >  1:60020/1‑219 (MQ=255)
                                                tggcggagggtgtggtggcgtagacgcttGAGAATTTCTGGTGTCTGGAAGGGGATGGGGTGGTGGTAGGGATAATCAGGGCGGCGCGCCTTGTCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGGCGAACCTTGGTCGAGAGTTCTCATCCTTCCCTGTATTGGCAGGAGTTTGGTTTGTGGATTCGCTTGAGAGTTTGGGGACTTTTGGAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGcc                                                                                                                                                                                                             <  1:126212/261‑1 (MQ=2)
                                                gggCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGTCAGCGCTGCGCGCTTCGCGCTTGAGGTCGTTGCTGGCGGCAACGCCCTCTCGCGGACGCCCGCGGCGGACATTCGTCGCAGCTTCACATCCTTCCCGCCATCGGCGGTATGTTTGGTTTCGAATTCGGTGGCAGGTTCAGGGACTTTTGAAAGATATGGTGGTGGGCGAAGGATTCGTCACTGCTGCGGGCGTAGCGTTTAGGgttgttttc                                                                                                                                                                                                             >  1:33070/1‑286 (MQ=255)
                                                   cAGAACATTTGAGTGCGGATTCGTTTGAGATTGCCGGGGCTTGTGAAAGTGATGGTGTTGGGGGAAGCATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTg                                                                                                                                                                                                           <  2:120613/288‑1 (MQ=255)
                                                        tcatggATTGCGGCGTCGTTTGAGGATTAGGGTGCTTTGGGAAGGGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGATCGAGACGAACCTTGGGCGTAGCGTCTGATCCTTCCCCGCTGGGGCAGAATATTTGATTGGGGATTCGCATGAGAGTTCGGGGACTTTTGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCaa                                                                                                                                                                                                     <  1:590292/284‑1 (MQ=11)
                                                                tGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTTGTTGCCTGCGGCAACGCTCTCCGGCTGTCGCTCGCGCCGAACCTTAGTCGAAGCGTCTCAACCTTCCCCGCATGGGCAGAATAGTGCATTGCGGATTCGCTTGGTAGTTGAGGGAGTTTTTAAAGTGATGGTGGTGGGGGAAGGCTTCCTCAGTGCTACGCCCGTCGCCtttactgccgtggcgttgccgacagcactc                                                                                                                                                                                              >  1:824639/1‑261 (MQ=14)
                                                                tGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGTGCTTCACCCTTCAGGTACTTTCTTGTGGCAACGcctct                                                                                                                                                                                              >  2:593763/1‑284 (MQ=255)
                                                                          tgtgtgaATGCGGGGGGTGTTGAAAGTGAGCGTGGTGGGGGAAGTGTTCGGCAGCGGGGGCGGTGTGTCCCTGAGGGTCGTGGGGTGGGGTATCGCTTTCTCGCCGCCGCTCGGGCCGAAACTGAGTGGTATCTTCGCAGCCGTCCCCGCATGGGCATCATATTTGATTGCGGATTAGTTAGAGGGTTCAGGTACTTTTGGAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCt                                                                                                                                                                                   <  1:467166/284‑1 (MQ=9)
                                                                                     cGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCCGCGCTTCGCGCTTCGCCCTTCGGGTCGTTTCCTTCGGCAACGCTCCCTCGCTGGCCCTCGAGTCGaa                                                                                                                                                                          >  2:690882/1‑287 (MQ=255)
                                                                                          tttaaggaaGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTTGCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCCTTTGAAAGTGATGTTGGTGGGGGAAGGATTAATCAGCGCTGGGCGCTTAGACCCTCCGGTCGTTGCCGGCGGCAGCACCTTCTCCCTGgcgtttcatctggcgctctg                                                                                                                                                                    >  1:405164/8‑271 (MQ=14)
                                                                                                       tggtggtggtggAAGGCGTCCGCAGCGCCGCACGCCTCGTCGTCGGCGTCGTCGCCTGCGGCTACGCTCCAGCGCTGGCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGAGCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAt                                                                                                                                                      <  1:652275/289‑1 (MQ=14)
                                                                                                             tggGGGAAGGATTCCGCAGCGCGGCGCGCTTCGCCGTGCGGGGCGTTGCCTGCGGCAACGCTTTGTTGCTGATGCTCAGGCCGAACCTTAGTCGAAGCCTCTCATCCTTCCCGGCATGGGCAGAAGATTTGATGGCGGATTAGTTTGAGAGTTCGGTGACTTTTGGAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTc                                                                                                                                                 <  1:643243/288‑1 (MQ=11)
                                                                                                                        ttACTTAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                                                                    <  2:671840/289‑1 (MQ=255)
                                                                                                                         ttacagggtcctggtcacgcgtactgtCGGGACGTGGCGTACGGCAACGCTTACAGACGGGCGATCCAGACGGAGCTAAGTCCAAGATTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTGGAGAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCGTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCGCGAGTCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                                                                    <  1:540310/263‑1 (MQ=11)
                                                                                                                           cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTGGGGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGTGGTTGGCGCTTCGGGTCGTTGCTGGGGCCAACGCTCTCTCGCGGACGCGCGTGTCGAACACTGGTGGCGCCTTCTCCGCCTTCCCGCGCTGGGCAgca                                                                                                                                 >  1:609452/1‑287 (MQ=14)
                                                                                                                              aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAatat                                                                                                                              >  2:310456/1‑289 (MQ=255)
                                                                                                                                   tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGTTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCCCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAAt                                                                                                                         >  2:769294/1‑289 (MQ=255)
                                                                                                                                     cgctcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTTCCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGTGCAGCCTATCTGCTTGCGGATTCGCTTCCGATTTCCGGGCCTTTTGTAACTGATGGTGGTGCGGCAAGGATTACTCATCTTTCATCCTGTCGTCACTCCCTTCGTTGCCTCCGGCTGCGCTTTCTCGCTGTCCCTCGACGCTAGCCTTGTCCGCCGCTTCTCATTCTTCCCTGcgtagcttgcacagttactt                                                                                                                        >  1:760503/1‑269 (MQ=2)
                                                                                                                                     cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGGGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCCTTGCCGGCGCCAATGCCCTCTCGCTGGCGCCCAGCTGGCACCTTGGTCGAAACTTCCACTCCTTCCCCGCGTGGGCCGCAGATTTaacta                                                                                                                       >  1:334211/1‑286 (MQ=14)
                                                                                                                                                   tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTggg                                                                                                          <  2:643205/289‑1 (MQ=255)
                                                                                                                                                            tGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCCCTCTCGCCGGCGCTCGGGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTTAgg                                                                                                 >  2:146825/1‑289 (MQ=255)
                                                                                                                                                             cccTGCGGAGACGCTCTCTTGGGCACGCTCAAGCTGATACATGGTCCAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGTTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGGGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAggg                                                                                                <  1:692572/288‑1 (MQ=14)
                                                                                                                                                                      aaCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGa                                                                                        <  2:670573/288‑1 (MQ=255)
                                                                                                                                                                           tctctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGa                                                                                  <  1:310456/289‑1 (MQ=255)
                                                                                                                                                                                 gCTGACGCTCGAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCCGAATATTTGGTTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtgg                                                                            <  2:114265/289‑1 (MQ=255)
                                                                                                                                                                                         tCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGAGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCCGGGACTTTTGAAAATGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                                                                    <  2:570612/289‑1 (MQ=255)
                                                                                                                                                                                                  aCCTTAGTCGAAGCTTCTCATACTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGa                                                           <  2:824639/288‑1 (MQ=255)
                                                                                                                                                                                                             aGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTACCCGCTACGTCCTTCGCGGTTCGGTTCGTGTCTTGCGCCTATGCTCGCTCGCTGACGCGTGTGTTATACCCTGTTTGACCTTTCGACTCCTTCCCCTGTTGTGCAGTATATTTTACCCCTGATCCTGTGGCAAGTGCAGGGACTCATGAACGTGACGCTGGGCGGGGCAGGGTTCGTAcgtacaggtc                                                >  1:133923/1‑279 (MQ=2)
                                                                                                                                                                                                              gCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTTGTGGGGGAAGGATTCGAACCTTCGAAGTCg                                              >  2:10928/1‑289 (MQ=255)
                                                                                                                                                                                                              gCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCCCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGTATTCGTTGTGAAGTTCAGTGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGCCg                                              >  2:206481/1‑289 (MQ=255)
                                                                                                                                                                                                                     gtgtataatggacagggaggAATCTTTGATTGCGGAGTCGCTTGATAGTTGAGGAGGTTGTGAAAGTGATGGGGGTGGGGGAAGGGTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCGCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACgg                                       <  1:483704/270‑1 (MQ=17)
                                                                                                                                                                                                                     aTCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGGTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACgg                                       <  2:305317/289‑1 (MQ=255)
                                                                                                                                                                                                                     aTCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACgg                                       <  2:819689/289‑1 (MQ=255)
                                                                                                                                                                                                                       gcTTCCCCGCATGGGCAGTATATTGTATTGTGGATTCGCTTGGTAGTTCAGGGACTTTTGAAAGGGATGTGGGTGGGGGAGGGATTATCAGGTGCTGCTCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCa                                     <  1:146825/288‑1 (MQ=17)
                                                                                                                                                                                                                             ccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGCGGTCCGGGACTTTTGTAAGGGAATGAGGTGGGGGGAGGGTTCGGAACCtcggaatccatcgtgggatgttttc                               >  2:265586/1‑265 (MQ=14)
                                                                                                                                                                                                                                  agggCAGGATTATTTGTTTGCGGATTCGTTTGGGAGTTCGGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACa                             <  1:593763/286‑1 (MQ=255)
                                                                                                                                                                                                                                     gCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTg                        <  2:334211/288‑1 (MQ=255)
                                                                                                                                                                                                                                     gCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTg                        <  2:270353/288‑1 (MQ=255)
                                                                                                                                                                                                                                                tGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGc             <  2:609452/289‑1 (MQ=255)
                                                                                                                                                                                                                                                tGATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTGCACCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGc             <  2:133923/289‑1 (MQ=255)
                                                                                                                                                                                                                                                tGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGTTCTGCTCGCTTCGCCCTTCGGGTCGTGGCCTGCGGCGATGCTCTCTCGCTGTCCCTCGAGTCGAAAATTGGTGGAATCTTCTCATCCTTCCACGCATGGGCAGAATATGTTATTGCGGGTCCGTTGGGAAGTTCAGAGACTGTTGAAGGTGGTGAGGGTGGAGAATGGAGTTGTCTCTTCGAGCTCTAACACTGGAGATTTACAGTATGCTCCATTGGGc             >  1:855772/1‑289 (MQ=9)
                                                                                                                                                                                                                                                tGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGACTATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGGGGGGGAAGGATTCGAACCTTTGAAGTCGATGATGGCAGGTTTACa                             >  2:483704/1‑273 (MQ=255)
                                                                                                                                                                                                                                                    tGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTGGAAAGTGATGGTGGGGGGGGAAGGAGTCGAACCTACGAAGTCGATGACGGGAGATTTACAGTCTGCTCCCTTTGGCCGc          >  2:502602/1‑288 (MQ=255)
                                                                                                                                                                                                                                                      cGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCCACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCCCATCCTTCCCCGCATGGGCCGAATATTTACTTGCGGCTTTTGTGGGAAGTTCGGGGACGTGTGAAAGTGATGGTGGTGGGGCAAGGCTTCCAACCTTCGAAGTCGATGACGTAAGATTAACAGTCTGCCCCTTTTTGCCGCTCg       >  1:765585/1‑289 (MQ=17)
                                                                                                                                                                                                                                                      cGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCg       >  2:42138/1‑289 (MQ=255)
                                                                                                                                                                                                                                                      cGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCCCg       >  2:139736/1‑289 (MQ=255)
                                                                                                                                                                                                                                                       ggATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTTCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGCTGGTGGTGTGGGAAGGATTCGAACCTTCGAAGTCGTTGACGGCAGCTTTACCGTCTGCTCCCCTTGGCCgaccgg      >  2:364269/1‑284 (MQ=255)
                                                                                                                                                                                                                                                       ggATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCgg      >  2:479307/1‑289 (MQ=255)
                                                                                                                                                                                                                                                         aTTGGCTTGGGAGTTCCGGGATGTTTGAAAGTGATGTTGGTGGGGGTAGGATTACTCAGCGCTGGGCGCTTCGGCCATCGAGTCGTTGCATGTGGTAACGCGCGCTCGCTGGCTATCGGGTCGAGCCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCggg     <  1:10928/288‑1 (MQ=18)
                                                                                                                                                                                                                                                         aTTAGCTTGAGAGTTCGGTGACTTGAGAAAGTGATAGTGGTGGGGGAGGGATTAATCAGGGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGGATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCggg     <  2:49055/288‑1 (MQ=18)
                                                                                                                                                                                                                                                           tCGCTTGAGAGTTCAGGGACTTTTGAAAGGGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAAc  <  2:674421/289‑1 (MQ=255)
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                
TGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAAC  >  NC_000913/1286785‑1287322

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