Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,394,935 C→T 5.2% R290R (CGC→CGT)  nnr → bifunctional NAD(P)H‑hydrate repair enzyme; C‑terminal domain ADP‑dependent (S)‑NAD(P)H‑hydrate dehydratase and N‑terminal domain NAD(P)H‑hydrate epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,394,9350CT 5.2% 252.5 / 3.0 77R290R (CGC→CGT) nnrbifunctional NAD(P)H‑hydrate repair enzyme; C‑terminal domain ADP‑dependent (S)‑NAD(P)H‑hydrate dehydratase and N‑terminal domain NAD(P)H‑hydrate epimerase
Reads supporting (aligned to +/- strand):  ref base C (37/36);  new base T (2/2);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

AAAACCGCGTCGCCCGACTTCGCATAAAGGCGATCACGGGCGGCTGGTAATTATCGGTGGCGATCACGGCACGGCGGGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGCCAGCAAGAGTGGGGGAAAAAAGCA  >  NC_000913/4394791‑4395082
                                                                                                                                                |                                                                                                                                                   
aaaaCCGCGTCGCCCGACTTCGCATAAAGGCGATCACGGGCGGCTGGTAATTATCGGTGGCGATCACGGCCCGGCGGGGGCTATTCGTATGACGGGGGAAGCGGGGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGt                                                                                                                                                  >  6:49492/1‑148 (MQ=255)
      gcgTCGCCCGACTTCGCATAAAGGCGATCACGGGCGGCTGGTAATTATCGGTGGCGATCACGGCACGGCGGGGGCTATTCGTAAGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAAc                                                                                                                                            <  1:151904/148‑1 (MQ=255)
        gTCGCCCGACTTCGCATAAAGGCGATCACGGGCGGCTGGTAATTATCGGTGGCGATCACGGCACGGCGGGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGCCCCGCAGTGAAAACAt                                                                                                                                          >  6:66683/1‑148 (MQ=255)
           gCCCGACTTCGCATAAAGGCGATCACGGGCGGCTGGTAATTATCGGTGGCGAGCACGGCACGGCGGGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTg                                                                                                                                        <  2:199633/147‑1 (MQ=255)
                   tCGCATAAAGGCGATCACGGGCGGCTGGTAATTATCGGTGGCGATCACGGCACGGCGGGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTGTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCgctgc                                                                                                                               >  2:217504/1‑148 (MQ=255)
                           aGGCGATCACGGGCGGCTGGTAATTATCGGTGGCGATCACGGCACGGCGGGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCa                                                                                                                       >  5:2906/1‑148 (MQ=255)
                             gCGATCACGGGCGGCTGGTAATTATCGGTGGCGATCACGGCACGGCGGGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCAc                                                                                                                      >  5:198622/1‑147 (MQ=255)
                                 tCACGGGCGGCTGGTAATTATCGGTGGCGATCACGGCACGGCGGGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGa                                                                                                                            >  3:224992/1‑137 (MQ=255)
                                 tCACGGGCGGCTGGTAATTATCGGTGGCGATCACGGCACGGCGGGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCg                                                                                                                 >  3:66610/1‑148 (MQ=255)
                                       gcggcTGGTAATTATCGGTGGCGATCACGGCACGGCGGGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTg                                                                                                           <  4:293001/148‑1 (MQ=255)
                                        cggcTGGTAATTATCGGTGGCGATCACGGCACGGCGGGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGa                                                                                                          <  6:198622/148‑1 (MQ=255)
                                         ggcTGGTAATTATCGGTGGCGATCACGGCACGGCGGGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGG‑TTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATg                                                                                                        <  1:248971/148‑1 (MQ=255)
                                            tGGTAATTATCGGTGGCGATCACGGCACGGCGGGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGt                                                                                                      >  5:219934/1‑148 (MQ=255)
                                              gTAATTTTCTGTGGCTATCACGGCTCGGCGGGGGCTATTCGTATTACGGGGTAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTg                                                                                                     <  4:84514/147‑1 (MQ=255)
                                               tAATTATCGGTGGCGATCACGGCACGGCGGGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCa                                                                                                   >  7:134275/1‑148 (MQ=255)
                                                 aTTCTCGGTGGCGATCACGGCACGGCGGGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATg                                                                                                 <  8:218823/148‑1 (MQ=255)
                                                  ttATCGGTGGCGATCACGGCACGGCGGGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGAATg                                                                                                 >  6:222674/1‑147 (MQ=255)
                                                   tATCGGTGGCGATCACGGCACGGCGGGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGaa                                                                                               >  8:115306/1‑148 (MQ=255)
                                                   tATCGGTGGCGATCACGGCACGGCGGGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGaa                                                                                               >  3:29617/1‑148 (MQ=255)
                                                    ttCGGGGGCGCTCACGGCTCGGCGGGGGCTATTCGTCTGACGGGGGCCGCGGCGCTGCGTGCTGGTGCGGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGaa                                                                                               <  3:92451/146‑1 (MQ=255)
                                                    aTCGGGGGCGATCACGGCTCGGAGTGGGCTATTCGTATGACGGGGGAAGCTGCGATGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGGGCATGaa                                                                                               <  4:66610/147‑1 (MQ=255)
                                                      cGGTGGCGATCACGGCACGGCGGGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACt                                                                                             >  1:226478/1‑147 (MQ=255)
                                                             gATCACGGCACGGCGGGGGATATACGTATGACGGGGCAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATgg                                                                                     <  4:29617/148‑1 (MQ=255)
                                                              aTCACGGCACGGCGGGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGa                                                                                    >  7:38662/1‑148 (MQ=255)
                                                              aTCACGGCACGGCGGGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGa                                                                                    >  8:174882/1‑148 (MQ=255)
                                                                  cggcacggcGGGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCAt                                                                                                  >  6:19022/1‑130 (MQ=255)
                                                                  cggcacggcGGGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCAt                                                                                                  <  5:19022/130‑1 (MQ=255)
                                                                  cggcacggcGGGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGActct                                                                                <  7:115306/148‑1 (MQ=255)
                                                                   ggcacggcGGGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGActctc                                                                               <  6:219934/148‑1 (MQ=255)
                                                                    gcacggcGGGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGTAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGActctct                                                                              >  8:17069/1‑148 (MQ=255)
                                                                    gcacggcGGGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGTAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGActctct                                                                              >  8:22034/1‑148 (MQ=255)
                                                                    gcacggcGGGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGActctct                                                                              >  8:71640/1‑148 (MQ=255)
                                                                      acggcGGGGGCTATGCGTATGACGGGGGCAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTAGTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTa                                                                            <  5:287350/148‑1 (MQ=255)
                                                                       cggcggGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTAc                                                                           >  7:114089/1‑148 (MQ=255)
                                                                          cggGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGa                                                                        <  5:188040/148‑1 (MQ=255)
                                                                               gCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGcc                                                                   <  2:226478/148‑1 (MQ=255)
                                                                                  aTTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTg                                                                 >  8:312670/1‑147 (MQ=255)
                                                                                    tCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGc                                                                                                                        <  8:186166/90‑1 (MQ=255)
                                                                                    tCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGc                                                                                                                        >  7:186166/1‑90 (MQ=255)
                                                                                         tGACGGGGGAAGCGGCGCTGCGTGCGGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATggg                                                          <  6:2906/147‑1 (MQ=255)
                                                                                          gACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGaaa                                                                                                                                              >  1:124193/1‑62 (MQ=255)
                                                                                          gACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGaaa                                                                                                                                              <  2:124193/62‑1 (MQ=255)
                                                                                             gggggAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGaa                                                                                               <  7:26686/106‑1 (MQ=255)
                                                                                             gggggAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGaa                                                                                               >  8:26686/1‑106 (MQ=255)
                                                                                                ggAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCt                                                                  >  8:81205/1‑132 (MQ=255)
                                                                                                ggAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCGGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCt                                                                  <  7:81205/132‑1 (MQ=255)
                                                                                                       gcgcTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGa                                           <  7:71640/148‑1 (MQ=255)
                                                                                                                     tgctggtTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATg                                                            >  6:99963/1‑117 (MQ=255)
                                                                                                                     tgctggtTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATg                                                            <  5:99963/117‑1 (MQ=255)
                                                                                                                     tgctggtTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTg                             >  3:291535/1‑148 (MQ=255)
                                                                                                                      gctggtTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATgtggtggtg                                            >  1:83868/1‑132 (MQ=255)
                                                                                                                      gctggtTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATgtggtggtg                                            <  2:83868/132‑1 (MQ=255)
                                                                                                                       ctggtTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTa                                                                            <  6:2461/99‑1 (MQ=255)
                                                                                                                       ctggtTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTa                                                                            >  5:2461/1‑99 (MQ=255)
                                                                                                                         ggtTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGc                          >  7:259850/1‑147 (MQ=255)
                                                                                                                           tttAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGCCAg                       >  1:219556/1‑148 (MQ=255)
                                                                                                                            ttAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTAcc                                                                          >  8:161837/1‑96 (MQ=255)
                                                                                                                            ttAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTAcc                                                                          <  7:161837/96‑1 (MQ=255)
                                                                                                                             tAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGt                                >  7:202691/1‑137 (MQ=255)
                                                                                                                             tAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGt                                <  8:202691/137‑1 (MQ=255)
                                                                                                                              aGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGCCAGCa                     <  4:291535/147‑1 (MQ=255)
                                                                                                                               gTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGCCAGCAag                   <  7:174882/148‑1 (MQ=255)
                                                                                                                               gTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGGCTGGGCCAGCAag                   <  5:49492/148‑1 (MQ=255)
                                                                                                                                 ccGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGCGAGCAagag                 >  3:128680/1‑148 (MQ=255)
                                                                                                                                 ccGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGCCAGCAagag                 >  4:70582/1‑148 (MQ=255)
                                                                                                                                 ccGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACAGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGCCAGCAagag                 >  6:89104/1‑148 (MQ=255)
                                                                                                                                   gAGTACTGACCCGCAGTGAAAACATAGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTCGGCCAGCAAGAGTg               >  7:158982/1‑148 (MQ=255)
                                                                                                                                     gTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGCCAGCAAGAGTggg             <  2:219556/148‑1 (MQ=255)
                                                                                                                                     gTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGCCAGCAAGAGTggg             <  8:134275/148‑1 (MQ=255)
                                                                                                                                        cTGACCCTCAGTGAAATCATTGCTCCGCTGCTGTCTGCACGACCGTAATTTATGGTGCATGCACTGTCGATGGACTCTCTTACCGAATGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGCCAGCAAGAGTGGGGGa          <  8:158982/148‑1 (MQ=255)
                                                                                                                                            cccGTAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGCCAGCAAGAGTGGGGGaaaaa      <  7:17069/148‑1 (MQ=255)
                                                                                                                                            cccGTAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGCCAGCAAGAGTGGGGGaaaaa      <  7:22034/148‑1 (MQ=255)
                                                                                                                                              cGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGCCAGCAag                   >  7:305368/1‑133 (MQ=255)
                                                                                                                                              cGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGCCAGCAag                   <  6:229234/133‑1 (MQ=255)
                                                                                                                                              cGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGCCAGCAag                   >  5:229234/1‑133 (MQ=255)
                                                                                                                                              cGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGGGGTGGTGATTGGTCCCGGTCTGGGCCAGCAag                   <  8:305368/133‑1 (MQ=255)
                                                                                                                                               gCAGAGAAAACATAGCGCCGCTGCTGACTCCACGACCGGAATTGATGGTGCATGAACTGACGATTGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGCCAGCAAGAGTGGGGGAAAAAAGc   <  4:128680/148‑1 (MQ=255)
                                                                                                                                               ccAGTGAAAACATTGCGCCGCTGCTGACTTCACTACCGGAATTGATGGTGCATGAACTGACGATGGAATCTCTTACCGAAAGCCTGGAATTGGCCGATGTTGTTGTGATTGGTCCCGGTCTGGGCCAGCAAGAGTGGGGGAAAAAAGc   <  5:222674/147‑1 (MQ=255)
                                                                                                                                                cAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGATTGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGCCAGCAAGAGTGGGGGAAAAAAGCa  >  5:134165/1‑148 (MQ=255)
                                                                                                                                                cAGTGAAAACATTGCGCCGCTGATGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGCCAGCAAGAGTGGGGGAAAAAAGCa  >  1:75774/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
AAAACCGCGTCGCCCGACTTCGCATAAAGGCGATCACGGGCGGCTGGTAATTATCGGTGGCGATCACGGCACGGCGGGGGCTATTCGTATGACGGGGGAAGCGGCGCTGCGTGCTGGTGCTGGTTTAGTCCGAGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGCCAGCAAGAGTGGGGGAAAAAAGCA  >  NC_000913/4394791‑4395082

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