Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,052,254 C→G 5.9% R87P (CGA→CCA)  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,2540CG 5.9% 215.4 / 3.7 68R87P (CGA→CCA) PP_4452NAD/NADP octopine/nopaline dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base C (33/30);  new base G (2/2);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.78e-01

GCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTCGAACACACCCGTTGCATGTAACTTGAAGCGCCCCAACGAAGGCTGATCCGA  >  NC_002947/5052107‑5052399
                                                                                                                                                   |                                                                                                                                                 
gCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCg                                                                                                                                                  >  4:354727/1‑149 (MQ=255)
   gTTCTGCCAAAACTGGCGACCACAAGATCTAACCCCAGCTTCCTCGCGCGCTCCCGAAGATAAAGAGAAGACAACGAGCCAATTGAACTAACGATAAACGTCTGGCCGCTTCTCAAGCTTGTGATCGCAGCATCAATTCCCTGTCGGTg                                                                                                                                               >  8:197721/1‑149 (MQ=255)
     tgtgCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTgc                                                                                                                                              <  2:64191/148‑1 (MQ=255)
       tgCCACAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCa                                                                                                                                           <  6:402806/149‑1 (MQ=255)
        gCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCa                                                                                                                                           <  8:9584/148‑1 (MQ=255)
            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                                                                                                               >  6:406256/1‑147 (MQ=255)
                                     ccAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGt                                                                                                               >  3:182222/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                                                                                               <  1:163555/139‑1 (MQ=255)
                                                             gATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCgt                                                                                               >  2:168973/1‑139 (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)
                                                                         aCAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTgc                                                                                                                                              <  1:103491/80‑1 (MQ=255)
                                                                         aCAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTgc                                                                                                                                              >  2:103491/1‑80 (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                                                                              <  1:276348/102‑1 (MQ=255)
                                                                                                                   tcAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGa                                                                              >  2:276348/1‑102 (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                      >  2:309245/2‑135 (MQ=255)
                                                                                                                                          cttGCCTGTGAATGCGCACCGACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTCGAACACACCCGTTGCATGTAACTTGAAgcgc                      <  1:309245/134‑1 (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                             >  8:316789/1‑122 (MQ=255)
                                                                                                                                                tGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTCGAACACACCCGTTGCATGTAACtt                             <  7:316789/122‑1 (MQ=255)
                                                                                                                                                tGTCGGTGCGCATTCACAGGCAAGCAGATAATTACAACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTCGAACACACCCGTTGCATGTAACTTGAAGCGCCCCAACGAAGGCTGATCCGa  >  1:60075/1‑149 (MQ=255)
                                                                                                                                                  tCGGTTCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTGTCAATGCTCACACTTCCCTCGAACACACCCGTTGCATGTAACTTGAAGCGCCCCAACGAAGGCTGATCCGa  <  8:119108/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGACAACGAGCCCATTGAACTTACGATAACCGTCTGGCCGCTTCTCAAGCTTGGGAGCGCAGCATCAATTGCCTGTCGGTGCGCATTCACAGGCAAGCAGATAATTACCACGTCAGCCGCTGAAACCGTGTCCTCAGCGTTATTGACTCTATCAATGCTCACACTTCCCTCGAACACACCCGTTGCATGTAACTTGAAGCGCCCCAACGAAGGCTGATCCGA  >  NC_002947/5052107‑5052399

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: ‑