Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,913,430 G→A 14.7% A215T (GCG→ACG)  spuI → glutamylpolyamine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,913,4300GA14.7% 244.0 / 30.8 102A215T (GCG→ACG) spuIglutamylpolyamine synthetase
Reads supporting (aligned to +/- strand):  ref base G (42/45);  new base A (7/8);  total (49/53)
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.06e-01

Only 100 of 104 total aligned reads displayed.

TCGCCAATTTCGGGCAAACGCCCGCAGTCGACCCAGGTGTACCTGATCGACGACCTCGACGAATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCGCCTACGACCATGAAATGGACACCACCTTCATGGCCAAGCCCTATCCGGGCCAGGCAG  >  NC_002947/5913283‑5913577
                                                                                                                                                    |                                                                                                                                                   
tCGCCAATTTCGGGCAAACGCCCGCAGTCGACCCAGGTGTACCTGATCGACGACCTCGACGAATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGgcg                                                                                                                                                   <  1:1843066/150‑1 (MQ=255)
tCGCCAATTTCGGGCAAACGCCCGCAGTCGACCCAGGTGTACCTGATCGACGACCTCGACGAATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGacg                                                                                                                                                   <  1:1563881/150‑1 (MQ=255)
  gCCAATTTCGGGCAAACGCCCGCAGTCGACCCAGGTGTACCTGATCGACGACCTCGACGAATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCa                                                                                                                                                 >  1:2003908/1‑150 (MQ=255)
   ccAATTTCGGGCAAACGCCCGCAGTCGACCCAGGTGTACCTGATCGACGACCTCGACGAATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAg                                                                                                                                                >  1:2241391/1‑150 (MQ=255)
    cAATTTCGGGCAAACGCCCGCAGTCGACCCAGGTGTACCTGATCGACGACCTCGACGAATATGCCGACTGCCTGCAGCCCACTATTCCACGCTGACTTA‑AATGAGGCATCTCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGt                                                                                                                                               >  1:284703/1‑150 (MQ=17)
         tCGGGCAAACGCCCGCAGTCGACCCAGGTGTACCTGATCGACGACCTCGACGAATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGaa                                                                                                                                          <  1:807968/150‑1 (MQ=255)
         tCGGGCAAACGCCCGCAGTCGACCCAGGTGTACCTGATCGACGACCTCGACGAATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGaa                                                                                                                                          <  2:2257153/150‑1 (MQ=255)
           gggCAAACGCCCGCAGTCGACCCAGGTGTACCTGATCGACGACCTCGACGAATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGt                                                                                                                                        <  1:1472092/150‑1 (MQ=255)
                   gCCCGCAGTCGACCCAGGTGTACCTGATCGACGACCTCGACGAATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGc                                                                                                                                <  2:2003908/150‑1 (MQ=255)
                     ccGCAGTCGACCCAGGTGTACCTGATCGACGACCTCGACGAATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGcac                                                                                                                              >  2:1045503/1‑150 (MQ=255)
                           tCGACCCAGGTGTACCTGATCGACGACCTCGACGAATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTg                                                                                                                        <  2:1321720/150‑1 (MQ=255)
                           tCGACCCAGGTGTACCTGATCGACGACCTCGACGAATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGACGCAGTTCGAAGTCAACCTGCACCACGTg                                                                                                                        >  1:1138427/1‑150 (MQ=255)
                            cGACCCAGGTGTACCGGATCGATGACCTCAACGAATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGTGGGTGGCATCCCGGCCGACGCCGTCGCCAATCAAATCCCCCCGGCGCAGTTCGAAGGCAACCTGCACCACGTgg                                                                                                                       <  1:1622650/150‑1 (MQ=255)
                             gACCCAGGTGTACCTGATCGACGACCTCGACGAATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGc                                                                                                                      >  2:1686699/1‑150 (MQ=255)
                               cccAGGTGTACCTGATCGACGACCTCGACGAATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCg                                                                                                                    >  1:336976/1‑150 (MQ=255)
                                 cAGGTGTACCTGATCGACGACCTCGACGAATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGAc                                                                                                                  <  2:1856256/150‑1 (MQ=255)
                                     tgtACCTGATCGACGACCTCGACGAATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGc                                                                                                              >  2:2216199/1‑150 (MQ=255)
                                         ccTGATCGACGACCTCGACGAATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCaa                                                                                                          >  2:561517/1‑150 (MQ=255)
                                          cTGATCGACGACCTCGACGAATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGA‑CCGCTCAAgg                                                                                                        >  2:1880068/1‑150 (MQ=255)
                                               cgacgaCCTCGACGAATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTg                                                                                                    >  2:1337573/1‑150 (MQ=255)
                                                 acgacCTCGAGGAATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGAc                                                                                                                  >  2:2274772/1‑134 (MQ=255)
                                                 acgacCTCGAGGAATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGAc                                                                                                                  <  1:2274772/134‑1 (MQ=255)
                                                     ccTCGACGAATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTa                                                                                              <  1:1686699/150‑1 (MQ=255)
                                                      cTCGACGAATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCgcttac                                                                                             <  2:656043/150‑6 (MQ=255)
                                                      cTCGACGAATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTAc                                                                                             <  1:1134748/150‑1 (MQ=255)
                                                       tCGACGAATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTAcg                                                                                            >  2:445692/1‑150 (MQ=255)
                                                        cgacgaATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTAcgc                                                                                           <  2:1958316/150‑1 (MQ=255)
                                                        cgacgaATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTAcgc                                                                                           >  2:1712701/1‑150 (MQ=255)
                                                          acgaATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCgg                                                                                         >  1:384167/1‑150 (MQ=255)
                                                                    cGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCaa                                                                               <  1:445692/150‑1 (MQ=255)
                                                                    cGACGGCCTGCAGGACA‑GCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCaa                                                                               <  1:1045503/150‑1 (MQ=255)
                                                                        tggctgcAGGACA‑TCCTCGAACGCGCGAAAGAACAACGCACCCCGGCCGAGGCCATCGTCAAGGAAAGCGCCCCGACGCAGTTCCAAGTCAACCTGCACCACGTGGCCGCCCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCgg                                                                           <  2:547615/147‑1 (MQ=25)
                                                                        tgcctgcAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCgg                                                                           <  1:561517/150‑1 (MQ=255)
                                                                        tgcctgcAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCgg                                                                           <  1:1712701/150‑1 (MQ=255)
                                                                         gcctgcAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGt                                                                          <  1:1337573/150‑1 (MQ=255)
                                                                          cctgcAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGcaccac                                                                                                                           >  1:1552432/1‑100 (MQ=255)
                                                                          cctgcAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGcaccac                                                                                                                           <  2:1552432/100‑1 (MQ=255)
                                                                          cctgcAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGtt                                                                         <  1:1009100/150‑1 (MQ=255)
                                                                          cctgcAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGtt                                                                         >  2:214402/1‑150 (MQ=255)
                                                                          cctgcAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGACGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGtt                                                                         <  1:2313489/150‑1 (MQ=255)
                                                                           ctgcAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTg                                                                        >  2:1830016/1‑150 (MQ=255)
                                                                           ctgcAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTg                                                                        >  2:1506542/1‑150 (MQ=255)
                                                                           ctgcAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTg                                                                        <  2:284703/150‑1 (MQ=255)
                                                                            tgcAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCt                                                                                                     >  2:437927/1‑120 (MQ=255)
                                                                            tgcAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCt                                                                                                     <  1:437927/120‑1 (MQ=255)
                                                                              cAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATc                                                                     >  1:2209703/1‑150 (MQ=255)
                                                                                ggACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCaa                                                                   <  1:1830016/150‑1 (MQ=255)
                                                                                 gACA‑TCCTCGAAGGCGCGAAATAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGaa                                                                                                                                          <  2:1566716/78‑1 (MQ=255)
                                                                                 gACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGaa                                                                                                                                          >  1:1566716/1‑78 (MQ=255)
                                                                                 gACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGGCGGCCCGCTCAAGGCCTGCTACTACGCGGTACTGCTCAAGCGGTTGATCAAg                                                                  <  1:1506542/150‑1 (MQ=255)
                                                                                 gACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCACGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAg                                                                  <  1:2293649/150‑1 (MQ=255)
                                                                                   cA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGaa                                                                <  2:384167/150‑1 (MQ=255)
                                                                                   cA‑TCCTCGAAGGCGCGAAAGAACAAGGCAGCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGaa                                                                >  1:167541/1‑150 (MQ=255)
                                                                                    a‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGACGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAAc                                                               >  1:31212/1‑150 (MQ=255)
                                                                                       ccTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACAt                                                             >  2:1086645/1‑150 (MQ=255)
                                                                                       ccTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCAGACCCGCTCAAGGCCTGCGACTACGCGGTACCGCTCAAGCGGTTGATCAAGAACAt                                                             >  2:1783529/1‑150 (MQ=255)
                                                                                         tCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGAc                                                                                                                  >  2:618378/1‑95 (MQ=255)
                                                                                         tCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGAc                                                                                                                  <  1:618378/95‑1 (MQ=255)
                                                                                         tCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGAt                                                                      >  2:115788/1‑139 (MQ=255)
                                                                                         tCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGAt                                                                      <  1:115788/139‑1 (MQ=255)
                                                                                         tCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGAGCCCCGACGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCg                                                           >  2:1070223/1‑150 (MQ=255)
                                                                                                cgcgAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCGCCTACGa                                                    <  1:1164295/150‑1 (MQ=255)
                                                                                                cgcgAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCGCCTACGa                                                    <  1:214402/150‑1 (MQ=255)
                                                                                                cgcgAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCACGGCCCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCGCCTACGa                                                    <  2:1772139/150‑1 (MQ=255)
                                                                                                 gcgAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGA‑CCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCGCCTACGAcc                                                  <  1:1880068/150‑1 (MQ=255)
                                                                                                   gAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATc                                                            >  2:1334286/1‑139 (MQ=255)
                                                                                                   gAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATc                                                            <  1:1334286/139‑1 (MQ=255)
                                                                                                   gAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCGCCTACGACCa                                                 >  1:2214799/1‑150 (MQ=255)
                                                                                                   gAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCGCCTACGACCa                                                 >  1:807720/1‑150 (MQ=255)
                                                                                                     aaGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCGCCTACGACCATg                                               >  2:2096047/1‑150 (MQ=255)
                                                                                                        aaCAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGACGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCGCCTACGACCATGaaa                                            <  1:1070223/150‑1 (MQ=255)
                                                                                                        aaCAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCACCGGCGCAGTTCGAAGTCAACCGACACCACGGGGCCGACCCGCTCAAGGCCTGCTACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCGCCTACGACCATGaaa                                            >  2:884035/1‑150 (MQ=21)
                                                                                                           aaGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGACGCAg                                                                                                                                                <  2:1362730/47‑1 (MQ=255)
                                                                                                           aaGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGACGCAg                                                                                                                                                >  1:1362730/1‑47 (MQ=255)
                                                                                                           aaGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGACGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCGCCTACGACCATGAAATgg                                         >  1:603268/1‑150 (MQ=255)
                                                                                                           aaGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGACGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCGCCTACGACCATGAAATgg                                         >  2:2284716/1‑150 (MQ=255)
                                                                                                             ggCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGACGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCGCCTACGACCATGAAATGGac                                       <  1:2284716/150‑1 (MQ=255)
                                                                                                                aTCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGa                                                                       <  1:1489751/115‑1 (MQ=255)
                                                                                                                aTCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGa                                                                       >  2:1489751/1‑115 (MQ=255)
                                                                                                                 tCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGc                                                                                                   >  2:588820/1‑86 (MQ=255)
                                                                                                                 tCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGc                                                                                                   <  1:588820/86‑1 (MQ=255)
                                                                                                                   ccagCCGACGGCATCGTCAAGGAAAGCGCCCCCGCGCAGTTCGAAGCCAACCTCCACCATGTGGACGACCCGGTCAGGTCCTGCGCCTATGCGGTACTGCTCGAGCGGTTGATCAAGAACATCGCCTACGACGATGAAATGGACaccacc                                 <  2:807720/147‑1 (MQ=17)
                                                                                                                    cGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCGCCTACGACCATGAAATGGACACCACCt                                <  1:95301/150‑1 (MQ=255)
                                                                                                                        cGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCGCCTACGACCATGAAATGGACACCACCTTCAt                            <  2:2241391/150‑1 (MQ=255)
                                                                                                                        cGACGCCATCGTCAAGGAAAGCGCCCCGACGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCGCCTACGACCATGAAATGGACACAACCTTcct                            <  2:603268/150‑3 (MQ=255)
                                                                                                                         gACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCGCCTACGACCATGAAATGGACACCACCTTCATg                           >  2:613861/1‑150 (MQ=255)
                                                                                                                         gACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCGCCTACGACCATGAAATGGACACCACCTTCATg                           >  2:1025638/1‑150 (MQ=255)
                                                                                                                              cATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTCCGCGGTACTGCTCAAGc                                                                             <  2:1231443/95‑1 (MQ=255)
                                                                                                                              cATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGc                                                                             >  1:1231443/1‑95 (MQ=255)
                                                                                                                               aTCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCGCCTACGACCATGAAATGGACACCACCTTCATGGCCAAg                     >  2:2379985/1‑150 (MQ=255)
                                                                                                                                 cGTCAAGGAAAGCGCCCCGACGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCGCCTACGACCATGAAATGGACACCACCTTCATGGCCAAGcc                   <  2:1138427/150‑1 (MQ=255)
                                                                                                                                   tCAAGGAAAGCGCCCCGACGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCGCCTACGACCATGAAATGGACACCACCTTCATg                           >  1:2247916/1‑140 (MQ=255)
                                                                                                                                   tCAAGGAAAGCGCCCCGACGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCGCCTACGACCATGAAATGGACACCACCTTCATg                           <  2:2247916/140‑1 (MQ=255)
                                                                                                                                     aaGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCGCCTACGACCATGAAATGGACACCACCTTCATGGCCAAGCCCTAt               >  2:495680/1‑150 (MQ=255)
                                                                                                                                            gcgcCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCGCCTACGACCATGAAATGGACACCACCTTCATGGCCAAGCCCTATCCGGGcc        <  2:2209703/150‑1 (MQ=255)
                                                                                                                                              gcCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCGCCTACGACCATGAAATGGACACCACCTTCATGGCCAAGCCCTATCCGGGCcag      <  1:2379985/150‑1 (MQ=255)
                                                                                                                                              gcCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCGCCTACGACCATGAAATGGACACCACCTTCATGGCCAAGCCCTATCCGGGCcag      <  1:2096047/150‑1 (MQ=255)
                                                                                                                                                cccGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGtt                                                                         <  2:1053277/81‑1 (MQ=255)
                                                                                                                                                cccGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGtt                                                                         >  1:1053277/1‑81 (MQ=255)
                                                                                                                                                  cGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCGCCTACGACCATGAAATGGACACCACCTTCATGGCCAAGCCCTATCCGGGCcaggcag  >  2:770791/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCGCCAATTTCGGGCAAACGCCCGCAGTCGACCCAGGTGTACCTGATCGACGACCTCGACGAATATGCCGACTGCCTGCAGGACA‑TCCTCGAAGGCGCGAAAGAACAAGGCATCCCGGCCGACGCCATCGTCAAGGAAAGCGCCCCGGCGCAGTTCGAAGTCAACCTGCACCACGTGGCCGACCCGCTCAAGGCCTGCGACTACGCGGTACTGCTCAAGCGGTTGATCAAGAACATCGCCTACGACCATGAAATGGACACCACCTTCATGGCCAAGCCCTATCCGGGCCAGGCAG  >  NC_002947/5913283‑5913577

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