Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,565,497 G→A 58.1% 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,4970GA58.1% ‑2.7 / 98.7 74S269L (TCG→TTG) gadBglutamate decarboxylase B, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (13/18);  new base A (28/15);  total (41/33)
Fisher's exact test for biased strand distribution p-value = 6.01e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

CCCAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGG  >  CP009273/1565455‑1565646
                                          |                                                                                                                                                     
cccAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTa                                            <  2:4866043/150‑1 (MQ=11)
 ccAGCGGAGCCAGACCGTATTTATGGCCTGAAGCACTGATCAATTTCATACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTAt                                           <  1:3002854/150‑1 (MQ=9)
 ccAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTAt                                           <  2:2339749/150‑1 (MQ=31)
  cAGCGGAGCCAGACCGAATTTCTGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATc                                          <  2:439563/150‑1 (MQ=17)
   aGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATTTCGATGTCGATACCGGTATCg                                         <  1:3192280/150‑1 (MQ=17)
   aGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCg                                         >  2:546255/1‑150 (MQ=17)
      ggAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGcc                                      >  1:3576380/1‑150 (MQ=32)
      ggAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGcc                                      >  1:80589/1‑150 (MQ=32)
      ggAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGACCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGcc                                      >  1:80614/1‑150 (MQ=17)
      ggAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGcc                                      >  1:123862/1‑150 (MQ=17)
      ggAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGcc                                      >  1:718437/1‑150 (MQ=17)
      ggAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGcc                                      >  1:2758518/1‑150 (MQ=17)
      ggAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGcc                                      >  1:2758741/1‑150 (MQ=17)
      ggAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGcc                                      >  1:2759028/1‑150 (MQ=17)
       gAGCCAGACCGAATTTATGGCCTGACGCACTGATCGATTTCACCCGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCt                                     <  1:3866043/150‑1 (MQ=17)
       gAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCt                                     >  1:2625248/1‑150 (MQ=17)
       gAGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCt                                     >  1:2625082/1‑150 (MQ=17)
        aGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTg                                    >  1:752033/1‑150 (MQ=32)
        aGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTg                                    >  1:752330/1‑150 (MQ=32)
        aGCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTg                                    <  2:4231906/150‑1 (MQ=17)
         gCCAGACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTgg                                   <  1:4421760/150‑1 (MQ=17)
          ccAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGa                                  >  2:3275333/1‑150 (MQ=32)
             gACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAtt                               >  1:2851167/1‑150 (MQ=17)
             gACCGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAAtt                               <  2:2834228/150‑1 (MQ=17)
               ccGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTa                             >  2:279522/1‑150 (MQ=32)
               ccGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTa                             <  2:3691351/150‑1 (MQ=32)
               ccGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTa                             >  2:2918492/1‑150 (MQ=32)
               ccGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTa                             <  2:3692791/150‑1 (MQ=32)
               ccGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCTTCGATGTGCATGACGATGTCGATACAGGTATCGGCCTGGAATTTa                             >  1:3504403/1‑150 (MQ=255)
               ccGAATTTATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTa                             >  2:1161346/1‑150 (MQ=17)
                cGAATTTAGGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGAGATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTAt                            <  2:3555104/150‑1 (MQ=17)
                cGAATGTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTAt                            <  2:3555124/150‑1 (MQ=17)
                 gAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATc                           <  1:3416122/150‑1 (MQ=32)
                  aaTTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATcc                          <  2:1610904/150‑1 (MQ=32)
                    tttATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAg                        <  1:3023711/150‑1 (MQ=32)
                     ttATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAgc                       <  2:3852343/150‑1 (MQ=17)
                     ttATGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAgc                       <  2:3852368/150‑1 (MQ=17)
                      tATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAgcg                      <  1:316911/150‑1 (MQ=32)
                        tGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCa                    >  1:738053/1‑150 (MQ=17)
                        tGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCa                    >  1:737963/1‑150 (MQ=17)
                        tGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCa                    <  2:3583711/150‑1 (MQ=17)
                        tGGCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCa                    <  2:3583717/150‑1 (MQ=17)
                          gCCTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATc                  <  1:3750123/150‑1 (MQ=17)
                           ccTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCg                 >  1:4257770/1‑150 (MQ=32)
                            cTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGt                >  2:5074515/1‑150 (MQ=32)
                            cTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGt                >  2:3297364/1‑150 (MQ=17)
                            cTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGt                >  2:670561/1‑150 (MQ=17)
                            cTGAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGt                >  2:670737/1‑150 (MQ=17)
                             tGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTg               <  2:2485344/150‑1 (MQ=32)
                              gAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGACTTTATCCAGCGCATCGTGc              >  2:3620595/1‑150 (MQ=17)
                              gAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGc              >  2:2482777/1‑150 (MQ=17)
                              gAAGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGc              >  2:4742215/1‑150 (MQ=17)
                               aaGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCa             >  2:3226190/1‑150 (MQ=17)
                               aaGCACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCa             <  1:5150572/150‑1 (MQ=17)
                                aGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAg            <  2:2178036/150‑1 (MQ=32)
                                  cACTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCg          <  1:1339087/150‑1 (MQ=17)
                                   aCTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCgg         >  2:2765070/1‑150 (MQ=17)
                                   aCTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCgg         >  1:4578517/1‑150 (MQ=17)
                                   aCTGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCgg         >  1:759288/1‑150 (MQ=17)
                                    cTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGt        <  1:4965944/150‑1 (MQ=32)
                                    cTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGt        <  1:2373054/150‑1 (MQ=32)
                                    cTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGt        <  1:4965922/150‑1 (MQ=32)
                                    cTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGt        <  1:424920/150‑1 (MQ=32)
                                     tGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGtt       >  2:1710520/1‑150 (MQ=32)
                                     tGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGtt       <  1:1457268/150‑1 (MQ=32)
                                     tGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGtt       >  2:1710572/1‑150 (MQ=32)
                                     tGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGtt       <  1:4375568/150‑1 (MQ=17)
                                     tGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGtt       >  2:2724082/1‑150 (MQ=17)
                                     tGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGtt       >  1:2403940/1‑150 (MQ=17)
                                     tGATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGtt       <  1:4374555/150‑1 (MQ=17)
                                      gATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTtg      >  1:5063334/1‑150 (MQ=32)
                                      gATCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTtg      >  1:2997281/1‑150 (MQ=17)
                                       aTCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCCGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTtgt     >  2:3691977/1‑150 (MQ=17)
                                       aTCAATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTtgt     >  2:3691962/1‑150 (MQ=17)
                                          gATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTggg  >  2:1827037/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: