Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,798,255 C→T 6.6% V215V (GTG→GTA PP_3358 ← hydroxycinnamoyl‑CoA hydratase‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,798,2550CT 6.6% 178.0 / 2.1 61V215V (GTG→GTAPP_3358hydroxycinnamoyl‑CoA hydratase‑lyase
Reads supporting (aligned to +/- strand):  ref base C (30/26);  new base T (2/2);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.62e-01

GAGGAACTGCTTCATGCCCTGCTCGCGGCCGCCGGTAGTGTCCAGCAGGCGCGACTGGTCGAGCTTGGCGTAGAGGTAGTCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACCAGGCCCATCTCGGCAGCCTTGCGACCATCGAAGGTCTTGCCGGTCATGATGTAGTACAG  >  NC_002947/3798114‑3798398
                                                                                                                                              |                                                                                                                                               
gAGGAACTGCTTCATGCCCTGCTCGCGGCCGCCGGTAGTGTCCAGCAGGCGCGACTGGTCGAGCTTGGCGTAGAGGTAGTCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCaccacc                                                                                                                                           <  3:134080/148‑1 (MQ=255)
gAGGAACTGCTTCATGCCCTGCTCGCGGCCGCCGGTAGTGTCCAGCAGGCGCGACTGGTCGAGCTTGGCGTAGAGGTAGTCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCg                                                                                                                                          >  2:233636/1‑149 (MQ=255)
     aCTGCTTCATGCCCTGCTCGCGGCCGCCGGTAGTGTCCAGCAGGCGCGACTGGTCGAGCTTGGCGTAGAGGTAGTCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACgcagc                                                                                                                                                 <  4:127126/137‑1 (MQ=255)
     aCTGCTTCATGCCCTGCTCGCGGCCGCCGGTAGTGTCCAGCAGGCGCGACTGGTCGAGCTTGGCGTAGAGGTAGTCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACgcagc                                                                                                                                                 >  3:127126/1‑137 (MQ=255)
      cTGCTTCATGCCCTGCTCGCGGCCGCCGGTAGTGTCCAGCAGGCGCGACTGGTCGAGCTTGGCGTAGAGGTAGTCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCt                                                                                                                                    >  7:235353/1‑149 (MQ=255)
            cctGCCCTGCTCGCGGCCGCCGGTAGTGTCCAGCAGGCGCGACTGGTCGAGCTTGGCGTAGAGGTAGTCCTCGTTCTGTTCCCAGGTCCGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCca                                                                                                                              <  7:122689/147‑1 (MQ=255)
               gCCCTGCTCGCGGCCGCCGGTAGTGTCCAGCAGGCGCGACTGGTCGAGCTTGGCGTAGAGGTAGTCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCcag                                                                                                                             >  4:6967/1‑147 (MQ=255)
                cccTGCTCGCGGCCGCCGGTAGTGTCCAGCAGGCGCGACTGGTCGAGCTTGGCGTAGAGGTAGTCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCcagca                                                                                                                           <  1:233636/149‑1 (MQ=255)
                      tCGCGGCCGCCGGTAGTGTCCAGCAGGCGCGACTGGTCGAGCTTGGCGTAGAGGTAGTCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCca                                                                                                                              <  8:87613/139‑1 (MQ=255)
                      tCGCGGCCGCCGGTAGTGTCCAGCAGGCGCGACTGGTCGAGCTTGGCGTAGAGGTAGTCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCca                                                                                                                              >  7:87613/1‑139 (MQ=255)
                              gccgGTAGTGTCCAGCAGGCGCGACTGGTCGAGCTTGGCGTAGAGGTAGTCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACt                                                                                                            >  6:144205/1‑149 (MQ=255)
                                  gTAGTGTCCAGCAGGCGCGACTGGTCGAGCTTGGCGTAGAGGTAGTCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTcgcg                                                                                                        >  6:140636/1‑149 (MQ=255)
                                      tgtCCAGCAGGCGCGACTGGTCGAGCTTGGCGTAGAGGTAGTCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGggtg                                                                                                    >  2:32615/1‑149 (MQ=255)
                                         ccAGCACGCGCGACTGGTCGAGCTTGACGTAGAGGTAGTCCTCGTTCACTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGACATTCTTCTCGGAACGCAGCACCACAGCGTTCTTTTCCAGCAGGTTCAGCGCCAACTCCCGGGTGttt                                                                                                 >  6:285471/1‑149 (MQ=255)
                                          cagcagGCGCGACTGGTCGATCTTGGCGTAGAGGTAGTCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCCCGCAGCACCACCGGGTTCTTTTCCAGCAGGTTAAGCGCCCACTCGAGGGTGGttt                                                                                                >  6:216331/1‑149 (MQ=255)
                                              agGCGCGACTGGTCGAGCTTGGCGTAGAGGTAGTCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACg                                                                                            <  1:32615/149‑1 (MQ=255)
                                                gcgcGACTGGTCGAGCTTGGCGTAGAGGTAGTCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGTACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCAc                                                                                             <  3:76316/146‑1 (MQ=255)
                                                gcgcGACTGGTCGAGCTTGGCGTAGAGGTAGTCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGTACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCAc                                                                                             <  3:80570/146‑1 (MQ=255)
                                                gcgcGACTGGTCGAGCTTGGCGTAGAGGTAGTCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGTACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCAc                                                                                             >  4:80570/1‑146 (MQ=255)
                                                gcgcGACTGGTCGAGCTTGGCGTAGAGGTAGTCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGTACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCAc                                                                                             >  4:76316/1‑146 (MQ=255)
                                                  gcgACTGGTCGAGCTTGGCGTAGAGGTAGTCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGc                                                                                        >  6:180010/1‑149 (MQ=255)
                                                     aCTGGTCGAGCTTGGCGTAGAGGTAGTCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCg                                                                                     <  8:235353/149‑1 (MQ=255)
                                                             aGCTTGGCGTAGAGGTAGTCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCgg                                                                             <  3:6967/149‑1 (MQ=255)
                                                                       agagGTAGTCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCg                                                                   >  1:243666/1‑149 (MQ=255)
                                                                       agagGTAGTCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCg                                                                   >  3:43414/1‑149 (MQ=255)
                                                                         agGTTGTCCTCTTTCTGTTCCGAGGTCAGATGGCGGCAACGGTTGATGCCATTGTTCGCGTCACGCTTGACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCATCTCGCTTGTGGTTTCACTCAGCTCGGCCAGCGGCACACTGTCGt                                                                  <  5:285471/148‑1 (MQ=255)
                                                                           gTAGTCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCa                                                               <  5:144205/149‑1 (MQ=255)
                                                                           gTAGTCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCAGGTGGTTTCACGCAGCTCg                                                                                     >  1:289021/1‑127 (MQ=255)
                                                                           gTAGTCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCAGGTGGTTTCACGCAGCTCg                                                                                     <  2:289021/127‑1 (MQ=255)
                                                                              ggccTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACCa                                                            <  8:150172/147‑1 (MQ=255)
                                                                               tCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACCAg                                                           >  6:285336/1‑149 (MQ=255)
                                                                               tCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACCAg                                                           <  5:84884/149‑1 (MQ=255)
                                                                                ccTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCcag                                                                                                                             >  4:290027/1‑82 (MQ=255)
                                                                                ccTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCcag                                                                                                                             <  3:290027/82‑1 (MQ=255)
                                                                                ccTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTATCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACCAgg                                                          >  7:223267/1‑149 (MQ=255)
                                                                                 cTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACCAGGc                                                         >  1:75502/1‑149 (MQ=255)
                                                                                 cTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACCAGGc                                                         <  5:285336/149‑1 (MQ=255)
                                                                                     ttCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACAACAGGGTTCTTTTCAAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACAAGGCCCa                                                      >  3:186933/1‑148 (MQ=255)
                                                                                      tCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACACGGCCCAtc                                                    >  6:222314/1‑149 (MQ=255)
                                                                                          ttCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCTCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACCAGGCCCATCTCgg                                                >  6:154691/1‑149 (MQ=255)
                                                                                             ccAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTcc                                                                                                                               <  4:87809/67‑1 (MQ=255)
                                                                                             ccAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTcc                                                                                                                               >  3:87809/1‑67 (MQ=255)
                                                                                             ccAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACCAGGCCCATCTCGGGAg                                             >  4:158666/1‑149 (MQ=255)
                                                                                              cAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTGCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACCAGGCCCATCTCGGCAGc                                            <  4:43999/149‑1 (MQ=255)
                                                                                                         gcgGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACCAGGCCCATCTCGGCAGCCTTGCGACCAt                                 <  3:158666/149‑1 (MQ=255)
                                                                                                              aaCGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACCAGGCCCATCTCGGCAGCCTTGCGACCATCGAAg                            <  5:140636/149‑1 (MQ=255)
                                                                                                                cGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACCAGGCCCATCTCGGCAGCCTTGCGACCATCGAAGGt                          <  5:154691/149‑1 (MQ=255)
                                                                                                                 gCTTGAAGCCAATCTTCGCGGCACGCAGCACCACCTGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACCAGGCCCATCTCGGCAGCCTTGCGACCATCGAAGGTc                         <  4:186933/149‑1 (MQ=255)
                                                                                                                    ttGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACCAGGCCCATCTCGGCAGCCTTGCGACCATCGAAg                            <  3:302405/144‑1 (MQ=255)
                                                                                                                    ttGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACCAGGCCCATCTCGGCAGCCTTGCGACCATCGAAg                            >  4:302405/1‑144 (MQ=255)
                                                                                                                            aTTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTGCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACCAGGCCCATCTCGGCAGCCTTGCGACCATCGAAGGTCTTGCCGGt                 >  8:181819/1‑147 (MQ=255)
                                                                                                                                ttcGCGGCACGCAGCCCCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTGCACCCGGCCCCTCTCGGCCGCCTTGCGACCAGCGCAGGTCTTg                      >  2:304408/1‑138 (MQ=255)
                                                                                                                                ttcGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTc                                                                    >  8:281923/1‑92 (MQ=255)
                                                                                                                                ttcGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTc                                                                    <  7:281923/92‑1 (MQ=255)
                                                                                                                                ttcGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACCAGGCCCATCTCGGCAGCCTTGCGACCATCGAAGGTCTTg                      <  1:309291/138‑1 (MQ=255)
                                                                                                                                ttcGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACCAGGCCCATCTCGGCAGCCTTGCGACCATCGAAGGTCTTg                      >  2:309291/1‑138 (MQ=255)
                                                                                                                                ttcGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACCAGGCCCATCTCGGCAGCCTTGCGACCATCGAAGGTCTTg                      <  1:304408/138‑1 (MQ=255)
                                                                                                                                ttcGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACCAGGCCCATCTCGGCAGCCTTGCGACCATCGAAGGTCTTGCCGGTCatgat           >  7:195277/1‑149 (MQ=255)
                                                                                                                                  cgcgGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACCAGGCCCATCTCGGCAGCCTTGCGACCATCGAAGGTCTTGCCGGTCatgatg          >  2:218116/1‑148 (MQ=255)
                                                                                                                                   gcgGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACCAGGCCCATCTCGGCAGCCTTGCGACCATCGAAGGTCTTGCCGGTCATGATgta        <  2:75502/149‑1 (MQ=255)
                                                                                                                                     ggCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACCAGGCCCATCTCGGCAGCCTTGCGACCATCGAAGGTCTTGCCGGTCATGATgtagt      <  5:222314/149‑1 (MQ=255)
                                                                                                                                     ggCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACCAGGCCCATCTCGGCAGCCTTGCGACCATCGAAGGTCTTGCCGGTCATGATgtag       >  6:139175/1‑148 (MQ=255)
                                                                                                                                         cGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACCAGGCCCATCTCGGCAGCCTTGCGACCATCGAAGGTCTTGCCGGTCATGATGTAGTACAg  <  8:195277/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                               
GAGGAACTGCTTCATGCCCTGCTCGCGGCCGCCGGTAGTGTCCAGCAGGCGCGACTGGTCGAGCTTGGCGTAGAGGTAGTCCTCGTTCTGTTCCCAGGTCAGCTCGCGGCAACGCTTGAAGCCATTCTTCGCGGCACGCAGCACCACCGGGTTCTTTTCCAGCAGGTTCAGCGCCAACTCGCGGGTGGTTTCACGCAGCTCGGCCAGCGGCACACTGTCGTTCACCAGGCCCATCTCGGCAGCCTTGCGACCATCGAAGGTCTTGCCGGTCATGATGTAGTACAG  >  NC_002947/3798114‑3798398

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: