Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,555,156 G→T 100% A171E (GCG→GAG)  yjiM ← putative 2‑hydroxyglutaryl‑CoA dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,555,1560GT99.7% 1086.2 / ‑4.5 323A171E (GCG→GAG) yjiMputative 2‑hydroxyglutaryl‑CoA dehydratase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (177/145);  total (178/145)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 323 total aligned reads displayed.

TCCAGTTCATTGATCAACGCCTCTTTATCGAACCGGAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCGCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGCAATTGCATCACATGAACAGGC  >  CP009273/4555009‑4555302
                                                                                                                                                   |                                                                                                                                                  
    gTTCATTGATCAACGCCTCTTTATCGAACCGGAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATgg                                                                                                                                              <  2:652395/150‑1 (MQ=255)
     ttCATTGATCAACGCCTCTTTATCGAACCGGAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGc                                                                                                                                             >  2:1040914/1‑150 (MQ=255)
       cATTGATCAACGCCTCTTTATCGAACCGGAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCAt                                                                                                                                           >  2:3721181/1‑150 (MQ=255)
         ttGATCAACGCCTCTTTATCGAACCGGAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATcg                                                                                                                                         <  1:281386/150‑1 (MQ=255)
           gATCAACGCCTCTTTATCGAACCGGAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATcgcg                                                                                                                                       <  2:75331/150‑1 (MQ=255)
              cAACGCCTCTTTATCGAACCGGAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCGCAATGGCATCGCGCag                                                                                                                                    >  1:4065642/1‑150 (MQ=255)
               aaCGCCTCTTTATCGAACCGGAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCaga                                                                                                                                   >  1:5077649/1‑150 (MQ=255)
                 cGCCTCTTTATCGAACCGGAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGc                                                                                                                                 >  1:3117886/1‑150 (MQ=255)
                 cGCCTCTTTATCGAACCGGAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGc                                                                                                                                 >  1:1631742/1‑150 (MQ=255)
                 cGCCTCTTTATCGAACCGGAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGc                                                                                                                                 <  1:3772322/150‑1 (MQ=255)
                  gCCTCTTTATCGAACCGGAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCa                                                                                                                                >  2:3327981/1‑150 (MQ=255)
                   ccTCTTTATCGAACCGGAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCAt                                                                                                                               >  2:2201058/1‑150 (MQ=255)
                   ccTCTTTATCGAACCGGAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCAt                                                                                                                               >  2:4876505/1‑150 (MQ=255)
                     tctTTATCGAACCGGAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCt                                                                                                                             <  1:3939123/150‑1 (MQ=255)
                       tttATCGAACCGGAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTc                                                                                                                           <  1:4838275/150‑1 (MQ=255)
                       tttATCGAACCGGAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTc                                                                                                                           <  1:3454069/150‑1 (MQ=255)
                         tATCGAACCGGAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGc                                                                                                                         <  2:2958843/150‑1 (MQ=255)
                          aTCGAACCGGAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCt                                                                                                                        >  1:2900612/1‑150 (MQ=255)
                             gAACCGGAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAAt                                                                                                                     >  2:3654904/1‑150 (MQ=255)
                             gAACCGGAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAAt                                                                                                                     >  2:2617215/1‑150 (MQ=255)
                             gAACCGGAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAAt                                                                                                                     <  1:3721643/150‑1 (MQ=255)
                             gAACCGGAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAAt                                                                                                                     >  1:1625782/1‑150 (MQ=255)
                             gAACCGGAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAAt                                                                                                                     <  2:1883304/150‑1 (MQ=255)
                               aCCGGAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAAtct                                                                                                                   <  1:4737595/150‑1 (MQ=255)
                               aCCGGAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAAtct                                                                                                                   >  1:552305/1‑150 (MQ=255)
                                 cGGAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCg                                                                                                                 >  1:247299/1‑150 (MQ=255)
                                  ggAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGt                                                                                                                >  2:2898643/1‑150 (MQ=255)
                                    aaGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGc                                                                                                              >  1:3042129/1‑150 (MQ=255)
                                      ggTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGccc                                                                                                            >  1:1790869/1‑150 (MQ=255)
                                       gTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCa                                                                                                           >  1:4006097/1‑150 (MQ=255)
                                          gcgcCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCaaaa                                                                                                        >  2:4531132/1‑150 (MQ=255)
                                            gcCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACg                                                                                                      <  2:4800423/150‑1 (MQ=255)
                                              cGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGtt                                                                                                    >  2:1235690/1‑150 (MQ=255)
                                                tAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGttct                                                                                                  <  2:127893/150‑1 (MQ=255)
                                                 aaaCCACTTTCAGACTGTCGCTGCCGCTAAGCGCCGGAGGATTTAACGGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTAGTGCCCAAAACGttctt                                                                                                 <  2:1017926/150‑1 (MQ=255)
                                                   aCCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTtcttct                                                                                               <  1:3212779/150‑1 (MQ=255)
                                                   aCCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTtcttct                                                                                               >  1:928211/1‑150 (MQ=255)
                                                    ccACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTa                                                                                              <  2:4123970/150‑1 (MQ=255)
                                                    ccACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTa                                                                                              <  1:3549736/150‑1 (MQ=255)
                                                    ccACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTa                                                                                              <  1:1587477/150‑1 (MQ=255)
                                                    ccACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAAGGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTa                                                                                              <  2:2372749/150‑1 (MQ=255)
                                                       cTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCg                                                                                           <  2:3501371/150‑1 (MQ=255)
                                                        tttCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGt                                                                                          <  2:1187788/150‑1 (MQ=255)
                                                        tttCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGt                                                                                          <  1:2704044/150‑1 (MQ=255)
                                                          tCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGttt                                                                                        >  2:2836066/1‑150 (MQ=255)
                                                                    gCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAgcg                                                                              <  1:3507453/150‑1 (MQ=255)
                                                                      tGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCgca                                                                            <  2:2443300/150‑1 (MQ=255)
                                                                         cGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCTCGACGTTCGCGGTTTTTCAGCTCAATGGAATCGCGCAGAGCATCTTCGCTAATCTAGTGCCCAAAACGTTATTCTACCGTTTTTTGCAAGCgaagca                                                                         >  1:1021876/1‑150 (MQ=255)
                                                                         cGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCgcagca                                                                         >  1:4175951/1‑150 (MQ=255)
                                                                             aaGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCAtctc                                                                     >  2:4374344/1‑150 (MQ=255)
                                                                              aGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCg                                                                    >  1:1143819/1‑150 (MQ=255)
                                                                                cgcCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGc                                                                  <  2:827181/150‑1 (MQ=255)
                                                                                     gaggaTTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTcc                                                             <  1:1751452/150‑1 (MQ=255)
                                                                                      aggaTTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCa                                                            <  1:3571080/150‑1 (MQ=255)
                                                                                            tAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAacgc                                                      >  2:4589840/1‑150 (MQ=255)
                                                                                            tAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAacgc                                                      >  2:3933641/1‑150 (MQ=255)
                                                                                            tAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAacgc                                                      >  2:571070/1‑150 (MQ=255)
                                                                                            tAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGATCGCGGTTATTCAGCTCAATGTCATCGAGCAGAGCATCTTCGCGAATCTCGTGCCCAATACTTTCTTCTACCGAATTATGCAAGCGCAGCAGCTCCTCTTTCCATAacgc                                                      >  2:1817877/1‑150 (MQ=255)
                                                                                               cTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAacgcacg                                                   <  1:2489522/150‑1 (MQ=255)
                                                                                               cTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAacgcacg                                                   <  1:2560069/150‑1 (MQ=255)
                                                                                               cTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTACTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAacgcacg                                                   >  2:1590582/1‑150 (MQ=255)
                                                                                                tGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAacgcacgc                                                  >  1:3381383/1‑150 (MQ=255)
                                                                                                     aaGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGc                                             <  2:2801313/150‑1 (MQ=255)
                                                                                                      aGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCa                                            <  1:1999224/150‑1 (MQ=255)
                                                                                                        atgatAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATc                                          >  2:2695465/1‑150 (MQ=255)
                                                                                                         tgatACAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCg                                         <  2:1726067/150‑1 (MQ=255)
                                                                                                          gatAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGt                                        <  1:2211255/150‑1 (MQ=255)
                                                                                                          gatAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGt                                        >  2:4942400/1‑150 (MQ=255)
                                                                                                          gatAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGt                                        <  1:3289059/150‑1 (MQ=255)
                                                                                                             aaaaaTTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCt                                     >  2:308820/1‑150 (MQ=255)
                                                                                                              aaaaTTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATCACGCACGCGAGGCATCGTCCtt                                    >  1:2634235/1‑150 (MQ=255)
                                                                                                                   tAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGc                               >  2:3625923/1‑150 (MQ=255)
                                                                                                                    aGCCAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCt                              >  1:4616611/1‑150 (MQ=255)
                                                                                                                       cAGTGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGtt                           <  1:2739445/150‑1 (MQ=255)
                                                                                                                          tGCGCGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTggg                        <  2:4018819/150‑1 (MQ=255)
                                                                                                                           gcgcGACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGc                       <  1:868264/150‑1 (MQ=255)
                                                                                                                            cgcgACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGCa                      <  1:2898643/150‑1 (MQ=255)
                                                                                                                            cgcgACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGCa                      <  1:126171/150‑1 (MQ=255)
                                                                                                                             gcgACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGCaa                     >  1:241944/1‑150 (MQ=255)
                                                                                                                               gACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGCAAtt                   <  1:1050886/150‑1 (MQ=255)
                                                                                                                               gACGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGCAAtt                   >  1:2859693/1‑150 (MQ=255)
                                                                                                                                aCGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGCAATTg                  <  2:1790869/150‑1 (MQ=255)
                                                                                                                                aCGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGCAATTg                  <  1:64100/150‑1 (MQ=255)
                                                                                                                                aCGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGCAATTg                  >  1:3264016/1‑150 (MQ=255)
                                                                                                                                 cGTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGCAATTGc                 >  1:1904084/1‑150 (MQ=255)
                                                                                                                                  gTTCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGCAATTGCa                >  2:5120060/1‑150 (MQ=255)
                                                                                                                                   ttCGCGGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGCAATTGCAt               >  2:877919/1‑150 (MQ=255)
                                                                                                                                     ctcgGTTTTTCAGCTCAATGGCATAGCGCAGAGCATCTTCGATAATCTCTTGCCAAAAACGTTCTTCTACCGTTTTTTTCAAGCGCAGCATCTCGGCTTTCCATAACGCACTCGAGGCATCGTCCTTAACGCTGTTGGGCAATTtcatca             <  1:1499762/148‑1 (MQ=255)
                                                                                                                                     cgcgGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGCAATTGCATca             <  1:451275/150‑1 (MQ=255)
                                                                                                                                     cgcgGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGCAATTGCATca             <  1:2285686/150‑1 (MQ=255)
                                                                                                                                       tggTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACACGAGGCATCGTCCTTAACGCTGTTGGGCAATTGCATcaca           >  1:2519817/2‑150 (MQ=255)
                                                                                                                                       cgGTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGCAATTGCATcaca           <  2:3015072/150‑1 (MQ=255)
                                                                                                                                         gTTTTTCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGCAATTGCATCACATg         >  2:4314825/1‑150 (MQ=255)
                                                                                                                                          tttttCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGCAATTGCATCACATGa        <  2:1980856/150‑1 (MQ=255)
                                                                                                                                           ttttCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGCAATTGCATCACATGaa       <  1:902927/150‑1 (MQ=255)
                                                                                                                                            tttCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGCAATTGCATCACATGAAc      <  2:501153/150‑1 (MQ=255)
                                                                                                                                              tCAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGCAATTGCATCACATGAACAg    >  1:3151534/1‑150 (MQ=255)
                                                                                                                                               cAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGCAATTGCATCACATGAACAgg   >  1:4977252/1‑150 (MQ=255)
                                                                                                                                               cAGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGCAAGTGCATCACATGAACAgg   >  2:1065578/1‑150 (MQ=255)
                                                                                                                                                aGCTCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGCAATTGCATCACATGAACAGGc  >  1:3237456/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCCAGTTCATTGATCAACGCCTCTTTATCGAACCGGAAGGTTGCGCCGTAAACCACTTTCAGAATGTCGCTGCCGCTAAGCGCCGGAGGATTTAACTGCCCAAGATGATAAAAATTAGCCAGTGCGCGACGTTCGCGGTTTTTCAGCGCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGCAATTGCATCACATGAACAGGC  >  CP009273/4555009‑4555302

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: