Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,711,636 A→C 18.0% T407P (ACC→CCC)  nadB → quinolinate synthase, L‑aspartate oxidase (B protein) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,711,6360AC18.0% 164.6 / 3.9 74T407P (ACC→CCC) nadBquinolinate synthase, L‑aspartate oxidase (B protein) subunit
Reads supporting (aligned to +/- strand):  ref base A (37/21);  new base C (5/11);  total (42/32)
Fisher's exact test for biased strand distribution p-value = 2.51e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCTCACACAAGAACCGGTACCGATTGTGCCTGCTGCACATTATACCTGCGGTGGTGTAATGGTTGATGATCATGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCGCCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCC  >  NC_000913/2711429‑2711855
                                                                                                                                                                                                                 |                                                                                                                                                                                                                           
gatttcccacgaagtccggttccggttGTGCGTGGAGGACATTATCCCTGGGGTGGTGTAATGGTTGGTGGTCATGGGCGTACGGACGTCGAGGGCTTGGATGCCATTGGCGAGGTGAGTTATACCGGCTTACCCGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATAAg                                                                                                                                                                                              <  2:1021308/224‑1 (MQ=255)
gATCTCACACAAGAACCGGTACCGATTGTGCCTGCTGCACATTATACCTGCGGTGGTGTAATGGTTGATGATCATGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGGGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTTTGCCCACAACATCAg                                                                                                                                                                                              >  1:1067032/1‑239 (MQ=255)
      acacaAGAACCGGTACCGATTGTGCCTGCTGCACATTATACCTGCGGTGGTGTAATGGTTGATGATCATGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGGGGCGGACGATATCCCCAGACGGAAGCCTTTTGcctgtctcttataaacat                                                                                                                                                                                         >  1:454263/1‑222 (MQ=255)
      acacaAGAACCGGTACCGATTGTGCCTGCTGCACATTATACCTGCGGTGGTGTAATGGTTGATGATCATGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCAAGGCCTCGAATTCATTGCTGGAGTGTCTGGGCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCCCGACATCCGTACGtt                                                                                                                                                                                        >  1:165353/1‑239 (MQ=255)
       cacaAGAACCGGTACCGATTGTGCCTGCTGCACATTATACCTGCGGTGGTGTAATGGTTGATGATCATGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGGCGGCGGCGGGAGATATCACCAGACGTATGCCTTTTGCCCCCGCCATCAGTACGtt                                                                                                                                                                                        >  1:547259/1‑238 (MQ=255)
           aGAACCGGTACCGATTGTGCCTGCTGCACATTATACCTGCGGTGGTGTAATGGTTGATGATCATGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGGCTCTGGCTGTTCGGCGGCGGAAGATATCACCAGACGTGTGCATTTTGCCCCAGACATCAGTACGTTAccg                                                                                                                                                                                    >  1:567157/1‑238 (MQ=255)
            tgaccggtgcgtgttGTGCCGGCTGCACTTTATACCGGGGGGGGTGTAATGGTTGATGATCGTGGGCGTACGGGCGTCGGGGGTTTTTATGCCATTGGCGAGGTGAGTTAGCCCGGTTTACACGGCGCTAACCGCAGGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCCCCAGACGTATGCCTTATGCCCACGACATCAGTACGTTAccgcc                                                                                                                                                                                  <  2:1067032/226‑1 (MQ=255)
            gAACCGGTACCGATTGTGCCTGCTGCACATTATACCTGCGGTGGTGTAATGGTTGATGATCATGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTAccgcc                                                                                                                                                                                  >  1:841874/1‑239 (MQ=255)
              cccGGTCCCGTTTGTGCCTGCTGCCAATTATCCCTGGGGTGGTGTAATGGTTGATGATCATGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGt                                                                                                                                                                                <  2:841874/238‑1 (MQ=255)
              acgggTACCTTTTGTGCCGGGTGCACATTACACCGGCGGTGGGGTAATGGTTGATGATCAGGGGGGTACGGACGGCGGGGGTTTTATGGCCAGTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGt                                                                                                                                                                                <  2:124534/236‑1 (MQ=255)
                  gTACCGATTGTGCCTGCTGCACATTATACCTGCGGTGGTGTAATGGTTGATGATCATGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGCCATCAGTACGTTACCGCCGTgggg                                                                                                                                                                            >  1:733374/1‑238 (MQ=255)
                        aTTGTGCCTGCTGCACATTATACCTGCGGTGGTGTAATGGTTGATGATCATGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGGAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCCTGGGATGaaa                                                                                                                                                                       >  1:577117/1‑238 (MQ=255)
                             gCCTGCTGCACATTATACCTGGGGTGGTGTAATGGTTGATGATCATGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCcgc                                                                                                                                                                  <  1:705179/238‑1 (MQ=255)
                                 gctgcACATTATCCCTGCGGTGGTGGAATGGTTGATGATCATGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTga                                                                                                                                                             <  1:725615/239‑1 (MQ=255)
                                    gacagttATACCTGCGGTGGTGTAAGGGTTGATGATCATGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATCCCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTgag                                                                                                                                                            <  1:133520/232‑1 (MQ=255)
                                            taCCTGCGGTGGGGGAATGGTTGATGATCAGGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCAc                                                                                                                                                                                                      <  2:270216/187‑1 (MQ=255)
                                                tGCGGTGGTGTAATGGTTGATGATCATGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGCCGAACg                                                                                                                                              >  1:156122/1‑239 (MQ=255)
                                                tGCGGTGGTGAAATGTTTGACGATCATGGGGGTACGGACGTCGAGGGCTTGTATGCCATGGGCGAGGTGAGTTAACCCGGCTAACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACg                                                                                                                                              <  1:685606/239‑1 (MQ=255)
                                                  cGGTGGTGTAATGGTTGATGATCATGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCAc                                                                                                                                                                                                      >  1:370368/1‑181 (MQ=255)
                                                  cGGTGGTGTAATGGTTGATGATCATGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCAc                                                                                                                                                                                                      <  2:370368/181‑1 (MQ=255)
                                                  cGGTGGTGTAATGGTTGATGATCATGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAACCGCGTTGAGAACCCTGCCGAACggg                                                                                                                                            >  1:73398/1‑239 (MQ=255)
                                                  cGGTGGTGTAATGGTTGATGATCATGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGGCTATGGCTGGTCGGCGGCGGAAGATCTCACCCGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGGGGGATGAAAGCCGCGTTGAAAACCCCGACGAAAggg                                                                                                                                            >  2:1099940/1‑239 (MQ=255)
                                                    gtggtgTAATGGTTGATGATCATGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTCTGGCTGGGTGGCGGGGGCAGCAATCACCGAACGAATGCCCTATGCCCACGAACTCAGTACTGAACCGCCGGGGAGTGAAAACCCCGTGGAAAACACCAACGAACggggg                                                                                                                                          >  2:408687/1‑237 (MQ=255)
                                                    gtggtgTAATGGTTGATGATCATGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACgggga                                                                                                                                          >  1:977828/1‑237 (MQ=255)
                                                    gtggtgTAATGGTTGATGATCATGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACgggaa                                                                                                                                          >  1:873082/1‑237 (MQ=255)
                                                    gtggtgTAATGGTTGATGATCATGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGAAATACCCAGACGGATCGCTTATGCCCCCGGCATCCGTACGTTTACCCCGGGGGATGAAAACCCCGCTGGGAACCCTCACGACCgcg                                                                                                                                            >  2:133520/1‑235 (MQ=255)
                                                    ggtgggtattggtgtagtatctgggggGTCCGGGCGTGGGGGGTTGGTAGGCCGTTGGCGAGGTGAGTTATACCGGCTGACAGTGCGCTAACCGCATGGCCTCGAATTCATTGCGGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATAGCACCAGACGTATGCCTCATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGgt                                                                                                                                           <  2:646336/215‑1 (MQ=255)
                                                      ggtgTAATGGTTGATGATCATGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTa                                                                                                                                                                                            >  1:270216/1‑187 (MQ=255)
                                                       gtgtAATGGTTGATGATCATGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGCCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGAACGTTACCGCCGTCGGATGAAAGCCGCGTTTAGAACCCTGACGAACGGttagta                                                                                                                                       >  1:348168/1‑239 (MQ=255)
                                                        tgtAATGGTTGATGATCATGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGGGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTCTGGCTGGTCGGCGGCGGAAGATATCCCCAGACGTATGCCTTATGCCCCCCGCCTCAGTACCTTACCCCCCGGGGCTGAAACCCGCGGGGAGAAACCAGACAAACGGGAAGaaa                                                                                                                                      >  2:574221/1‑239 (MQ=255)
                                                             ttgTTGAGGGTCATGGGGGTGCGGGGGTGGCGGGGTGGTATTGCGTTGGCGGGGTTGGTTATACCGGTTTTACCGGCCGTACCCGCATGGCCTCGAATTCATTGCGGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATTTCCCCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAg                                                                                                                                 <  2:577117/237‑1 (MQ=255)
                                                             tGGTTGATGATCCGGGGCGTACGGGCGTCGGGGGTTTGTATGCCATTGGCGAGGGAAGTTATGCCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAg                                                                                                                                 <  1:712617/239‑1 (MQ=255)
                                                                 tgatgaTCATGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGCCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGTACCCTGACGAACGGGAAGTAATTACGCATa                                                                                                                             >  1:499263/1‑239 (MQ=255)
                                                                      atCATGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACAt                                                                                                                                                                                                 <  1:452539/166‑1 (MQ=255)
                                                                        cgggggCGTTCGGAAGGGGGGGGCTTGTTTTCCATTGGCGAGGTGGGTTAGACCGGTTTACACGGCGTAAACCGCATGGCCTGGAATTGATGGCTGGAGTGTCTGTTCTATGGCTGGTCGGGGGCGGAAGATATCACCAGAGGTATGCCTTATGCCCAGGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCa                                                                                                                      <  2:73398/236‑1 (MQ=255)
                                                                         aTGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATGCCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCCTTGCTGGGGTGTCTGGTCTATGGCTGGTCGGCGGCGGACGATATCCACAGACGGATGCCTTTTGCCCACGACATCAGGACGTTACCGCCGTGGGAGTCAGGCCGCGTTGCGAACCCTGACCACCGGGGAGTAATTCAGCaaaaaagggca                                                                                                                     >  2:191620/1‑232 (MQ=255)
                                                                         aTGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGGGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGCCTATGGCTGGCCGGCGGCGGAAGATATCACCAGACGTATTCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGGAGAAATTCAGCATAACTGGCAc                                                                                                                     >  1:167959/1‑239 (MQ=255)
                                                                             gagtgaggccgccgggggCTGTTATGCCATTGGCGAGGGGGTTTATCCCGGTTTACACGGCCCTACCCGCATGGCCTCGAATACATGGCTGGAGGGTCTGGCCTATGGCTTGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTACGCCCACGACGTCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAg                                                                                                                  <  2:671371/223‑1 (MQ=255)
                                                                               ttgcggcgtggggagggttttttGCCATTGGGGGGGTGAGTTATGCCGGCTTCGCGGGCGCTAACGGGAGGGCCTCGAATCCATTGATGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCCCCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTa                                                                                                               <  2:733374/217‑1 (MQ=255)
                                                                                tACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGCCGAACGGGTAGTAATTCAGCATAACTGGCACGAGcacg                                                                                                              >  1:1041234/1‑236 (MQ=255)
                                                                                  cGGGCGTCGGGGGTTTTTATTCCTTGGGGGGGGGTAGTTATACCGGGTTACCCGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGt                                                                                                            <  2:977828/239‑1 (MQ=255)
                                                                                       gTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGAGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGCAGTAATTAAGCATAACTGGCACGAGCTACGCCTGtt                                                                                                       >  1:242387/1‑239 (MQ=255)
                                                                                               cTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGAAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTggg                                                                                               >  1:545968/1‑239 (MQ=255)
                                                                                                 tGTATGCCATTGGCGAGGTGAGTTATACCGGTTTACACGGTGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCCCCAGACGTa                                                                                                                                                                                                                    <  1:782384/120‑1 (MQ=255)
                                                                                                 tGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCCCCAGACGTa                                                                                                                                                                                                                    >  2:782384/1‑120 (MQ=255)
                                                                                                 tGTATGCCATGGGCGAGGGGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGAt                                                                                             <  1:628555/239‑1 (MQ=255)
                                                                                                         cATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCCTAACTGGCCCGAGCTACGTCTGTTTATTTGGGATTACGTTg                                                                                      >  1:567707/1‑239 (MQ=255)
                                                                                                           gtGGCGAGGTGAGTTATACCGGCTTACGCGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGc                                                                                    <  1:722618/238‑1 (MQ=255)
                                                                                                                   gTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGAAATTCAGCATACCTGGCACGAGCAACGTCTGTTTATGTGGGATTACGTTGGCATTGGgcg                                                                            >  1:1102820/1‑239 (MQ=255)
                                                                                                                    aaaagggacagggCTTACCCGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTgcgc                                                                           <  1:184797/228‑1 (MQ=255)
                                                                                                                           taCCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCCCAACGa                                                                    >  1:832352/1‑239 (MQ=255)
                                                                                                                                          cgcTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCCCGTTGAGAACCCTGGCGAACGGGTAGTAAATCAGCATAACTGGCCCGAGCCACGTCGGTTTTTGTGGGTTTACGTTGGGCTTTTTCGCAAAAAAAAGCCCCGGGAAggcg                                                     >  2:287165/1‑239 (MQ=255)
                                                                                                                                          cgcTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCCCGTTGAGAACCCTGACGAACCGGTGAGAATTCAGCCTAACTGGCCCGAGCTACGTCTGTTTATGGGGGTTTACGTTGGACTTGTGCCGCCAACAAAgc                                                                 >  2:184797/1‑227 (MQ=255)
                                                                                                                                           gcaaaCCGAATGGCCTCGAATTCATGCGTGGAGTGCCTGGTCTAGGGCTGGTCGGCGGCGGAAGATATCCCCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAAcgcg                                                     <  1:999382/235‑1 (MQ=255)
                                                                                                                                               aCCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGt                                                                                                            <  2:127854/179‑1 (MQ=255)
                                                                                                                                               aCCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGt                                                                                                            >  1:127854/1‑179 (MQ=255)
                                                                                                                                                ccGCATGGCGTCGAATTCTTTGCTGGGGTGTCTGGTCTGTGGTTGGTGGGCGGCGGAAGATAGCACCGGACTTATGCCTTATGCCGACGACAGTAGTTCGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCGCCCTg                                                <  2:156122/238‑1 (MQ=255)
                                                                                                                                                 cGCTTGGCCTTGAATTCATTGCTGGGATGTCGGGGTTATGGCTGGGCGGCGGCGGAAGATATCCCCAGACGTGTGCCTTATGCCCACGACATCAGTACGTTACCGACGTGGGAGGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCGCCCTgcg                                              <  2:567707/239‑1 (MQ=255)
                                                                                                                                                    aTGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGCCGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCCCAACGAAGCGCCtctgtcccttatacacat                                            >  1:883132/1‑221 (MQ=255)
                                                                                                                                                    aTGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCCCAACGAAGCGCCTGGAACGCGCCCTgcggcg                                           >  1:323291/1‑239 (MQ=255)
                                                                                                                                                                       tGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATTTGGGATTACGTTGGCATTGTGCGCCCAACGAAGCGCCTGGAACGCGCCCTGcggcgggaaacaccgccccaaaaa                        >  2:660495/1‑221 (MQ=255)
                                                                                                                                                                           gggtgtgttcttgggcggggcggGGGCGGAAGATTCCCCCAGACGTTTGCCTTATGCCCACGAAATGGGTGCGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCGCCCTGCGGCGGATAACCATGCTCCAACAAGaaa                    <  2:952075/226‑1 (MQ=255)
                                                                                                                                                                            ttgttGGTTTATGGCGGGTCGGGGGCCGAAGTTATCCCCAGACGTTTGCCTTATGCCCACGGACTCATTCCTTTACCGGCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCGCCCTGCGGCGGATAACCATGCTCCAACAAGaaa                    <  1:1081393/234‑1 (MQ=255)
                                                                                                                                                                              tCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCCGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCGCCCTGCGGCGGGTAACCATGCTCCAACAAGAAATaa                 >  1:126576/1‑238 (MQ=255)
                                                                                                                                                                              tCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTACGCACAACGAAGCGCCTGGAACGCGCCCTGCGGCGGATAACCATGGTCCAACCAGAAATAg                 >  1:517739/1‑239 (MQ=255)
                                                                                                                                                                              tCTGGTCTAGGGCGGGTCGGGGGCGGAAGATATCCCCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCGCCCTGCGGCGGATAACCATGCTCCAACAAGAAATa                  <  1:660495/238‑1 (MQ=255)
                                                                                                                                                                                tGGTCTATGGGTGGTCGGCGGCGGAAGTAACCCCCAGACGTATGCCTTATGCCCACGACATCAGTACGGTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCGCCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGAc               <  1:1099940/239‑1 (MQ=255)
                                                                                                                                                                                 ggTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGACGCGTTGAGAACCCTGACGACCGGGTAGTAATTCAGCATAACTGGCACGAGCTACGGCGCTTTATGTGGGATTAGCTGGGCATTGTGCGCACAACGAAGCCCCTGGAACGCGCCCTGCGGCGGATAAACAAGCCCCCACaaaaaaaaaaaa              >  2:539975/1‑233 (MQ=255)
                                                                                                                                                                                       tGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCGCCCTGCGGCGGATAACCATGCTCCCACAAGAAATAGACGAATAtt        >  1:716747/1‑239 (MQ=255)
                                                                                                                                                                                       tGGCTGGGCGGCGGCGCAAGATATCACCAGACGTTTGCCTTATGCCACGGAAATCAGTATGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCGCCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATAtt        <  1:852913/239‑1 (MQ=255)
                                                                                                                                                                                       ggggtggtggGCGGCGGAAATATTCCCAAGGCGTATGCCTTATGCCCACCACCTAATTACGTTACCGCCGTGGGATAAAACCCGCGTTTAGAACCCTGACGAACGGGTAGTAAGTCAGCATAACTGGCACGAGCTACGTCTGGTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCGCCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATAtt        <  2:242387/235‑1 (MQ=255)
                                                                                                                                                                                        ggccagTCGCCGGAGGAAGATATCACCCGACGTATGCCTTATGCCCACGACCTCCGGACGTTACCGCCGTGGGGTGGAAGCCGCGTAGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGTAACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCCCAACGAAGCGCCTGGAACGGGCCCTGGTGCGGATAACCATGCTCCAAAAAAAAAAAAACGCACAtt        >  2:652721/6‑238 (MQ=255)
                                                                                                                                                                                        ggCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCGCCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAatatct       >  1:952075/1‑237 (MQ=255)
                                                                                                                                                                                             gTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCGCCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGccc  >  1:298771/1‑239 (MQ=255)
                                                                                                                                                                                                                 |                                                                                                                                                                                                                           
GATCTCACACAAGAACCGGTACCGATTGTGCCTGCTGCACATTATACCTGCGGTGGTGTAATGGTTGATGATCATGGGCGTACGGACGTCGAGGGCTTGTATGCCATTGGCGAGGTGAGTTATACCGGCTTACACGGCGCTAACCGCATGGCCTCGAATTCATTGCTGGAGTGTCTGGTCTATGGCTGGTCGGCGGCGGAAGATATCACCAGACGTATGCCTTATGCCCACGACATCAGTACGTTACCGCCGTGGGATGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAAGCGCCTGGAACGCGCCCTGCGGCGGATAACCATGCTCCAACAAGAAATAGACGAATATTACGCCC  >  NC_000913/2711429‑2711855

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: