Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,052,265 C→G 6.7% V83V (GTG→GTC)  PP_4452 ← NAD/NADP octopine/nopaline dehydrogenase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,052,2650CG 6.7% 183.0 / 4.0 60V83V (GTG→GTC) PP_4452NAD/NADP octopine/nopaline dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base C (29/26);  new base G (2/2);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01

AAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTCGAACACACCCGTTGCATGTAACTTGAAGCGCCCCAACGAAGGCTGATCCGATGAGCGCCCATTC  >  NC_002947/5052119‑5052412
                                                                                                                                                  |                                                                                                                                                   
aaaCTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTcac                                                                                                                                                   <  5:307652/149‑1 (MQ=255)
           gacAAGATCTACCCCCAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTGAATGCGCACCGACAGGCaattgat                                                                                                                                        <  8:217920/149‑6 (MQ=255)
           gacAAGATCTACCCCCAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAg                                                                                                                                            >  6:337996/1‑145 (MQ=255)
                      cccccAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTAc                                                                                                                                 <  3:354727/145‑1 (MQ=255)
                      cccccAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCAc                                                                                                                              <  8:36474/148‑1 (MQ=255)
                      cccccAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCAc                                                                                                                              <  3:149126/148‑1 (MQ=255)
                        cccAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCGATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGt                                                                                                                            >  8:55479/1‑148 (MQ=255)
                         ccAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGt                                                                                                                            >  3:182222/1‑147 (MQ=255)
                         ccAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGt                                                                                                                            >  6:406256/1‑147 (MQ=255)
                          cAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAg                                                                                                                         >  4:397802/1‑149 (MQ=255)
                                     cgcgcTCCCGAAGATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAAcc                                                                                                              >  8:213096/1‑149 (MQ=255)
                                        gcTCCCGAAGATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCgtg                                                                                                           <  8:231541/149‑1 (MQ=255)
                                         cTCCCGAAGATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCgtgt                                                                                                          <  4:251127/149‑1 (MQ=255)
                                            ccGAAGATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTcc                                                                                                        >  6:258421/1‑148 (MQ=255)
                                             cGAAGATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCt                                                                                                       >  7:83760/1‑148 (MQ=255)
                                             cGAAGATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAAACGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCt                                                                                                       >  8:82252/1‑148 (MQ=255)
                                                 gATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAAt                                                                                                                                    >  2:163555/1‑115 (MQ=255)
                                                 gATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCgt                                                                                                            >  2:168973/1‑139 (MQ=255)
                                                 gATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCgt                                                                                                            <  1:163555/139‑1 (MQ=255)
                                                 gATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCgt                                                                                                            <  1:168973/139‑1 (MQ=255)
                                                      aGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCGATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTAtt                                                                                             <  7:55479/149‑1 (MQ=255)
                                                         gcCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGAc                                                                                          <  8:83760/149‑1 (MQ=255)
                                                            gACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGActc                                                                                        >  5:32008/1‑148 (MQ=255)
                                                             aCAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTa                                                                                      >  8:34060/1‑149 (MQ=255)
                                                                      ccATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCt                                                                              <  5:406256/148‑1 (MQ=255)
                                                                            aaCTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGAGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACtt                                                                       >  5:85900/1‑149 (MQ=255)
                                                                               ttACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTGTCAATGCTCACACTTccc                                                                    >  7:119108/1‑149 (MQ=255)
                                                                               ttACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTca                                                                            <  6:301632/141‑1 (MQ=255)
                                                                               ttACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTccc                                                                    <  5:258421/149‑1 (MQ=255)
                                                                               ttACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGAGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTca                                                                            >  5:301632/1‑141 (MQ=255)
                                                                                  cGATACCCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCACTTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTc                                                                  <  8:188017/148‑1 (MQ=255)
                                                                                  cGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTc                                                                  <  3:397802/148‑1 (MQ=255)
                                                                                  cGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTCg                                                                 <  7:34060/149‑1 (MQ=255)
                                                                                      aaCCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCAAAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTCGAaca                                                             >  5:291189/1‑149 (MQ=255)
                                                                                       aCCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTCGAacac                                                            <  5:337996/149‑1 (MQ=255)
                                                                                        ccGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCATCAGCGTTATTGACTCTATCAATGCTCACACTTCCCACGAacaca                                                           >  3:379959/1‑149 (MQ=255)
                                                                                         cttCTGGCCGCTTCTCAAGCTTGGTAGCGCAGCATGAATTGCCTGTCGTTGCTCATTCACAGGCAAGCAGATAATTACCACGTCAGCAGCTGAAACCGTGTCCTCAGCGTTATTGAATCTATCAATGCTCACACTTCCCTCGAacacac                                                          <  7:197721/147‑1 (MQ=255)
                                                                                                cagCTTCTCAAGCTGGGGAGCGCAGCATCAATTGCCGGTCGGTGCCCATTCACAGCCAAGCAGATAATTACCACGTCAGGCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTCGAACACACCCGTTGc                                                   <  8:360860/147‑1 (MQ=255)
                                                                                                  gCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTCGAACACACCCGTTGCAt                                                 <  7:82252/149‑1 (MQ=255)
                                                                                                       tcAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGa                                                                                           >  2:276348/1‑102 (MQ=255)
                                                                                                       tcAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGa                                                                                           <  1:276348/102‑1 (MQ=255)
                                                                                                        cAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGa                                                                                           >  5:163000/1‑101 (MQ=255)
                                                                                                        cAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGa                                                                                           <  6:163000/101‑1 (MQ=255)
                                                                                                         aaGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTCGaa                                                               <  5:105414/128‑1 (MQ=255)
                                                                                                         aaGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGAGTTATTGACTCTATCAATGCTCACACTTCCCTCGaa                                                               >  6:105414/1‑128 (MQ=255)
                                                                                                                ggAGCGCAGCATCAATTGCCTGTCGGTGCGCATTAACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTCGAACACACCCGTTGCATGTAACTTGAAgcgc                                   >  6:182255/1‑149 (MQ=255)
                                                                                                                           tCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTCGAACACACCCGTTGCATGTAACTTGAAGCGCCCCAACGAAgg                        >  5:149577/1‑149 (MQ=255)
                                                                                                                            cAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTCGAACACACCCGTTGCAt                                                 >  6:357656/1‑123 (MQ=255)
                                                                                                                            cAATTGCCTGTCGGTCCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTCGAACACACCCGTTGCAt                                                 <  5:357656/123‑1 (MQ=255)
                                                                                                                             aaTTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTCGAACACACCCGTTGCATGTAACTTGAAGCGCCCCAACGAAGGCt                      >  8:142261/1‑149 (MQ=255)
                                                                                                                             aaTTGCCTGTAGGTTCTCATTCACAGGAAAGGAGATAATTACCAGGTCAGCCGCTGAAACCGTGTCCTCAGCTTTATTGACTCTATCAATGCTCACACTTCCCTAGAACACACCCGTTGAATGTAACTTGAAGCGCCCCAACGAAGGc                       <  4:379959/148‑1 (MQ=255)
                                                                                                                              cttGCCTGTGAATGCGCACCGACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTCGAACACACCCGTTGCATGTAACTTGAAgcgc                                   <  1:309245/134‑1 (MQ=255)
                                                                                                                              cttGCCTGTGAATGCGCACCGACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTCGAACACACCCGTTGCATGTAACTTGAAgcgc                                   >  2:309245/2‑135 (MQ=255)
                                                                                                                              aTTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTCGAACACACCCGTTGCATGTAACTTGAAGCGCCCCAACGAAGGCt                      >  4:371931/1‑148 (MQ=255)
                                                                                                                               ttGCCTGTGAATGCGCACCGACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTCGAACACACCCGTTGCATGTAACTTGAAGCGCCCCAACGAAGGCTGa                    >  8:263389/1‑149 (MQ=255)
                                                                                                                               ttGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTCGAACACACCCGTTGCATGTAACTTGAAGCGCCCCAACGAAGGCTg                     >  8:166069/1‑148 (MQ=255)
                                                                                                                                  ccTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTCGAACACACCCGTTGCATGTAACTTGAAGCGCCCCAACGAAGGCTGATcc                 >  3:273270/1‑149 (MQ=255)
                                                                                                                                    tGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTCGAACACACCCGTTGCATGTAACtt                                          <  7:316789/122‑1 (MQ=255)
                                                                                                                                    tGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTCGAACACACCCGTTGCATGTAACtt                                          >  8:316789/1‑122 (MQ=255)
                                                                                                                                    tGTCGGTGCGCATTCACAGGCAAGCAGATAATTACAACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTCGAACACACCCGTTGCATGTAACTTGAAGCGCCCCAACGAAGGCTGATCCGa               >  1:60075/1‑149 (MQ=255)
                                                                                                                                      tCGGTTCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTGTCAATGCTCACACTTCCCTCGAACACACCCGTTGCATGTAACTTGAAGCGCCCCAACGAAGGCTGATCCGa               <  8:119108/147‑1 (MQ=255)
                                                                                                                                                 acacaGGCAAGCAGATAATTACAACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTATATCAATGCTCACACTTCCCTCGAACACACCCGTAGCATGTAACTTGAAGCGCCCCAACGAAGGCTGATCCGATGACCGCACAtt   >  8:113841/2‑148 (MQ=255)
                                                                                                                                                  cacaGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTCGAACACACCCGTTGCATGTAACTTGAAGCGCCCCAACGAAGGCTGATCCGATGAGCGCCCATTc  <  6:32008/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTCGAACACACCCGTTGCATGTAACTTGAAGCGCCCCAACGAAGGCTGATCCGATGAGCGCCCATTC  >  NC_002947/5052119‑5052412

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