Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 2,910,161 T→C intergenic (+65/‑594) gudD → / → RPPX_12965 glucarate dehydratase/GABA permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099742,910,1610TC85.9% 979.6 / 96.6 394intergenic (+65/‑594)gudD/RPPX_12965glucarate dehydratase/GABA permease
Reads supporting (aligned to +/- strand):  ref base T (12/43);  new base C (126/210);  total (138/256)
Fisher's exact test for biased strand distribution p-value = 3.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 414 total aligned reads displayed.

GGTGCGCGGGATGACAGCGTGGCGATGCAGTTCCTGATTCCGGGTTGGAGCTTCGATAACAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCCCCAGTCATTCACCGCCGCACCTGCCACCGCATTCATCAATCGGCACCATCGCCACGCACCCTTGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCC  >  CP009974/2910013‑2910311
                                                                                                                                                    |                                                                                                                                                      
 gTGCGCGGGATGACAGCGTGGCGATGCAGTTCCTGATTCCGGGTTGGAGCTTCGATAACAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACcgtcgtc                                                                                                                                                     >  2:6222763/1‑150 (MQ=255)
  ggcgcGGGATGACAGCGTGGCGATGCAGTTCCTGATTCCGGGTTGGAGCTTCGATAACAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAg                                                                                                                                                    <  1:5892755/150‑1 (MQ=255)
     gcgGGATGACAGCGTGGCGATGCAGTTCCTGATTCCGGGTTGGAGCTTCGATAACAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAg                                                                                                                                                 >  2:2714943/1‑151 (MQ=255)
      cgGGATGACAGCGTGGCGATGCAGTTCCTGATTCCGGGTTGGAGCTTCGATAACAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGc                                                                                                                                                <  2:5450691/151‑1 (MQ=255)
      cgGGATGACAGCGTGGCGATGCAGTTCCTGATTCCGGGTTGGAGCTTCGATAACAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGc                                                                                                                                                >  1:467903/1‑151 (MQ=255)
       gggATGACAGCGTGGCGATGCATTTCCTGATTCCGGGTTGGAGCTTCGATAACAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCGCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCg                                                                                                                                               <  2:4500116/151‑1 (MQ=255)
       gggATGACAGCGTGGCGATGCAGTTCCTGATTCCGGGTTGGAGCTTCGATAACAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCg                                                                                                                                               >  2:3959598/1‑151 (MQ=255)
         gATGACAGCGTGGCGATGCAGTTCCTGATTCCGGGTTGGAGCTTCGATAACAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCggg                                                                                                                                             >  2:3296365/1‑151 (MQ=255)
         gATGACAGCGTGGCGATGCAGTTCCTGATTCCGGGTTGGAGCTTCGATAACAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCggg                                                                                                                                             >  1:813439/1‑151 (MQ=255)
           tGACAGCGTGGCGATGCAGTTCCTGATTCCGGGTTGGAGCTTCGATAACAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCTGGCAAGCCCGCTCCTACATCAACCTGTGGTGTCGCCCCACACCGGCGTAGGAGCggggt                                                                                                                                           >  2:2100953/1‑149 (MQ=255)
           tGACAGCGTGGCGATGCAGTTCCTGATTCCGGGTTGGAGCTTCGATAACAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGtt                                                                                                                                           >  2:6100574/1‑151 (MQ=255)
             aCAGCGTGGCGATGCAGTTCCTGATTCCGGGTTGGAGCTTCGATAACAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCa                                                                                                                                         <  1:3870626/151‑1 (MQ=255)
               aGCGTGGCGATGCAGTTCCTGATTCCGGGTTGGAGCTTCGATAACAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCAcc                                                                                                                                       >  2:282048/1‑151 (MQ=255)
               aGCGTGGCGATGCAGTTCCTGATTCCGGGTTGGAGCTTCGATAACAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCGCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCAcc                                                                                                                                       <  2:3771125/151‑1 (MQ=255)
                 ggTGGCGATGCAGTTCCTGATTCCGGGTTGGAGCTTCGATAACAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCcg                                                                                                                                     <  2:1961682/150‑1 (MQ=255)
                  gTGGCGATGCAGTTCCTGATTCCGGGTTGGAGCTTCGATAACAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCcgc                                                                                                                                    <  1:1027250/151‑1 (MQ=255)
                    ggCGATGCAGTTCCTGATTCCGGGTTGGAGCTTCGATAACAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGa                                                                                                                                  >  2:2955040/1‑151 (MQ=255)
                    ggCGATGCAGTTCCTGATTCCGGGTTGGAGCTTCGATAACAAACGGCCTTGCCGGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCGACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGa                                                                                                                                  <  1:4163814/151‑1 (MQ=255)
                      cGATGCAGTTCCTGATTCCGGGTTGGAGCTTCGATAACAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAAt                                                                                                                                >  1:1715260/1‑151 (MQ=255)
                       gATGCAGTTCCTGATTCCGGGTTGGAGCTTCGATAACAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGAATCGCCACGCACCGTCGTAGAAGCCGGTTCACCCGCGaaaa                                                                                                                               >  1:4977088/1‑149 (MQ=25)
                           cAGTCCCGGATTCCGGGTTGGACCTTCGATAACAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGcc                                                                                                                           <  2:53873/151‑1 (MQ=255)
                            aGTTCCTGATTCCGGGTTGGAGCTTCGATAACAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCg                                                                                                                          >  2:5554234/1‑151 (MQ=255)
                               tCCTGATTCCGGGTTAGAGCTTCGATAACAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGCGCGGGTTCACCCGCGAATAGGCCGGaa                                                                                                                       >  1:1188924/1‑150 (MQ=255)
                                    aTTCCGGGTTGGAGCTTCGATAACAACCGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTgcc                                                                                                                  <  2:3210735/151‑1 (MQ=255)
                                       ccGGGTTGGAGCTTCGATAACAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTgccgcc                                                                                                               >  1:2556934/1‑151 (MQ=255)
                                            ttGGAGCTTCGATAACAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATa                                                                                                          <  1:6265313/151‑1 (MQ=255)
                                               gAGCTTCGATAACAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATAcct                                                                                                       <  1:6073056/151‑2 (MQ=255)
                                                    tCGATAACAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGt                                                                                                  >  2:5941519/1‑151 (MQ=255)
                                                         aaCAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTc                                                                                             <  1:576474/151‑1 (MQ=255)
                                                         aaCAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTc                                                                                             >  1:6147828/1‑151 (MQ=255)
                                                         aaCAAACGGCCTTGCCTGGTGCGTTGAGGCTACGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGAGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTc                                                                                             <  2:5104063/151‑1 (MQ=255)
                                                              aCGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTg                                                                                                                    >  1:6387554/1‑123 (MQ=255)
                                                              aCGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCAccgc                                                                                        <  1:536168/151‑1 (MQ=255)
                                                                gtccTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCCCCAGTCATTCAccgccg                                                                                      >  1:5626203/3‑151 (MQ=255)
                                                                ggCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCAccgccg                                                                                      <  2:3099458/151‑1 (MQ=255)
                                                                  ccTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCa                                                                                    <  2:4441871/151‑1 (MQ=255)
                                                                   cTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCAc                                                                                   <  1:3760808/151‑1 (MQ=255)
                                                                       ccTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCCCCCGc                                                                               >  2:2612229/1‑151 (MQ=255)
                                                                         tGGGGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCa                                                                             <  1:2745955/151‑1 (MQ=255)
                                                                           gTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCCCCGCCGCACCCGCCAcc                                                                           <  2:3543915/151‑1 (MQ=255)
                                                                                  gAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCCCCCGCCACCGCCTtca                                                                    >  1:3945746/1‑151 (MQ=25)
                                                                                  gAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATtca                                                                    <  1:3175092/151‑1 (MQ=255)
                                                                                      cTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCaa                                                                <  1:4035586/151‑1 (MQ=255)
                                                                                         cGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCg                                                             <  2:3659512/151‑1 (MQ=255)
                                                                                            gCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCa                                                          >  1:770271/1‑151 (MQ=255)
                                                                                             cgtctTCGCGGGCAAGCCCGCTCCTACAGCGACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCAc                                                         <  2:1577713/149‑1 (MQ=21)
                                                                                              ctctTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTg                                                                                                                    <  1:3132707/91‑1 (MQ=255)
                                                                                                ctTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCAccgt                                                      >  1:1008650/1‑149 (MQ=255)
                                                                                                 ttCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCAccgtc                                                     >  1:3452426/1‑148 (MQ=255)
                                                                                                 ttCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCACCGTc                                                     <  1:514455/151‑1 (MQ=255)
                                                                                                 ttCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCACCGTc                                                     <  1:3367203/151‑1 (MQ=255)
                                                                                                 ttCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCCCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCACCGTc                                                     <  2:4712783/151‑1 (MQ=255)
                                                                                                    gcgGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATAcc                                                                                                        <  2:6501940/97‑1 (MQ=255)
                                                                                                    gcgGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATc                                                              >  1:218849/1‑139 (MQ=255)
                                                                                                    ccgGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCACCGTCGcc                                                  <  1:3267727/150‑1 (MQ=255)
                                                                                                     cgGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCACCGTCGCCa                                                 <  1:425046/151‑1 (MQ=25)
                                                                                                          aaGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCCCCAGTCATTCACCGCCGCACCTGCCACCGCATTCATCAATCGGCACCATCGCCACGCAc                                            <  2:426221/151‑1 (MQ=255)
                                                                                                             cccGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATGCCCCCAGTCATTCACCGTCGCACCCGCCACCGCATtc                                                                     <  1:3029898/123‑1 (MQ=17)
                                                                                                             cccGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCCCCAGTCATTCACCGTCGCACCCGCCACCGCATtc                                                                     >  2:3029898/1‑123 (MQ=255)
                                                                                                              ccGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGGGCAGCCGCGCATGGGCCGGACCTGCCGCCCAATACCTGGAGGCATTCACCGCCGCACCCGCCACCGCATTCGCCAACCAGCACGGTCGCCACGCAccgtg                                        >  2:849546/1‑148 (MQ=18)
                                                                                                               cGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCACCGTCGCCACGTACCGTTg                                       <  2:2331327/151‑1 (MQ=21)
                                                                                                               cGCCCCTACATCAGCCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCACCGTCGCCACGTACCGTTg                                       <  2:4105397/151‑1 (MQ=21)
                                                                                                                   ccTACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCACCGTCGCCACGTACCGTTGTAgg                                   <  1:475616/151‑1 (MQ=21)
                                                                                                                   cattcATCGACCGGCACTGTCGCCCCGCACCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCCCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCACCGTCGCCACGCCCCCTTGTAgg                                   <  2:4831179/147‑1 (MQ=17)
                                                                                                                     tACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCACCGTCGCCACGTACCGTTGTAGGAg                                 <  2:5292543/151‑1 (MQ=21)
                                                                                                                     tACATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCACCGTCGCCACGTACCGTTGTAGGAg                                 <  1:6422254/151‑1 (MQ=21)
                                                                                                                      aCATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCACCGTCGCCACGTACCGTTGTAGGAGc                                >  1:5527017/1‑151 (MQ=21)
                                                                                                                       cATCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCACCGTCGCCACGTACCGTTGTAGGAGCg                               <  2:5021109/151‑1 (MQ=21)
                                                                                                                        aTCAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGCCACCGTCGCCACGTCCCGTTGTAGGAGCgg                              >  1:1595338/1‑151 (MQ=21)
                                                                                                                          cAACCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCACCGTCGCCACGTACCGTTGTAGGAGCGGGt                            <  1:5664366/151‑1 (MQ=21)
                                                                                                                           caCCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCACCGTCGCCACGTACCGTTGTAGGAGCGGGtt                           <  1:2911059/150‑1 (MQ=21)
                                                                                                                           aaCCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCACCGTCGCCACGTACCGTTGTAGGAGCGGGtt                           <  2:5616447/151‑1 (MQ=21)
                                                                                                                           aaCCGGTGGTGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCACCGTCGCCACGCCCCCTTGTAGGAGc                                >  1:43871/1‑146 (MQ=21)
                                                                                                                               ggtggtGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCTGCCACCGCATTCATCAATCGGCACCATCGCCACGCACCCTTGTAGGAGCgg                              >  2:1360484/1‑144 (MQ=255)
                                                                                                                               ggtggtGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCACCGTCGCCACGTACCGTTGTAGGAGCGGGTTCAcc                       <  2:119162/151‑1 (MQ=21)
                                                                                                                               ggtggtGTCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCCCCAGTCATTCACCGCCGCACCTGCCACCGCATTCATCAATCGGCACCATCGCCACGCACCCTTGTAGGAGCgg                              <  1:1360484/144‑1 (MQ=255)
                                                                                                                                gcggcgtCGCCACGCCCCCTTGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCCCCAGTCATTCACCGCCGCACCTGCCACCGCATTCATCAATCGGCACCATCGCCACGCACCCTTGTAGGAGCGGGTTCAccc                      <  2:6387679/146‑1 (MQ=18)
                                                                                                                                 cggtgCCGCCCCGCACCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCGATACCCCCAGTCATTCACCGCCGCACCTGCCACCGCATTCATCAATCGGCACCGTCGCCACGCACCCTTGTAGGAGCGGGTTCACCcg                     <  2:3453609/150‑1 (MQ=255)
                                                                                                                                 caccgtCGCCACGCCCCCTTGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCTCAATACCCCCAGTCATTCACCGCCGCACCTGCCCCCGCATTCACCAATCGGCACCCTCGTCACGCACCCTCGGAGGCGCCGGTGCACCcg                     >  1:189740/5‑151 (MQ=14)
                                                                                                                                  cccgtCGCCCCGCCCCCTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCCCCAGTCATTCACCGCCGCACCTGCCACCGCATTCATCGATCGGCACCATCGCCACGCACCCTTGTAGGAGCGGGTTCACCcgc                    <  1:4913362/148‑1 (MQ=21)
                                                                                                                                   gtgtCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCACCGTCGCCACGTACCGTTGTAGGAGCGGGTTCACCcgcg                   <  1:2421496/151‑1 (MQ=21)
                                                                                                                                    cgtCGCCACGCCCCCTTGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCCCCAGTCATTCACCGCCGCACCTGCCACCGCATTCATCAATCGGCACCGTCGCCACGCACCCTTGTAGGAGCGGGTTCACCCGCGa                  <  1:1925323/150‑1 (MQ=21)
                                                                                                                                     gtCGCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCACCGTCGCCACGTACCGTTGTAGGAGCGGGTTCACCCGCGaa                 <  1:6481424/151‑1 (MQ=21)
                                                                                                                                      tCGCCCCGCACCGTCGCAGGAGCGGGTTCCCCCGCGAATAGGCCGGCCCTGCCGCCCCACACCCCCAGTCATTCACCGCCGCACCCGCCCCCGCATTCATCAATCGGCACCATCGCCACGCACCCTTGTAGGAGCGGGTTCACCCGCGAAt                <  2:6346961/151‑1 (MQ=37)
                                                                                                                                        gccacgcacCCTTGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATGCCCCCCGTCGTTAACCGCCCCCCCTGCCACCGCATCCATCATCCGGCCCCGTCGCCACGCACCCTTGTAGGAGCGGGTTCACCCGCGAATAg              <  2:1766987/145‑1 (MQ=18)
                                                                                                                                        gCCCCACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCCCCAGTCATTCAcc                                                                                          <  2:746307/75‑1 (MQ=255)
                                                                                                                                         ccccACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCCCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCACCGTCGCCACGCCCCCTTGTAGGAGCGGGTTCACCCGCGAATAgg             <  1:900631/151‑1 (MQ=255)
                                                                                                                                          cccACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGccc                                                                                                              >  1:5062473/1‑53 (MQ=255)
                                                                                                                                          cccACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCACCGTCGCCACGTACCGTTGTAGGAGCGGGTTCACCCGCGAATAGGc            <  1:4538721/151‑1 (MQ=21)
                                                                                                                                          cacgcacCCTTGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCCCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCACCGTc                                                     <  1:2067603/106‑1 (MQ=11)
                                                                                                                                           ccggacCGTGGTGGGAGCGGGTTCACCCGCGAATAGGCCGAACCTGCCGCCCAATACCCCCGGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCACCGTCGCCACGCCCCCTTGTAGGAGCGGGTTCACCCGCGAATAGGcc           <  2:3351764/147‑1 (MQ=21)
                                                                                                                                           ccACACCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCACCGTCGCCACGTACCGTTGTAGGAGCGGGTTCACCCGCGAATAGGcc           <  2:4636919/151‑1 (MQ=21)
                                                                                                                                            cgtacCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCCCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCACCGTCGCCACGCCCCCTTGTAGGAGCGGGTTCACCCGCGAATAGGCCg          <  1:4070780/148‑1 (MQ=255)
                                                                                                                                             gtacCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCCCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCAccgtc                                                     >  2:5778901/3‑104 (MQ=255)
                                                                                                                                             gcacCCTTGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCCCCAGTCATTCACCGCCGCACCCGCCACCGCATTCGTCAATCGGCACCGTCACCACGCACCCTTGTAGGAGCGGGGTCACCCGCGAATAGGCCgg         >  1:3654660/2‑151 (MQ=18)
                                                                                                                                             acacCGTCGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCCCCAGTCATTCACCGCCGCACCCGCCa                                                                             <  2:4038022/83‑1 (MQ=255)
                                                                                                                                               acCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCCCCAGTc                                                                                                 <  2:14404/61‑1 (MQ=35)
                                                                                                                                               acCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCCCCAGTc                                                                                                 >  1:14404/1‑61 (MQ=35)
                                                                                                                                                 cgtcgtAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCTCCAGTCATTCACCGCCGCACCCGCCACCGCATTCATCAATCGGCACCGTCGCCACGCCCCCTTGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCt     <  2:1013763/151‑1 (MQ=25)
                                                                                                                                                  cttGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCCCCAGTCATTCACCGCCGCACCTGCCACCGCATTCATCAATCGGCACCATCGCCACGCACCCTTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGACCTg    >  1:2658917/2‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GGTGCGCGGGATGACAGCGTGGCGATGCAGTTCCTGATTCCGGGTTGGAGCTTCGATAACAAACGGCCTTGCCTGGTGCGTTGAGGCTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACATCAACCGGTGGTGTCGCCCCACACCGTTGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCCGCCCAATACCCCCAGTCATTCACCGCCGCACCTGCCACCGCATTCATCAATCGGCACCATCGCCACGCACCCTTGTAGGAGCGGGTTCACCCGCGAATAGGCCGGACCTGCC  >  CP009974/2910013‑2910311

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