Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,565,497 G→A 90.0% 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,4970GA90.0% 189.0 / 24.6 110S269L (TCG→TTG) gadBglutamate decarboxylase B, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (2/9);  new base A (51/48);  total (53/57)
Fisher's exact test for biased strand distribution p-value = 5.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 110 total aligned reads displayed.

CCAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTG  >  CP009273/1565456‑1565642
                                         |                                                                                                                                                 
ccAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTa                                        >  2:3682105/1‑149 (MQ=11)
ccAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTa                                        >  2:2988801/1‑149 (MQ=11)
 cAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTAt                                       <  1:3527967/149‑1 (MQ=31)
 cAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTAt                                       <  2:3325361/149‑1 (MQ=14)
  aGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATc                                      >  2:4649748/1‑149 (MQ=17)
   gCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCg                                     >  2:1587826/1‑149 (MQ=17)
   gCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCg                                     >  2:3922666/1‑149 (MQ=17)
     ggAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGc                                   >  2:1120602/1‑149 (MQ=17)
     ggAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGc                                   >  2:1120521/1‑149 (MQ=17)
     ggAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGc                                   >  2:1717715/1‑149 (MQ=17)
      gAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGcc                                  <  2:3297269/149‑1 (MQ=32)
      gAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGcc                                  <  2:3297287/149‑1 (MQ=32)
      gAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGcc                                  <  2:868667/149‑1 (MQ=32)
      gAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGcc                                  <  2:868942/149‑1 (MQ=32)
      gAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGcc                                  <  1:4404785/149‑1 (MQ=17)
      gAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGcc                                  <  1:2787270/149‑1 (MQ=17)
        gCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTg                                >  2:381387/1‑149 (MQ=32)
         ccAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTgg                               <  2:619669/149‑1 (MQ=17)
         ccAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTgg                               <  2:3894617/149‑1 (MQ=17)
         ccAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTgg                               <  2:619347/149‑1 (MQ=17)
         ccAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTgg                               <  2:619654/149‑1 (MQ=17)
         ccAGACCGAATTTAGGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTgg                               <  2:619614/149‑1 (MQ=17)
          cAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGa                              <  1:609944/149‑1 (MQ=32)
          cAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGa                              <  1:3103067/149‑1 (MQ=17)
          cAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGa                              <  1:1749461/149‑1 (MQ=17)
           aGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGaa                             >  2:2299249/1‑149 (MQ=17)
           aGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGaa                             <  2:1316320/149‑1 (MQ=17)
           aGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGaa                             <  2:4553384/149‑1 (MQ=17)
           aGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGaa                             >  2:292281/1‑149 (MQ=17)
            gACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAt                            >  1:3815069/1‑149 (MQ=17)
            gACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAt                            >  2:1327174/1‑149 (MQ=17)
            gACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAt                            >  2:1327217/1‑149 (MQ=17)
            gACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAt                            >  2:2739530/1‑149 (MQ=17)
            gACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAt                            >  1:54493/1‑149 (MQ=17)
             aCCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAtt                           >  1:3875384/1‑149 (MQ=17)
             aCCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAtt                           <  2:1750596/149‑1 (MQ=17)
             aCCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAtt                           >  1:3463681/1‑149 (MQ=17)
              ccGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAttt                          <  1:2316405/149‑1 (MQ=32)
              ccGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCTGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAttt                          >  2:1010524/1‑149 (MQ=17)
              ccGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAttt                          <  2:20285/149‑1 (MQ=17)
              ccGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAttt                          >  2:1010538/1‑149 (MQ=17)
               cGAATTTATGGCCTGAAGCACTGATCACTTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTa                         <  2:4290372/149‑1 (MQ=17)
               cGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTa                         <  2:4290472/149‑1 (MQ=17)
               cGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTa                         <  2:4290378/149‑1 (MQ=17)
               cGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTa                         >  2:735159/1‑149 (MQ=17)
               cGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTa                         <  1:4203394/149‑1 (MQ=17)
               cGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTa                         <  2:4290374/149‑1 (MQ=17)
               cGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTa                         >  1:3475029/1‑149 (MQ=17)
                gAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTAt                        >  2:3629685/1‑149 (MQ=17)
                gAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTAt                        >  1:1947339/1‑149 (MQ=17)
                gAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCCGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTAt                        <  2:3054317/149‑1 (MQ=17)
                gAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCCGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTAt                        <  2:3054171/149‑1 (MQ=17)
                gAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCCGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTAt                        <  2:2941595/149‑1 (MQ=17)
                 aaTTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATc                       >  1:89106/1‑149 (MQ=32)
                 aaTTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATc                       >  2:2027643/1‑149 (MQ=17)
                 aaTTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATc                       >  2:2027688/1‑149 (MQ=17)
                  aTTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATcc                      <  2:1099356/149‑1 (MQ=17)
                  aTTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATcc                      >  2:4097675/1‑149 (MQ=17)
                  aTTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATcc                      <  2:1795674/149‑1 (MQ=17)
                   tttCTGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCa                     <  2:2188174/149‑1 (MQ=17)
                   tttATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCa                     <  1:808946/149‑1 (MQ=17)
                   tttATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCa                     >  1:2997443/1‑149 (MQ=17)
                     tATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAgc                   >  1:2693453/1‑149 (MQ=17)
                     tATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAgc                   <  2:1116156/149‑1 (MQ=17)
                     tATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAgc                   >  1:351732/1‑149 (MQ=17)
                      aTGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAgcg                  >  1:3906148/1‑149 (MQ=17)
                      aTGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAgcg                  >  1:1514464/1‑149 (MQ=17)
                         gCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCAt               >  2:2853671/1‑149 (MQ=17)
                         gCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCAt               >  1:3973245/1‑149 (MQ=17)
                         gCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCAt               >  2:180/1‑149 (MQ=17)
                           cTGAAGCTCTGATCAATTTCACCCGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGGAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCg             <  2:4222462/149‑1 (MQ=255)
                             gAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTg           >  1:2326383/1‑149 (MQ=17)
                             gAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTg           >  1:2326009/1‑149 (MQ=17)
                             gAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTg           >  1:1420232/1‑149 (MQ=17)
                              aaGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGc          <  2:1962344/149‑1 (MQ=32)
                              aaGCACTGATCGATGTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGc          <  2:4297382/149‑1 (MQ=17)
                              aaGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGc          <  2:2954009/149‑1 (MQ=17)
                              aaGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGc          <  2:2176309/149‑1 (MQ=17)
                              aaGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGc          <  2:2176295/149‑1 (MQ=17)
                              aaGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGc          <  1:1041613/149‑1 (MQ=17)
                               aGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCGGCGCATCGTGCa         >  1:400614/1‑149 (MQ=17)
                               aGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCa         >  2:2860888/1‑149 (MQ=17)
                               aGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCa         >  2:2860869/1‑149 (MQ=17)
                               aGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCa         <  2:406019/149‑1 (MQ=17)
                               aGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCa         <  2:406291/149‑1 (MQ=17)
                               aGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCa         >  1:1610337/1‑149 (MQ=17)
                                 cACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGc       <  1:3185691/149‑1 (MQ=17)
                                  aCTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCg      >  1:3648886/1‑149 (MQ=17)
                                   cTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCgg     >  1:1149444/1‑149 (MQ=17)
                                   cTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCgg     <  1:60551/149‑1 (MQ=17)
                                   cTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCgg     <  1:3561358/149‑1 (MQ=17)
                                   cTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCgg     <  1:2172961/149‑1 (MQ=17)
                                    tGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGt    <  2:3273833/149‑1 (MQ=17)
                                    tGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGt    >  2:895360/1‑149 (MQ=17)
                                    tGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGt    >  2:895396/1‑149 (MQ=17)
                                     gATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGtt   >  2:3519200/1‑149 (MQ=17)
                                      aTCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTtg  >  2:1732548/1‑149 (MQ=17)
                                      aTCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTtg  >  1:2639239/1‑149 (MQ=17)
                                      aTCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTtg  <  1:1688710/149‑1 (MQ=17)
                                      aTCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTtg  <  1:1688162/149‑1 (MQ=17)
                                         |                                                                                                                                                 
CCAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTG  >  CP009273/1565456‑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: