Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,565,497 G→A 100% S269L (TCG→TTG)  gadB ← glutamate decarboxylase B, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,565,4970GA98.4% 135.6 / ‑4.0 64S269L (TCG→TTG) gadBglutamate decarboxylase B, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (28/35);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCCAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGT  >  CP009273/1565455‑1565643
                                          |                                                                                                                                                  
cccAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTa                                         <  1:160628/150‑1 (MQ=11)
cccAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTa                                         <  1:161028/150‑1 (MQ=11)
  cAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATc                                       >  2:1694446/1‑150 (MQ=17)
  cAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATc                                       >  2:1694403/1‑150 (MQ=17)
  cAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATc                                       <  1:30840/150‑1 (MQ=17)
  cAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATc                                       <  1:30916/150‑1 (MQ=17)
   cgCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCg                                      <  1:3227787/149‑1 (MQ=17)
   aGCGGAGCCAGACCGAATTTCTGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCg                                      <  1:953727/150‑1 (MQ=17)
   aGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCg                                      <  1:1711464/150‑1 (MQ=17)
   aGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCg                                      <  1:3238150/150‑1 (MQ=17)
    gCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACCCGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCgg                                     <  2:3556252/150‑1 (MQ=17)
      ggAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGcc                                   >  2:3645033/1‑150 (MQ=17)
      ggAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGcc                                   >  2:3642990/1‑150 (MQ=17)
      ggAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGcc                                   <  1:4081541/150‑1 (MQ=17)
          ccAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGa                               <  1:1737742/150‑1 (MQ=17)
           cAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGac                              >  2:1492943/1‑149 (MQ=17)
             gACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAtt                            >  2:3920168/1‑150 (MQ=17)
              aCCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAttt                           <  1:2943990/150‑1 (MQ=17)
               ccGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCCGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTa                          >  2:2820027/1‑150 (MQ=17)
               ccGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTa                          >  2:462209/1‑150 (MQ=17)
               ccGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTa                          >  2:2820021/1‑150 (MQ=17)
                cGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTAt                         >  2:3739053/1‑150 (MQ=17)
                cGAATTTAGGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTAt                         <  2:888363/150‑1 (MQ=17)
                  aaTTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATcc                       <  1:4031791/150‑1 (MQ=17)
                  aaTTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATcc                       <  1:4031800/150‑1 (MQ=17)
                  aaTTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATcc                       <  1:4031828/150‑1 (MQ=17)
                  aaTTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATcc                       >  2:3588715/1‑150 (MQ=17)
                  aaTTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATcc                       >  2:3587718/1‑150 (MQ=17)
                   aTTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGACTTTATCCa                      >  1:3395261/1‑150 (MQ=17)
                   aTTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCa                      <  2:991775/150‑1 (MQ=17)
                   aTTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCa                      >  2:4140222/1‑150 (MQ=17)
                   aTTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCa                      >  1:4536074/1‑150 (MQ=17)
                    tttATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAg                     <  1:864215/150‑1 (MQ=17)
                     ttATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGACCGGGGCCAGGACGCCACCGCTGGCAGCGTCGATGTGCAGGTCGATGTCGATACCGGTATCGGCCTGGACTTTATCCAgc                    >  2:2971236/1‑150 (MQ=11)
                     ttATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAgc                    >  2:2971261/1‑150 (MQ=17)
                      tATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAgcg                   <  2:2214705/150‑1 (MQ=17)
                       aTGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAgcgc                  <  1:2768096/150‑1 (MQ=17)
                       aTGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAgcgc                  <  1:2265932/150‑1 (MQ=17)
                       aTGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAgcgc                  <  1:2265907/150‑1 (MQ=17)
                         ggCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCAt                >  1:1316700/1‑150 (MQ=17)
                         ggCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCAt                >  1:1624903/1‑150 (MQ=17)
                         ggCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCAt                >  1:2662960/1‑150 (MQ=17)
                          gCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGACTTTATCCAGCGCATc               >  2:936785/1‑150 (MQ=17)
                          gCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATc               <  2:783562/150‑1 (MQ=17)
                          gCCTGAAGCACTGATCAATGTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATc               <  2:783536/150‑1 (MQ=17)
                           ccTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACAATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCg              <  2:4184979/150‑1 (MQ=17)
                           ccTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACAATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCg              <  2:4184964/150‑1 (MQ=17)
                            cTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGt             >  2:4085661/1‑150 (MQ=17)
                              gAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGc           <  2:2681531/150‑1 (MQ=17)
                              gAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGc           <  2:415126/150‑1 (MQ=17)
                               aaGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCa          <  2:1081203/150‑1 (MQ=32)
                               aaGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCa          <  1:3421392/150‑1 (MQ=17)
                               aaGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCa          <  1:4244225/150‑1 (MQ=17)
                               aaGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCa          <  2:3797643/150‑1 (MQ=17)
                                    cTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGt     >  1:2252485/1‑150 (MQ=17)
                                     tGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGtt    >  2:4246513/1‑150 (MQ=17)
                                     tGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGtt    >  2:2923677/1‑150 (MQ=17)
                                      gATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTtg   >  1:590675/1‑150 (MQ=17)
                                      gATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTtg   <  2:653166/150‑1 (MQ=17)
                                       aTCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCggtggt  <  2:1345087/150‑1 (MQ=17)
                                       aTCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTtgt  >  2:1783438/1‑150 (MQ=17)
                                       aTCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTtgt  >  1:2073334/1‑150 (MQ=17)
                                       aTCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTtgt  <  2:1345047/150‑1 (MQ=17)
                                       aTCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTtgt  <  2:1345303/150‑1 (MQ=17)
                                          |                                                                                                                                                  
CCCAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGT  >  CP009273/1565455‑1565643

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