Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 710,187 A→T 100% S392C (AGT→TGT)  pgm → phosphoglucomutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273710,1870AT97.6% 529.4 / 2.5 163S392C (AGT→TGT) pgmphosphoglucomutase
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base T (73/86);  total (74/89)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
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 165 total aligned reads displayed.

CCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGAGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACGCT  >  CP009273/710044‑710331
                                                                                                                                               |                                                                                                                                                
  gTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAg                                                                                                                                            <  1:2780943/148‑1 (MQ=255)
         aaccGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGcttc                                                                                                                                     <  2:1229821/145‑1 (MQ=255)
         aaaCGCTGGTTTCATCGGCGATGATCGACCGTGTGGTCACCGACGTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGcttc                                                                                                                                     <  2:3336273/148‑1 (MQ=255)
           aCGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGcttcct                                                                                                                                   <  2:625225/148‑1 (MQ=255)
              cTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCttccttcc                                                                                                                                >  2:1069838/1‑148 (MQ=255)
              cTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCttccttcc                                                                                                                                >  2:580535/1‑148 (MQ=255)
                ggTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTg                                                                                                                              >  1:3093747/1‑148 (MQ=255)
                    tCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGtt                                                                                                                          <  2:859546/148‑1 (MQ=255)
                    tCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGtt                                                                                                                          >  1:2038423/1‑148 (MQ=255)
                      aTCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTc                                                                                                                        >  2:2510771/1‑148 (MQ=255)
                      aTCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTc                                                                                                                        >  1:778908/1‑148 (MQ=255)
                        cTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGa                                                                                                                      >  2:1659688/1‑148 (MQ=255)
                          gCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACg                                                                                                                    >  1:1012808/1‑148 (MQ=255)
                           cGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACgg                                                                                                                   >  1:3247760/1‑148 (MQ=255)
                             atgatCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCa                                                                                                                 <  1:2322798/148‑1 (MQ=255)
                                  cGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCATATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCg                                                                                                            <  2:1801545/148‑1 (MQ=255)
                                  cGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCg                                                                                                            >  2:1742938/1‑148 (MQ=255)
                                    aCCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTg                                                                                                          >  1:2058168/1‑148 (MQ=255)
                                     ccGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTgg                                                                                                         >  1:1191068/1‑148 (MQ=255)
                                     ccGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTgg                                                                                                         >  1:1377814/1‑148 (MQ=255)
                                      cGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGt                                                                                                        >  1:3005791/1‑148 (MQ=255)
                                            gTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCg                                                                                                  <  1:3186462/148‑1 (MQ=255)
                                            gTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCg                                                                                                  <  1:1364953/148‑1 (MQ=255)
                                             tCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGa                                                                                                 >  2:2224253/1‑148 (MQ=255)
                                               aaCGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACa                                                                                               <  1:210710/148‑1 (MQ=255)
                                                aCGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACaa                                                                                              <  2:279822/148‑1 (MQ=255)
                                                 cGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACaaa                                                                                             >  1:996952/1‑148 (MQ=255)
                                                  gACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAg                                                                                            <  2:130115/148‑1 (MQ=255)
                                                    cTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGAc                                                                                          >  1:503123/1‑148 (MQ=255)
                                                       gggCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGc                                                                                       <  1:510295/148‑1 (MQ=255)
                                                          ccGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGcatc                                                                                    <  1:286400/148‑1 (MQ=255)
                                                           cGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGcatca                                                                                   >  2:2653614/1‑148 (MQ=255)
                                                            gTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGcatcat                                                                                  >  1:2366046/1‑148 (MQ=255)
                                                             tAAGCTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTTGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGAGCGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGcatcatc                                                                                 <  2:1384858/148‑1 (MQ=255)
                                                             tAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGcatcatc                                                                                 >  1:2931692/1‑148 (MQ=255)
                                                             tAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGcatcatc                                                                                 <  2:3135566/148‑1 (MQ=255)
                                                             tAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGcatcatc                                                                                 >  1:1753815/1‑148 (MQ=255)
                                                               aaCTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGcatcatcat                                                                               <  2:460986/148‑1 (MQ=255)
                                                                aCTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCAtg                                                                              >  1:2679493/1‑148 (MQ=255)
                                                                aCTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCAtg                                                                              >  2:3284345/1‑148 (MQ=255)
                                                                 cTGGTAGACGTCCCGGTAGGTTTCAAATGGTTTGTCGTTGGTCTGTTCGACGGCAGCTTCGTCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCAtgt                                                                             <  2:2101474/148‑1 (MQ=255)
                                                                 cTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCAtgt                                                                             <  1:1875050/148‑1 (MQ=255)
                                                                 cTGGGAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCAtgt                                                                             <  2:2101986/148‑1 (MQ=255)
                                                                   ggTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTTTGAAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAATACGGCATCATCATctgt                                                                           >  1:360629/1‑148 (MQ=255)
                                                                   ggTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATgtgt                                                                           <  2:1617586/148‑1 (MQ=255)
                                                                   ggTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATgtgt                                                                           <  1:563643/148‑1 (MQ=255)
                                                                        aaGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTctgct                                                                      >  2:54888/1‑148 (MQ=255)
                                                                         aGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTctgctg                                                                     <  1:1083964/148‑1 (MQ=255)
                                                                         aGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTctgctg                                                                     <  1:3035847/148‑1 (MQ=255)
                                                                         aGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTctgctg                                                                     <  1:3036399/148‑1 (MQ=255)
                                                                          gTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTgg                                                                    <  2:3116898/148‑1 (MQ=255)
                                                                           tCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTggc                                                                   >  1:3283626/1‑148 (MQ=255)
                                                                            cccGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTggcg                                                                  >  1:3496495/1‑148 (MQ=255)
                                                                              cGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGGTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTggcggc                                                                >  2:1865626/1‑148 (MQ=255)
                                                                                 tAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGa                                                             >  2:3505183/1‑148 (MQ=255)
                                                                                 tAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGa                                                             >  2:1779965/1‑148 (MQ=255)
                                                                                  aGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGaa                                                            <  2:2042302/148‑1 (MQ=255)
                                                                                  aGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGaa                                                            <  1:3388526/148‑1 (MQ=255)
                                                                                  aGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCGTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGaa                                                            <  2:269671/148‑1 (MQ=255)
                                                                                   ggTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGaaa                                                           >  1:374540/1‑148 (MQ=255)
                                                                                   ggTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGaaa                                                           <  2:166111/148‑1 (MQ=255)
                                                                                    gTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAAt                                                          <  1:642753/148‑1 (MQ=255)
                                                                                     tttCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATc                                                         >  2:841906/1‑148 (MQ=255)
                                                                                        cAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCAcc                                                      >  2:3426584/1‑148 (MQ=255)
                                                                                         aaaTGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCg                                                     <  2:3015642/148‑1 (MQ=255)
                                                                                            tGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTg                                                  >  1:3211956/1‑148 (MQ=255)
                                                                                            tGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTg                                                  >  1:361847/1‑148 (MQ=255)
                                                                                                ttGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCAc                                              >  1:2580127/1‑148 (MQ=255)
                                                                                                 tGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCAcc                                             <  1:928040/148‑1 (MQ=255)
                                                                                                    cGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGt                                          <  2:1463555/148‑1 (MQ=255)
                                                                                                    cGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGt                                          <  2:2384411/148‑1 (MQ=255)
                                                                                                       tGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAg                                       >  1:2774752/1‑148 (MQ=255)
                                                                                                            tGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAAccc                                  >  1:2294195/1‑148 (MQ=255)
                                                                                                             gTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCg                                 >  1:1984438/1‑148 (MQ=255)
                                                                                                              ttCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGc                                <  1:2156466/148‑1 (MQ=255)
                                                                                                              ttCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGc                                <  1:2434319/148‑1 (MQ=255)
                                                                                                              ttCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGc                                <  1:885000/148‑1 (MQ=255)
                                                                                                              ttCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGc                                <  1:2156178/148‑1 (MQ=255)
                                                                                                                   cGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGaa                           >  1:576389/1‑148 (MQ=255)
                                                                                                                   cGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGaa                           <  2:2879561/148‑1 (MQ=255)
                                                                                                                     gCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAaca                         <  1:1704816/148‑1 (MQ=255)
                                                                                                                       aGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACt                       <  2:2877699/148‑1 (MQ=255)
                                                                                                                        gCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTa                      <  2:361847/148‑1 (MQ=255)
                                                                                                                        gCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTa                      <  1:3397911/148‑1 (MQ=255)
                                                                                                                        gCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTa                      <  1:3397888/148‑1 (MQ=255)
                                                                                                                        gCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTa                      <  2:3211956/148‑1 (MQ=255)
                                                                                                                         cTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTAc                     >  2:1740527/1‑148 (MQ=255)
                                                                                                                          ttCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACa                    <  1:2518495/148‑1 (MQ=255)
                                                                                                                           tCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACaa                   <  2:2371832/148‑1 (MQ=255)
                                                                                                                              gCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGa                <  2:2577346/148‑1 (MQ=255)
                                                                                                                              gCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGa                <  2:2576831/148‑1 (MQ=255)
                                                                                                                                tttGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAAc              >  2:1894856/1‑148 (MQ=255)
                                                                                                                                 ttGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACt             <  1:1199241/148‑1 (MQ=255)
                                                                                                                                 ttGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACt             >  2:955068/1‑148 (MQ=255)
                                                                                                                                   ggcggcGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTgg           <  1:2895104/148‑1 (MQ=255)
                                                                                                                                   ggcggcGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGCGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTgg           >  2:1136525/1‑148 (MQ=255)
                                                                                                                                          aaGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACg    >  1:1255831/1‑148 (MQ=255)
                                                                                                                                          aaGAGAGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACg    >  1:2590333/1‑148 (MQ=255)
                                                                                                                                           agagtgTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACGc   >  1:1767245/1‑148 (MQ=255)
                                                                                                                                            gagtgTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACGCt  >  2:1610537/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                
CCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGAGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACGCT  >  CP009273/710044‑710331

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: