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,4970GA88.9% 129.3 / 16.3 81S269L (TCG→TTG) gadBglutamate decarboxylase B, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (1/8);  new base A (28/44);  total (29/52)
Fisher's exact test for biased strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCCAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTG  >  CP009273/1565455‑1565642
                                          |                                                                                                                                                 
cccAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTa                                        <  1:1273895/150‑1 (MQ=11)
cccAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTa                                        <  2:4258769/150‑1 (MQ=11)
cccAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTa                                        <  1:1434382/150‑1 (MQ=11)
cccAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTa                                        <  1:1434398/150‑1 (MQ=11)
 ccAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTAt                                       <  2:176844/150‑1 (MQ=14)
 ccAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTAt                                       >  2:4853694/1‑150 (MQ=14)
  cAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATc                                      <  1:1863003/150‑1 (MQ=32)
  cAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATc                                      <  2:2473776/150‑1 (MQ=17)
  cAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATc                                      <  2:2473689/150‑1 (MQ=17)
   aGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCg                                     <  1:533793/150‑1 (MQ=17)
   aGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCg                                     <  1:533922/150‑1 (MQ=17)
    gCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACCCGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCgg                                    <  2:5166240/150‑1 (MQ=17)
    gCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCgg                                    >  2:1691806/1‑150 (MQ=17)
     cGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGc                                   <  1:175185/150‑1 (MQ=17)
     cGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGc                                   >  2:830673/1‑150 (MQ=17)
     cGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGc                                   >  2:3950990/1‑150 (MQ=17)
      ggAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGcc                                  >  1:4749432/1‑150 (MQ=17)
      ggAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGcc                                  >  1:1351221/1‑150 (MQ=17)
      ggAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGcc                                  <  2:1302381/150‑1 (MQ=17)
       gAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCt                                 >  2:3087769/1‑150 (MQ=17)
       gAGCCAGACCGAATTGATGGCCTGCAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCt                                 <  2:1934628/150‑1 (MQ=12)
        aGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTg                                >  1:3622536/1‑150 (MQ=17)
        aGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTg                                <  1:2885586/150‑1 (MQ=17)
        aGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTg                                <  1:2886070/150‑1 (MQ=17)
        aGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTg                                <  1:2885548/150‑1 (MQ=17)
         gCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTgg                               <  2:3034173/150‑1 (MQ=17)
         gCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTgg                               <  2:4118575/150‑1 (MQ=17)
          ccAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGa                              <  2:1972576/150‑1 (MQ=17)
          ccAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGa                              >  2:896813/1‑150 (MQ=17)
           cAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGaa                             >  2:3990830/1‑150 (MQ=17)
           cAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGaa                             >  1:3771895/1‑150 (MQ=17)
           cAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGaa                             >  2:5118062/1‑150 (MQ=17)
           cAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGaa                             >  2:5118066/1‑150 (MQ=17)
            aGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAt                            <  2:1444973/150‑1 (MQ=32)
            aGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAt                            <  2:4044823/150‑1 (MQ=32)
            aGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAt                            <  1:2982611/150‑1 (MQ=17)
            aGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCAAGTCGATGTCGATACCGGTATCGGCCTGGAAt                            >  2:894765/1‑150 (MQ=17)
             gACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAtt                           <  1:1385676/150‑1 (MQ=17)
              aCCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGACttt                          >  2:3250147/1‑150 (MQ=17)
              aCCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAttt                          <  1:95420/150‑1 (MQ=17)
              aCCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAttt                          <  2:1668578/150‑1 (MQ=17)
                cGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTAt                        <  2:63545/150‑1 (MQ=17)
                  caTTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATcc                      <  1:830963/149‑1 (MQ=17)
                  aaTTTATGGCCTGAAGCACTGATCGATTTCACCCGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATcc                      <  2:2752112/150‑1 (MQ=17)
                  aaTTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATcc                      <  2:2752105/150‑1 (MQ=32)
                  aaTTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCGAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATcc                      <  1:3128907/150‑1 (MQ=17)
                  aaTTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATcc                      <  1:3128875/150‑1 (MQ=17)
                   aTTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCa                     <  1:5182977/150‑1 (MQ=17)
                    tttATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAg                    >  1:3034996/1‑150 (MQ=17)
                    tttATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAg                    >  1:1577648/1‑150 (MQ=17)
                      tATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGACTTTATCCAgcg                  >  1:362218/1‑150 (MQ=17)
                      tATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAgcg                  >  1:1268115/1‑150 (MQ=17)
                       aTGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAgcgc                 >  1:141590/1‑150 (MQ=17)
                       aTGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAgcgc                 >  2:3485692/1‑150 (MQ=17)
                         ggCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCAt               <  1:1837443/150‑1 (MQ=32)
                         ggCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCAt               <  2:4097089/150‑1 (MQ=17)
                         ggCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCAt               <  2:1501064/150‑1 (MQ=17)
                         ggCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCAt               <  2:529180/150‑1 (MQ=17)
                         ggCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCAt               <  2:529181/150‑1 (MQ=17)
                          gCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATc              <  2:872347/150‑1 (MQ=17)
                          gCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATc              <  1:4549937/150‑1 (MQ=17)
                          gCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACAAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATc              >  2:1001221/1‑150 (MQ=17)
                           ccTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCg             <  1:2075051/150‑1 (MQ=32)
                             tGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTg           <  2:4617021/150‑1 (MQ=17)
                             tGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTg           <  2:2558752/150‑1 (MQ=17)
                              gAAGCACTGATCAATTTCACCCGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGc          <  1:4743858/150‑1 (MQ=17)
                              gAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGc          >  1:2671648/1‑150 (MQ=17)
                              gAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGc          <  1:4045494/150‑1 (MQ=17)
                               aaGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCa         <  2:462815/150‑1 (MQ=17)
                                 gCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGc       >  1:993359/1‑150 (MQ=32)
                                 gCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCCGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAGTTATCCAGCGCATCGTGCAGc       >  1:2436425/1‑150 (MQ=12)
                                    cTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGt    >  1:1755332/1‑150 (MQ=17)
                                    cTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGt    >  1:1012471/1‑150 (MQ=17)
                                    cTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGt    >  1:4722626/1‑150 (MQ=17)
                                     tGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGtt   <  2:1310718/150‑1 (MQ=32)
                                     tGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGtt   <  2:5054894/150‑1 (MQ=17)
                                     tGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGtt   <  2:464459/150‑1 (MQ=17)
                                     tGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGtt   <  2:2210588/150‑1 (MQ=17)
                                      gATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTtg  <  2:551905/150‑1 (MQ=17)
                                      gATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTtg  <  2:552009/150‑1 (MQ=17)
                                      gATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTtg  >  2:211127/1‑150 (MQ=17)
                                          |                                                                                                                                                 
CCCAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTG  >  CP009273/1565455‑1565642

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: