Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,565,497 G→A 81.5% 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,4970GA81.5% 163.5 / 69.2 137S269L (TCG→TTG) gadBglutamate decarboxylase B, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (11/14);  new base A (63/48);  total (74/63)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 140 total aligned reads displayed.

CCCAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGG  >  CP009273/1565455‑1565646
                                          |                                                                                                                                                     
cccAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTa                                            <  1:1305190/150‑1 (MQ=11)
cccAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTa                                            <  1:3153904/150‑1 (MQ=11)
 ccAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTAt                                           <  1:1142329/150‑1 (MQ=31)
 ccAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTAt                                           <  1:670518/150‑1 (MQ=31)
  cAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATc                                          >  1:2596612/1‑150 (MQ=32)
  cAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATc                                          >  1:2596777/1‑150 (MQ=32)
   aGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCg                                         <  2:4425476/150‑1 (MQ=32)
   aGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCg                                         <  2:4164444/150‑1 (MQ=17)
   aGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCg                                         <  2:5384779/150‑1 (MQ=17)
   aGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCg                                         <  2:4164501/150‑1 (MQ=17)
      ggAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGcc                                      >  1:3392291/1‑150 (MQ=32)
       gAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCt                                     <  1:3534839/150‑1 (MQ=32)
       gAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCt                                     <  2:3545165/150‑1 (MQ=17)
       gAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCAGGCAGCGTCGATGTGCATGACGATGTCGATACCGGTATCGGCCt                                     >  1:367362/1‑150 (MQ=12)
        aGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTg                                    <  2:549929/150‑1 (MQ=17)
        aGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTg                                    <  2:549975/150‑1 (MQ=17)
         gCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTgg                                   <  2:3743291/150‑1 (MQ=17)
         gCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTgg                                   <  2:3713382/150‑1 (MQ=17)
          ccAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGa                                  >  2:2389503/1‑150 (MQ=32)
          ccAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGa                                  <  2:3214763/150‑1 (MQ=17)
           cAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGTAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGaa                                 <  2:3470690/150‑1 (MQ=17)
           cAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGaa                                 >  2:2718477/1‑150 (MQ=17)
           cAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGaa                                 >  1:4782813/1‑150 (MQ=17)
           cAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGAGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGaa                                 <  1:1289531/150‑1 (MQ=17)
            aGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAt                                <  2:5438514/150‑1 (MQ=17)
            aGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAt                                >  1:5358925/1‑150 (MQ=17)
            aGACCGAATGTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAt                                <  2:851872/150‑1 (MQ=17)
             gACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAtt                               >  2:4913669/1‑150 (MQ=17)
             gACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAtt                               <  1:2837341/150‑1 (MQ=17)
             gACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAtt                               >  2:2167030/1‑150 (MQ=17)
             gACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAtt                               >  2:2084554/1‑150 (MQ=17)
             gACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAtt                               >  2:1176180/1‑150 (MQ=17)
             gACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAtt                               >  1:1373491/1‑150 (MQ=17)
              aCCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAttt                              <  1:3826898/150‑1 (MQ=17)
              aCCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCAGGAAttt                              >  1:5527540/1‑150 (MQ=17)
               ccGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTa                             <  2:659491/150‑1 (MQ=17)
               ccGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTa                             <  2:659506/150‑1 (MQ=17)
               ccGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTa                             <  2:659528/150‑1 (MQ=17)
                cGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTAt                            >  1:1053842/1‑150 (MQ=32)
                cGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTAt                            >  2:3197707/1‑150 (MQ=17)
                cGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTAt                            >  1:3774851/1‑150 (MQ=17)
                    tttATGGCCTGAATCACTGATCTATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAg                        <  2:581729/150‑1 (MQ=17)
                    tttATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAg                        <  2:2855905/150‑1 (MQ=32)
                    tttATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAg                        <  2:581804/150‑1 (MQ=32)
                    tttATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAg                        <  1:1050215/150‑1 (MQ=17)
                    tttATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAg                        <  1:4484542/150‑1 (MQ=17)
                    tttATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAg                        >  1:4226300/1‑150 (MQ=17)
                      tATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAgcg                      >  2:277523/1‑150 (MQ=32)
                      tATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAgcg                      >  1:4903771/1‑150 (MQ=17)
                      tATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAgcg                      >  1:4358608/1‑150 (MQ=17)
                      tATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATAACCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAgcg                      >  1:1523560/1‑150 (MQ=17)
                       aTGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAgcgc                     >  1:292893/1‑150 (MQ=17)
                        tGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCa                    >  1:547778/1‑150 (MQ=17)
                         ggCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCAt                   >  1:3704722/1‑150 (MQ=17)
                         ggCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCAt                   >  2:4035316/1‑150 (MQ=17)
                         ggCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCAt                   <  2:5161517/150‑1 (MQ=17)
                          gCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATc                  >  2:4939825/1‑150 (MQ=32)
                          gCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATc                  <  2:348436/150‑1 (MQ=17)
                           ccTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCg                 <  2:1808470/150‑1 (MQ=32)
                           ccTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCg                 <  2:1808499/150‑1 (MQ=32)
                             tGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTg               >  2:4634953/1‑150 (MQ=17)
                             tGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTg               <  2:2598068/150‑1 (MQ=17)
                             tGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTg               >  2:4634976/1‑150 (MQ=17)
                             tGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTg               <  2:2317502/150‑1 (MQ=17)
                              gAAGCACTGATCAATTTCCCACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGc              <  1:4286666/150‑1 (MQ=17)
                              gAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGACTTTATCCAGCGCATCGTGc              <  1:3302696/150‑1 (MQ=17)
                              gAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGc              >  1:4974386/1‑150 (MQ=17)
                              gAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGc              <  1:4286965/150‑1 (MQ=17)
                              gAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGc              <  1:4655520/150‑1 (MQ=17)
                              gAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGc              >  2:2046531/1‑150 (MQ=17)
                              gAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGc              >  2:2874653/1‑150 (MQ=17)
                              gAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGc              <  1:4286849/150‑1 (MQ=17)
                              gAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGAGc              >  2:2375948/1‑150 (MQ=17)
                               aaGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCa             <  2:1900585/150‑1 (MQ=17)
                               aaGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCa             >  2:4887145/1‑150 (MQ=17)
                               aaGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCa             >  1:4496015/1‑150 (MQ=17)
                               aaGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCa             >  1:2563550/1‑150 (MQ=17)
                               aaGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTAACGGCCTGGAATGTATCCAGCGCATCGTGCa             >  2:4887194/1‑150 (MQ=12)
                                aGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAg            <  2:2697199/150‑1 (MQ=32)
                                aGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAg            >  1:250134/1‑150 (MQ=17)
                                 gCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGc           <  1:5426888/150‑1 (MQ=17)
                                  cACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCg          <  1:4800040/150‑1 (MQ=17)
                                  cACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCg          >  1:710599/1‑150 (MQ=17)
                                   aCTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACTGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCgg         >  2:4387403/1‑150 (MQ=17)
                                   aCTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCgg         >  1:5581642/1‑150 (MQ=17)
                                    cTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGt        >  1:3836795/1‑150 (MQ=17)
                                    cTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGt        >  2:3524524/1‑150 (MQ=17)
                                    cTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGtt       >  1:1415199/1‑150 (MQ=12)
                                    cTGATCAATTACACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGt        >  1:3836085/1‑150 (MQ=17)
                                     tGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGtt       <  1:356049/150‑1 (MQ=32)
                                     tGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGtt       >  1:5048949/1‑150 (MQ=17)
                                      gATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTtg      >  2:4711436/1‑150 (MQ=32)
                                      gATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTtg      >  1:5331438/1‑150 (MQ=17)
                                       aTCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTtgt     <  1:2021505/150‑1 (MQ=32)
                                       aTCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTtgt     <  1:3636899/150‑1 (MQ=17)
                                       aTCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTtgt     <  1:3760525/150‑1 (MQ=17)
                                         cGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGTCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTgg   >  2:3573533/1‑150 (MQ=17)
                                         cGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTgg   >  2:726450/1‑150 (MQ=32)
                                          gATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTggg  >  2:539483/1‑150 (MQ=32)
                                          gATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTggg  >  2:2568613/1‑150 (MQ=32)
                                          |                                                                                                                                                     
CCCAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGG  >  CP009273/1565455‑1565646

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: